IL 


NUOVO CIMENTO 


ORGANO DELLA SOCIETA ITALIANA DI FISICA 
SOTTO GLI AUSPICI DEL CONSIGLIO NAZIONALE DELLE RICERCHE 


Von, HAN Serie decima 1° Agosto 1960 


On the Inelastic Scattering by Deformed Nuclei. 


Z. JANKOVIG 


Faculty of Science, University of Zagreb and Institute « Rudjer Boskovié » - Zagreb 


(ricevuto il 21 Marzo 1960) 


Summary. A more general treatment of the nuclear inelastic scattering 
on even-even target nuclei is given particularly by introducing other col- 
lective target states besides rotational ones, the deformed diffuse nuclear 
and spin-orbit potentials and the Coulomb potential for the deformed 
target nucleus. In the first order approximation explicit expressions for 
the differential cross-section (18) and for the total cross-section (19) for 
the transition 0-2 are deduced. 


In the last few years, the collective model description of the target nucleus 
was also applied in the attempts to explain the inelastic scattering of nuclear 
projectiles on nuclei with the excitation of rotational states and a direct inter- 
action between the target nucleus and the projectile has been assumed (1). 
We try to generalize the above mentioned problem of inelastic scattering and 
particularly by introducing also other collective target states (*), the deformed 


(1) D. M. Brink: Proc. Phys. Soc., 68, 994 (1955); S. HAYAKAWA and §. YOSHIDA: 
Progr. Theor. Phys., 14, 1 (1955). 

(2) S. YosHIpa: Proc. Phys. Soc., 69, 668 (1956). 

(3) D. M. Case, L. WiLers and A. R. Epmonps: Phys. Rev., 110, 1080 (1958). 

(4) Z. JANKOVIG: Nuovo Cimento, 14, 1174 (1959). Another paper, with more 
details, to be published. soon. 


19 - Il Nuovo Cimento. 
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smooth tail nuclear and spin-orbit potentials and the Coulomb potential of 
the deformed target nucleus for charged projectiles. 

Let the projectile a and the target nucleus X be described, in the center 
of mass system, by the sets of parameters & and & respectively. Then a Hamil- 
tonian of the form 


(1) H eta) se H(a) ir He È) 


describes the scattering of projectiles a on the targets X. Here, the first part 
H(é) is simply the kinetic energy operator of the projectile. 

The second term H(x) is the target nucleus Hamiltonian and we assume 
that the complete set of its orthonormalized eigenfunctions is at our disposal 


(2) (H(c) — E,)B,(x) = 0, fora Dax) dt, = dn. 


For the collective model even-even target nucleus wave functions we have 
the general form (4): 


(3) D, (x) = > Feng(B) Jxny(?) Dix (0%) 9) K (even) < I. 


We quote the following special cases: a) for the rotational case, 6 and y 
remain fixed and the sum (3) contains only one normalized Dix(d) function; 
b) for harmonic f stabilized vibrations the f wave function does not depend 
on K, unstabilized » vibrations not having any n, quantum; €) for f and 
y Stabilized (an) harmonie vibrations (with the small or vanishing equilibrium 
values fo, yo) the sum (3) reduces to one K term. 

The third term in the Hamiltonian (1) represents the interaction between 
the projectile and the target nucleus and is composed of three terms originating 
from nuclear forces, Coulomb force and spin-orbit coupling, all being dependent 
on both sets of parameters « and é: 


(4) Hilo €) = Vala, €) + Vela, E) + Vila, é)U-s. 


For the collective model the sharp deformed surface has been originally 
introduced by 


(5) k=k, (1 oF > a,Ys(9", g')) ’ Ry = TVA. 


a, being known functions of the deformation parameters B and y, &, y 
meaning the polar angles in the principal axes system. The nuclear forces 
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potential as well as the spin-orbit potential are functions of the distance 
r— R and one can expand them in the following way 


: . = 1 df(r—R) Seen 
Se i n tie | Pat, oh 


v 


where 7, is a constant (may be depending on target atomic number and pro- 
jectile’s energy) and À (and R,) has to be redefined as a mean value regarding 
the diffuse potential V. Thus, in the expression (6), we have a diffuse and 
deformed potential expanded in a power series of deformation parameters a,. 
Allowing a change of these parameters with time, the whole deformed and 
diffuse edge target, represented by the potential (6), executes vibrations and 
rotations. 

The optical model suggests a complex nuclear potential, in order to take 
eare of the target’s absorption of the projectiles, which leads in the most 
simple case to put 


(7) Vo=Vot two, 


where in general we could assume different r dependences both for real and 
imaginary parts, each part being then expanded like (6). 

Thus, by choosing a smooth complex potential (e.g. of the Woods-Saxon 
type) we have in expression (6) its series expansion, the zero order term being 
of spherical symmetry. 

For the square well potential we see from (6) that approximately, the first 
order term is represented by a spherical delta surface function which has been 
applied in (*) and (?). 

By assuming a spatial distribution of protons the Coulomb potential can 
also be expanded in a power series of deformation parameters a,. The expan- 
sion problem being here more difficult, we assume as a good approximation 
the uniform charge distribution in the deformed target nucleus with the sharp 
surface (5), in which case the expansion can be made easily. For further ap- 
plication we quote the first order approximation: 


47 RS STRO 7 Ne 
Vola, 6) = el | (+) 2 V2", 9), r> Ro 


i : 45 id r Ul 
Vola, €) = 2 Roo È — (| | 0 a r<Ro, 


where again the zero order terms are the known expressions for the Coulomb 
potential of a uniform spherically distributed charge. 
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To solve the scattering problem we expand the stationary solution of the 
Schrôdinger equation 


(9) (H — E)P(a, £) = 0, 


with H given by expression (1) and Æ representing the total energy of the 
projectile plus the target nucleus, into a series of the complete set of target 
wave functions 


(10) È D,(x)un(€) . 


Inserting the expression (10) in the Schrédinger equation (9) we find that 
the projectile’s wave functions have to satisfy this system of linear coupled 
second order differential equations 


he 


27 


(11) né) + (E — E,,) ZH mlt CEA 


m being the projectile’s reduced mass, 7 the set of quantum numbers specifying | 
completely the target quantum state and H7,, the following interaction Hamil- | 
tonian matrix elements 


(12) Hux(é) forma int(%, È) D,(2) dar La 


By the known procedure (°) the differential cross-section for a transition | 
î > is expressed in the following way: | 


(13) dors Q) = |r (UE) — UAC) do.) l'homme AQ - 


The system (11) contains in principle infinite number of differential equa- 
tions and cannot be solved without approximations with regard first to the 
number of equations (i.e. the number of target states taken into account) and 
second to the order of approximation in calculating the matrix elements (12). 
The most simple case represents the two states case where the spinless pro- 
jectile induces the target transition from the i to the f state, all other tran- 
sitions neglected, and only the lowest order target deformation taken into 


(5) J. M. Buarr and L. C. Brrpennarn: Rev. Mod. Phys., 24, 261 (1953); see 
also ref. (6). I 


2492 


ON THE INELASTIC SCATTERING BY DEFORMED NUCLEI 


bo 
co 
Qt 


account when calculating the matrix elements. Under these assumptions the 
system (11) reduces to the following two equations: 


|P Au(6) + (BB, — V(r) —V PU) wile) = 0, 


| 2m 
Lt) h2 
| Dr Au,(é) + (E-E;— Vr) — VO) u,(E) = Hfi(E) (E), 


where £ means the projectile’s spatial co-ordinate’s and the matrix element (12) 
calculated with the first order target deformation a, terms. For the more 
specified transition: i (fundamental state I=0, u — K — 0), f (first excited 
state J = 2, u; K=0) of an even-even target nucleus, one finds for the matrix 
element with wave functions (3): 


(15) H},(E) = MOVIE, 9), 


where the constant 0 takes on different values according to the different spe- 
cifications of the wave functions: 


yO =e 
a Lg = — 
WG 
b) c 36 h 
(16) je 
VT(@) 13) Bo 
¢) CCI) Ne Ty =1. 


I representing an integral of the corresponding f and g wave functions. Hach 
coefficient 0 can be numerically expressed when the collective model, i.e. the 
corresponding collective Hamiltonian, is completely specified. 

After performing a series expansion in spherical harmonics for both fun- 
ctions u,(7, 0, ~) and u,(r, Ÿ, 9) 


as: 


(17) u(r, Ÿ, 0) = S Vn(2t + 1) — u,(kr)Y7(0, 9), 
1=0 


(77 


one gets from (14) the radial wave equations: homogeneous for w,,, and in- 
homogeneous u,,. There is a known procedure (°) to find the asymptotic form 
of the solution w;,’? (and then by (17) also for u\) needed for the calculation 
of the differential cross-section (13) do%,(Q). 


(6) N. F. Morr and H. 8. W. Massey: The Theory of Atomic Collisions (Oxford, 
1952), p. 108. 
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A summation over the final state magnetic quantum number w gives the 
differential cross-section for the transition à (1 = 0) to f (I= 2): 


| 
il | 


2m ,\? ; PR 
8 dr = C avt. 
(18) do;,(9) (7 | 642k, ke 1, Du | 


(21, + 1)(24 + 1)(21, + 1264 DE, lux | 
(LU, 2:00) C(UI 2:00) Ol, L:00) CH, L:00)W (LULL £2) P,(cos 9) dQ 


and the integration gives the total cross-section 


2: totti 
(19) Cn = & 0) ICE D (24 + 1)(20; + 1)! Ti, [2 OU 2:00)? . 
SIT UTS Lt 


Here 7, , means the integral | 
Ft = 


œ 


(20) Lors, fr (kr) Ua, (kur) (V(r) + Ve (r)) dr, 


0 


both dashed % functions meaning the stationary solutions of the elastic scat- 
tering problem (homogeneous differential equations in (14)) vanishing at the 
origin and behaving as «7 — mu, at the infinity, «being the known asymp- 
totic forms exp [+ i(kr — mn 2kr — l(x/2))] for charged particles and with 
n= for neutral particles. 7, is the factor describing nuclear and Coulomb 
scattering. Clebsch-Gordan coefficients take care of the conservation of the 
total angular momentum, only even or odd J values can combine because 
of the conservation of parity. 

We see in this simple example, under assumptions which are essentially 
equivalent to the first order perturbation calculation, how the diffuseness, 
deformability and charge distribution of the target nucleus enter through the 
coefficient (20) into the differential and total cross-sections (18) and (19). The 
nature of the states è and f (k; and k, assumed to be the same) different under 
assumptions a), b) and c) influences only the absolute values of the differential 
and total cross-sections through the different coefficients C? given in (16) and 
has no influence on the angular distribution of the inelastically scattered pro- 
jectiles. 


This work has been done during the author’s stay at the Centre d'Etudes 
Nucléaires de Saclay. The author wishes to express his appreciation of the 
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arrangements which made possible his stay at the Centre d’Etudes Nucléaires 
de Saclay and is very much indebted to Prof. J. Turron for the hospitality 
at the Section des Réactions Nucléaires & Moyenne Energie. A grant from 
the Institute « Rudjer Boëkovié », Zagreb, is gratefullb acknowledged. 


‘ 


RIASSUNTO 


Si espone un trattamento più generale dello scattering nucleare anelastico su nuclei 
targhetta pari-pari, particolarmente introducendo altri stati collettivi della targhetta 
oltre a quelli rotazionali, i potenziali nucleari diffusi deformati e di spin-orbita e il 
potenziale coulombiano del nucleo deformato della targhetta. Con approssimazione di 
primo ordine si deducono espressioni esplicite per la sezione trasversale differenziale (18) 
e per la sezione trasversale totale (19) per la transizione 0 — 2, 


(*) Traduzione a cura della Redazione. 
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A Method for Trapping Muons in Magnetic Fields, 


and Its Application to a Redetermination of the EDM of the Muon. 


G. CHARPAK (*), L. M. LEDERMAN (*), J. C. SENS and A. ZICHICHI 
CERN - Geneva 


(ricevuto il 4 Aprile 1960) 


Summary. — A method for trapping muons in magnetic fields is described 
together with its application to higher accuracy determination of the 
electric dipole moment of the muon. The value found is: 


EDM < e(5 +5) :10-17 em 


consistent with time reversal invariance. 


1. — Introduction. 


The availability of beams of polarized muons with asymmetric decay has 
stimulated many groups to consider the extension to the case of the muon of 
the methods used by LovrsELL, Prop and CRANE (1), to measure the anom- 
alous magnetic moment of the electron. 

In this method, one can measure directly (g — 2) by storing the particles 
for some time in a magnetic field and then measuring the relative precession 
angle between the spin and the momentum, which serves as a reference vector. 
As in the electron experiments, the prime requirement is to succeed in injecting 
the muons into a magnetic field so that they circulate on essentially periodic 
orbits, and to trap them in this field for as large a number of « turns », de. 
orbit periods, as possible. The low fluxes, high momenta and large extensions 


(*) On leave from Centre National de la Recherche Scientifique, Paris. 

(**) Ford Foundation Fellow, on leave from Columbia University. Also J. S. Gug- 
genheim Fellow, 1958-1959. 

(1) W. K. H. Louisezz, R. W. Prop and H. R. CRANE: Phys. Rev., 94, 7 (1954). 
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in position and momentum space (low density in phase space) which the pres- 
ently available muon beams exhibit in comparison with conventional electron 
beams present, however, many new difficulties in meeting this requirement. 
On the other hand, there is one advantage of the muons—it is not necessary 
to analyze the spin polarization by scattering, since the asymmetric electron 
decay reveals the spin deviation. 

Several practical methods for overcoming these difficulties have been in- 
vestigated experimentally and/or theoretically over the past year at CERN. 
We have now succeeded in trapping muons of 85 MeV/c momentum for 28 turns, 
and our results indicate that minor modifications in the method used should 
enable one to achieve storage for several hundreds of turns. The main pur- 
pose of this note is to describe this method. 

As a first application, we have used equipment that was primarily devel- 
oped for trapping studies, for a redetermination of the electric dipole moment. 
(EDM) of the muon. It is known (?) that elementary particles can possess 
an EDM only if they have some interaction not invariant under time-reversal. 
In this sense, a search for an EDM constitutes a test of the validity of this 
invariance principle. Like previous workers (*) we have found a value for 
the muon EDM that is zero within the experimental errors. These errors were, 
however, considerably reduced with respect to those of refs. (*)—mainly in 
virtue of the much longer time which the muons spend in the field in our 
arrangement. The second purpose of this note is to discuss our EDM deter- 
mination and its results. 


2. — Injection and trapping. 


The principle of our method, to be described with reference to Fig. 1, is. 
to inject a muon beam into the median (X — Y) plane of a flat magnet gap 
(say, along the Y-axis), and to contain it therein by reducing the momentum 
by means of a moderator M, to a value p, the latter being suited to the mean 
vertical field B, and the lateral extension of the gap. The magnetic field 
configuration is so chosen as to produce a drift of the muon orbits along the 
X-axis away from M, and to focus them in the median plane. The trapping 
efficiency obtainable is of course limited by the Liouville theorem. 


(2) L. Lanpau: Nucl. Phys., 3, 127 (1957); I. B. ZeLpovicH: Zurn. Éksp. Teor. 
Piz., 33, 1488 (1957). Also T. D. Lee and C. N. Yane: B.N.L., 443 (1957). 

(3) D. Bertey, R. L. GARWIN, G. Gipat and L. M. LEDERMAN: Phys. Rev. Lett., 1, 
144 (1958); Annual International Conference on High Energy Physics, CERN (1958), 
p. 332; D. BerLEY, G. Gipaz and L. M. LEDERMAN: Kiev Conf., 1959 (to be publi- 
shed). See also R. L. GARWIN and L. M. LEDERMAN: Nuovo Cimento, 11, 776 (1959). 
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Our field had the form: 


(1) B, = By(1+ ay — by?) 


in the median plane (the exit point of M being taken as the origin). With | 


ar <1 and br? € 1, the muons emerging from M will move on nearly circular 


D ZAN 


M Là d 
2 | 
3 
ul 
:9 Scintillator C 
| Lucite 
Section E-E 
Fig. 1. — Experimental layout. A) magnetic channel; 8B) scintillator; D) lucite block; | 


M) moderator; 1, 2, 3) scintillators. 


orbits of radius r. As is well known, a linear gradient alone leads to a drift 
of these orbits in the X-direction by an amount (« step-size »): 


(2) s= ar (grad B,/B,> = xr°a/turn, 


where <> denotes average over one orbit loop. This drift enables some of 
the muons to clear M after their first turn, whereupon they proceed on a tro- 
choidal orbit. The linear gradient produces also a feeble vertical focusing, 
partly because the field index (n) is larger in the weak field region (vertically 
focusing) than in the strong field region (vertically defocusing) and partly 
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because of the «strong focusing » effects of even equal alternating n-values. 
An accurate analysis (*) gives for the vertical focusing wavelength: 


Oe 


This focusing is, however, very weak (A,/2ar>>1), because one wants r/s>>1 
in order to store as large as possible a number of turns in a magnet of given 
size. With r/s —16 for example, one has one vertical oscillation per 38 turns. 
Thus, as we could verify experimentally, the vertical divergence, due to mul- 
tiple scattering of the muons emerging from M, leads to poor storage as long 
as one provides only a linear gradient. In addition, it is very difficult in such 
a field to keep the magnetic median plane near the centre of the gap. Should 
the median plane touch the poles, all particles are lost. 

As was pointed out by GARWIN, one can improve the vertical focusing 
while maintaining a given large value of r/s by the addition of the by? term. 
In fact, for the field given by (eq. (1)), one has: 


Ay eer Di 
(b + 1.74a2)-} NL 


9 
a) 


while the step-size is still given by eq. (2) for orbits which start out parallel 
to the Y-axis. By taking, e.g. b = 50a?, one can, while maintaining the same 
r/s above (for such orbits), improve the focusing to 1 oscillation/7 turns. Thus 
the intensity of stored muons is increased by a factor 38/7~5 by the ad- 
dition of the quadratic term to the field. The by? term serves also the useful 
function of fixing more firmly the magnetic median plane in the centre of 
the magnet gap. An additional consideration is the location of the edge of 
the stability bands for orbits in a field of the form (1). The result quoted for 
/,/2% is of course correct only for small b and a. For large b and a the ver- 
tical focusing first improves and then becomes unstable. These questions are 
treated in reference (*). It is, however, not permissible to choose b arbitrarily 
large, because this would lead to a spread in step-sizes and hence in storage 
times—even if a mono-energetic beam could be made to emerge from M. 
Orbits emerging at an angle ® with respect to the Y-axis would have a 
step-size: 


(3) s(D) = ar?(a — 2br) 


(4) R. L. Garwin: Numerical calculations of the stability bands and solutions of a 
Hill differential equation, CERN internal report (October 1959); W. K. H. PANOFSKY: 
Orbits in the linear magnet, CERN internal report (October 1959). 
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as can be seen by introducing eq. (1) into eq. (2). A spread in @ is in practice 
unavoidable, and is reflected in a spread in s according to eq. (3). For a fixed 
@ one has, on the other hand, an additional momentum dependence as com- 
pared to the case with b=0. Both of these effects depend on b, and its 
magnitude must be chosen so as to maximize the number of particles stored 
for a given number of turns. 

The field of eq. (1) (with B, =13.1kG, a=10-* cm, b — 65-107 cm*) was 
obtained by modifying an existing magnet of gap volume 83x 53x10 cm’. 
The linear gradient was produced by inserting wedge-shaped iron plates behind 

the pole caps. Additional shims 
served to provide the quadratic 
term and to correct for field 
irregularities. Iron bars of suit- 
able profile, mounted along all 
the edges of the gap, served to 
increase its useful width and 
length by cancelling the field 
fall-off originally present. Out 


to a distance of 2.5 em from 
the edges, the so resulting field 
3 CATA, agreed to better than 0.2% 
Lae ee ) with the prescribed dependence 

A (eq. (1). 
In the actual experiments, 
amt, wt beam of momentum 
Ale 150 MeV/c was focused, by 
means of a pair of quadrupole 
= lenses (not indicated in Fig. 1), 
CR NO ns Duc ARS through an iron channel, which 
EE, bei DS PORT EN. served to eliminate the fringing 
200 160 120 


e orata field of the magnet, onto the 
Fig. 2. - Time distribution of the stored muons as SOU EDs M. The moderator 
displayed on the pulse height analyser, Distance Consisted of a Be block of 
between moderator M and counter O,7,=13em. 14.7 g/em? thickness (Be was 


C detects muons having completed an integer used because of its low Z and 


UGLY QE ne high density), of a few graphite 


plates for minor range adjust- 
ments, and of the three scintillators of counters 1, 2, 3 constituting an injec- 
tion telescope. The curved shape of this assembly, indicated in Fig. 1, enabled 
the muons to clear it after their first turn. M reduced the initial momentum 
to pa 85 MeV/c. A polyethylene bag filled with helium gas at atmospheric 
pressure served to reduce ionization loss and scattering below that due to air. 
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Even so, the muon population can be shown to decrease exponentially due to 
multiple scattering with a decay length of ~ 42 turns. 

The storage was studied by placing a scintillator C on the X-axis at a 
distance X, from M, and measuring the distribution of the times of arrival 
of the muons striking ©. This distribution is determined by means of a 
time-to-pulse-height converter (T.C.) system, which had (123) and C for inputs, 
and was displayed on a many- 
channel pulse-height analyser 

A ‘og 2 ar ar -  4| Muons 
EE Figs. 2 and 3a are ty per second 
pical examples of two distri- per channel. 

È - : ; . (arbitrary units) 
butions obtained in this way, 


PS È A 3F 
differing in the position at 
which counter 0 was placed. 
For X,< 2r (*) one de- oe b) 


tects in 0 mostly muons 
which have completed an 
integral number of turns (*), 42 


Fig. 3. — a) Time distribution 300 
of the stored muons displayed 
on the pluse-height analyser. 
Distance between moderator M 
and counter 0, æ,—66 cm. C 4F 
detects muons having complet- 
ed n+4 turns, where » is an 
integer. Group A: particles with 3 
mean turn number <n>) = 1, 
reaching © after excursion 
through the fringing field. 
Groups B,, B,, Bs: particles 2F 
with <n),=11.5; <n),=16.5; 
<n>, = 19.5. The arrows indi- 
cate the single channel analyser ;| 
window limits. - b) Calculated 


time distribution. The area of | ele 

È Ò È ; a 4 eran eni st 1 
this distribution is normalized 30% x re ‘pla SIAE Lion a 
to the experimentally obtained 200 160 120 40 40 (o) 


area from channel 7 to 20. —+——'me of arrival (mus) 


(*) From here on we indicate with a bar (b) all quantities (as for instance, radius, 
step-size, etc.) corresponding to the momentum p. 

(**) It is possible for muons to cross the X-axis at X > X, because counter O was 
built with an « air light-pipe ». The scintillator C was glued to two thin lucite strips 
on top and bottom of the scintillator itself. This system was wrapped in Al-foil of 
4 mg/em? thickness. The beam traverses therefore — 8 mg/cm? per turn. 
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for X,> 27 one detects muons mostly of a half-odd integral number of turns. 
The essential features of such distributions, viz. a) the occurrence of a main 
peak accompanied by smaller «satellite » peaks, b) the time separation bet- 
ween peaks, ¢) the widths of each peak, d) the width of the entire group of 
peaks, and e) the relative populations of the peaks, are readily understood on 
the basis of the following considerations: 

Let 0 be placed at the point x, where the muons of momentum 7 cross 
the Y-axis, after having completed an integral number, 7, of turns: 


(4) 0, == 0,=18. 

For such an #,, the T.C. distribution will exhibit a mam peak at: 

(5) tu = NT, 

where T is the orbital period corresponding to p, i.e. & 7.25 ns in our case. 
Since the time of arrival is a function of p, ® and w,: t=t(p, B,x,), à 

spread in momentum and angle of the beam emerging from M, as well as the 


finite width «w» of C (dv,< w< 23) will result in a spread dt in time of 
arrival. In fact, since: 


(6) t=nT and vy = NS 

one has: 

(7) dt/t = dun + dT/T = dn/n + B? dp}p ; (where 6 = r/c) 
(8) dx, /v, == dn/n + ds/s . 


Substituting into eq. (7) the value of dn/n obtained from eq. (8) and making 
use of eq. (3), one obtains: 


(9) dt/t = (6? — A) dp/p — Bd® + da, |x 


wu? 


where 4A and B are defined by the relation: 


(10) ds/s = Adp/p + Bd®, 


with 


1m 
a — 2br D) ? 


and 


= — 2br/(a — 2brD). 
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For example, for a beam with p = 85 MeV/c, d= 0, v, = 8.5 em, |dx, |= 
= w/2 — 1 em, dp/p =10%, AD = 0.1 rad, one finds for the relative spread on 
the T.C. time scale di/t = 50%. The distribution of times of arrival, however, 
is not a smooth function but exhibits a peak structure, as can be seen from 
eq. (8), by noting that dn must be an integer, say dn = m, where m indicates 
the order of the satellite peak. With eq. (9), one can evaluate the separation 
between satellites by substituting into eq. (9) the values of dp/p and d@ that 
satisfy eq. (8) and eq. (10) for a given dn = m, taking into account the cor- 
related distributions in p and ® of the stored beam. 

The width of each satellite is obtained in the same way, by computing 
\dt/t for a fixed n; i.e. dn — 0. From eqs. (8) and (10) we can write: 


(11) dx, >. dn A dp sa 
Ty n p 


If the available dp are large enough, this equation always can be satisfied 
for a given da, or AP, and we have from eq. (9): 


di f?dx, (A_f:)dn BB 


2 (AMATEUR A n A 


d@. 


Thus, the width of the n-th turn (dn = 0) is: 


di & aida 055 ao, 
a A Ly, A 


(the terms to be added in quadrature, of course). For dx, <s 


we have 
Ro 2 o 
se (5 + (haw), 


[1 + (0.1nb):} = ja [1 + (0.22n)?}. 


and d@ = 0.1, 


max 


per 


di = si 


Numerically the widths of the peaks should be: 


die = aly = 0.20.7, 
dt(n = 10) = 0.487, 
di (n = 20) = 0.90 7'. 


The separation between peaks is essentially constant and equal to T because 
of the weak dependence of T on p (eq. (7)). 


2503 


ro 


sani 


296 G. CHARPAK, L. M. LEDERMAN, J. C. SENS and A. ZICHICHI 


To find the overall width of the distribution in t, we use eq. (9) with: 


Di uo, 
? 
dut“ 1 cm 


to find 


dt f "4 (0.22) + (= 
7) AIN ER (3) 


Thus far we have considered the case where the counter C was placed such | 
as to accept muons that have completed an integer number, n, of turns. An| 
analogous derivation of the previous formulae can be obtained for the case | 
where © is set to accept muons after (n-+3) number of turns. In this case, 
one has instead of eq. (6): 


| 
| 
| 
0, = ©, = (m+ 9)8 + 2r ; t=(n+4)7. | 
One finds instead of eq. (9): | 

| 


(13) di/t = (B*— A— 2r/(n + 4)s) dp/p — BAD + da J(n + hs. 


By an argument similar to that following eq. (12), we obtain for the width 
of particular turns: 


der ) == 0.004 7 


role 


dt(n =10 +4) =0.17 7 


di (nm = 20 + 4) =0.36 T 


di(n = 30+ 4) = 0.55 T 


for the counter position of Fig. 3a. The improved visibility of the peak struc- 
ture of Fig. 3a over that of Fig. 2 is thus explained and is, in fact, due to the 
dispersion of muon orbits according to momentum, for the counter in the 
«a, > 2r» position. To find the overall width of the distribution in «t», we 
use eq. (13) and find: 


aD ; 
— = + 0.43 (itn a 0.1; (dp—01 dr, dem, +, — 2r = 92 cm) 


? 
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in good agreement with the observed spread. As we can see from Figs. 2 
and 3a, the overall width of the two distributions is quite different in the two 
eases, this difference being well explained by eqs. (9) and (13) respectively, 
as the two numerical examples quoted above show. Fig. 3a shows an exper- 
imentally observed T.C. picture; Fig. 3b presents a distribution caleulated on 
the basis of the measured momentum distribution and of the expected angular 
Spread due to multiple scattering of the beam in the moderator. 

The stored rate is measured by the area under the peaks in the T.C. distri- 
butions. At an injection (123) rate of 2-10* per second, of which about 600 are 
muons, about 16 muons/s are stored. This ratio 16/600 is in reasonable agree- 
ment with the calculated acceptance. The maximum number of turns de- 
tected with this arrangement was 28. The most important parameter entering 
the quantitative predictions of step-size, period, position of 0 to accept a given 
turn, ete., is the mean momentum p of the circulating muon beam. It can, 
in principle, be calculated from the value of the initial momentum and the 
properties of M. It is, however, useful to measure it directly. For this pur- 
pose we have measured the mean diameter 27 of the orbits by intercepting 
them by means of a lucite block B (Fig. 1) and observing the rate of muons 
arriving in ©. When this block is moved towards M, the rate in C drops for 
I,< 2, = 27 (for the definition of x,, see Fig. 1); after intercepting the first 
turn (at x,=%, =2Y—%;), the ratio increases again until for x,<0, the 
block cuts into late turns once more. It is from the observed Dy, and dy, that 
we conclude that p — 85 MeV/c; this is in agreement with the value expected 
from the momentum loss in M. 

It might be thought to be advantageous to slow the muons to a smaller 
momentum than 85 MeV/c. Because of the rapid variation of ionization loss 
with muon energy, however, the density of muons per unit dp/p falls like 
(p/p)* in the non-relativistic range. In addition multiple scattering increases 
and the resulting loss in intensity precludes the use of such slow muons. 


3. — Electric dipole moment of the muon. 


It has been first remarked by LANDAU (?) that the presence of a static 
electric dipole moment (EDM) for elementary particles can be interpreted only 
as the simultaneous violation of space (P) and time reversal (7) invariance. 
Under a P-operation, the electric dipole moment vector D which can be put 
equal to ko, changes sign whereas o does not; under a T-operation, o changes 
sign but not D. Therefore, in a P-conserving interaction, the absence of an 
EDM would allow no conclusion as to the invariance with respect to time 
reversal and so it is desirable to look for electric dipole moments in particles 
known to have characteristics non-invariant with reflection (?). The deter- 


= 20 - Il Nuovo Cimento. 
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mination of a finite value for the EDM of the muon would provide a lower 
limit for the non-invariance under T-reversal. 

T-non-invariance has been sought in a number of experiments: «) the 
EDM of the neutron (5), was found to be D < 5-10-*°-¢ (cm). b) in the }-decay 
of the neutron, any phaseshift due to T-non-invariance has been measured (°) 
to be less than 0.15 rad. ¢) Recently a series of experiments has established 
an upper limit to the EDM of the muon (*). Our measurements were based 
on the same principle as these latter, as follows: 

Assume that longitudinally polarized muons possessing an EDM given by 


h (4 h 


= Gf 
MC 


RE 2me 2’ 

move in a magnetic field B, in the plane perpendicular to the latter. In their 
instantaneous rest-frame, they experience an electric field E'= yBxB (B = v/e; 
y = (1—f?)*), which causes a precession of D. In the laboratory the spin 
recesses about vx B, i.e. out of the orbit plane by an angle: 


O: PE 2f By O, 7) 
when the orbit has gone through an angle: 


O,= oh t 4) 


MCV 


The polarization (perpendicular to the orbit) thus produced is detected by 
stopping the muons after a known 0; and measuring the up-downa symmetry (*) 
of the decay-positron with respect to the orbit plane. It is important to note 
that the determination of this component of the spin can still be achieved 
if the muons stop in a strong magnetic field. The present determination dif- 
fers from that described in ref. (*) only in the magnitude of ©,: While 
in ref. (*) O,< 2a, we used the storage device described here to obtain 
0, = 2nx, with n up to 28. The arrangement used to detect the decay positron 
in our measurements is schematically illustrated in Fig. 4. This arrangement 
consisted of three counters 4, 5 and 6; the scintillator of counter 4, thick 
enough (3 g/cm®) to stop the muons incident upon it, lies in the median plane 
of the magnet, while scintillators 5 and 6 are arranged symmetrically with 


(5) J. H. Suivre, E. M. Purcert and N. F. RAMSEY: Phys. Rev., 108, 120 (1957). 

(9) M. T. Burney, V. E. Kroun, T. B. Novey, G. R. Ringo and V. L. TELEGDI: 
Phys. Rev. Lett., 1, 324 (1958). 

(7) V. BARGMANN, L. MICHEL and V. L. TELEGDI: Phys. Rev. Lett., 2, 435 (1959). 

() The up-down asymmetry is given by a = EL NS 


down 
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respect to 4, respectively above (5) and below (6) it. Coincidences (45) or (46) 
are produced by the positron from muon decay in 4. The assembly (456) was 
mounted with centre at a distance r=66cm from M. The injection tele- 
scope (123) was as indicated in Fig. 1, and produced a gate of 300ns du- 
ration. 

A (123 gated 4) coincidence corresponds to a muon arriving in 4, while 
we refer to a (45) coincidence as « Top-electron » and to a (46) coincidence as 
« Bottom-electron ». A system of 
gate circuits (logic box) was trig- 
gered by the arriving muon and 
enabled one to count for « Top » 
and « Bottom » separately, « elec- 
trons + accidentals » and « acci- 
dentals ». Upon arrival of a (123 
gated 4)-signal at the input of 
the logic box, two identical sets 
of two gate pulses were produced; 
each set consisted of one prompt 
pulse (called G,; duration 5 us) 
and one pulse delayed by 6 us 
(called G,; duration also 5 us). 
The box had two more inputs: 
one fed by the « Top-electron 


Fig. 4. — Counter layout for the measurement 


: of the electric dipole moment. 
-electron signal». Each of these È oe 


two signals was delayed by 5.5 us 

and then brought into coincidence with one set (G, and 6G,). In this 
way a «prompt» event (with respect to the 4 pulse) in Top or Bottom 
gives no output, since it falls between G, and G,; events occurring be- 
tween +0.5 and 5.5 us (electrons+accidentals) pass through G,; events oc- 
curring between — 5.5 and — 0.5 us (accidentals) pass through G,: The out- 
puts Top-G,, Top-G,, Bottom-G, and Bottom-(, are split into two: each 
of them is brought into coincidence with two single channel pulse-height 
analyser output pulses (A or B). A and B are set to cover different ranges 
of the pulse-height of the output pulse from the T.C. In this way we select 
two ranges of time of flight: group A, particles which have made few turns 
in the magnet (*) and are used for calibration, and group B, particles which 
have made many turns. By subtracting the two asymmetries thus obtained, 


(*) In spite of the large d.stance (x—66 cm) of the EDM telescope from M, some 
muons reach counter 4 in half turn by going through the fringing field of the magnet. 
See Fig. 3a. 
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instrumental asymmetry and an up-down asymmetry in the incident beam, 
are eliminated. One thus measures simultaneously in 8 sealers (see Fig. 5): 


__ from Early muons 
electrons 
+ accidentals 
~~ from Late muons 


Top 
events 
__ from Early muons 
accidentals 
~~ from Late muons 
__ from Early muons 
electron 
+ accidentals — 
uom ~ from Late muons 
events 
-from Early muons 
accidentals 
-from Late muons 
5 ; 


Logic box 


eee 


Coincidence 


8 outputs 
A B 
N pu 
“4 = Wh 
DA N Zi 
Top Bottom Top Bottom 
AMB i. 


Total Random Total Random Total Random Total Random 


Time zero 
a Bau ee 3 

Fig. 5. — Simplified block diagram of the electronics for the electric dipole moment 
measurement. 
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Parallel to this the Top and Bottom events are identified as electron from muon 
decay in no. 4 by measuring the muon life-time in a separate « life-time con- 
verter » whose «startside » is fed by the gated 4 pulse and « stopside» by the 
(4 5) or (4 6) coincidence. 

The results can be divided into 3 groups: Group I in which the up-down 
asymmetry is measured for muons which have made on the mean, either one 
[A («early ») set on <n> —1] or 11.5 [B( late») set on <n> =: 11.5] turns. 
Group IT with mean turn index 1 and 16.5. Group ITI with mean turn index 1 
and 19.5. For each group the difference in asymmetry A between «early » 
and «late » is evaluated. The angles ©; through which the spin has rotated 
out of the median plane is given by ©, = A/a,,., where 4, is the maximum 
obtainable asymmetry in this (150 MeV/c) beam (*). The mean spatial angle 
On = <o(t):t) ( is the time of arrival in 4 enclosed by the windows of A 
resp. B). However, since d7/T=|da/wm|= f*(dp/p), one finds using eq. (13) 
that the dependence of © on t is weak and thus @,,,,,~ @<t>. One therefore 
takes the centre of gravity of the PHA range enclosed by A resp. B. The 
results are listed in Table I. The errors are statistical. 


TaBLe I. — Results obtained for particles which have made, in the mean 11.5, 16.5 and 
19.5 turns (Group I, II and IIT respectively). Geary refers to the asymmetry for muons, 
which have made, in the mean, one turn (ranging from } to $ turns). 


Group Late turns | early ate ate Es Aearly GO) <Qort? 
(%) | (%) (%) | (rad.) | (rad.) 
| | 
to} 64-18 | 8.01 20:84 | 69840.40) —1.0440.94 | 0.062-£0.047 66.0 
ee Cr 281681 07-0422 0:8 7a 0.242.005 0.014260-100 97.0 
TSE 264 | 7:41 41.24 | 7.19+0.94 | —0.22+1.55 | 0.013+0.080 | 116 


A systematic error could be introduced if the circulating muon beam does 
not enter 4 symmetrically with respect to the median plane, or if counters 5 
and 6 have different efficiencies for detection of the decay electrons. By taking 
the difference in up-down asymmetry between «early » and «late », both errors 
are eliminated provided the «early » and «late » fractions of the beam on 4 
can be shown to have the same deviations from symmetric incidence. A ver- 
tical scan on «early » and «late » turns, arriving in 4, shows no difference in 
the centre of the intensity profile to 2mm, while both are ~3 mm above 


() This was measured in Carbon as stopping material. Although the stopping 
medium in this case is scintillating material which is partly depolarizing, it can be 
inferred from data on the re-establishment of the polarization in strong magnetic fields 
that in 13 kG the polarization is restored to 100%. See: J. C. Sens, R. A. SWANSON, 
V. L. TeLEGDI and D. D. Yovanovitcu: Phys. Rev., 107, 1465 (1957). 
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the median plane. No error is thus introduced when taking the differences, 
but a higher top than bottom rate is to be expected. A better check is obtained 
by taking data with the EDM telescope displaced with respect to the median 
plane as shown in Table II. Although the «early » and «late» asymmetries 


Tage II. — Subdivision of Table I, Group I into asymmetries obtained by placing the 
EDM assembly at different height with respect to the median plane of the magnet. 


Height above | Cearly Mate Mate — Dearly 
median plane (mm) (% (92) (05) 

=) 10.14+2.40 9.42 = 1.06 — 0.72 +2.60 

0 10.25 +2.10 7.23 +0.80 — 3.02 42.25 

+ 2 6.85 +1.20 7.12+0.69 0.27 +1.39 

+ 3.5 | 9.37 42.10 5.49 +0.87 — 3.88+2.28 

ae) 4.00 +3.10 5.34+1.38 + 1.34+3.40 


each vary considerably, the variation in the difference has been shown to be 
insignificant. A second systematic error could be introduced by the longitu- 
dinal electric field, which is produced in the rest frame of the particle when 
it slows down in 4. The error is, however, negligible since it depends on the 
EDM to second order (*). As a separate check on the performance of the ap- 
paratus, the time distribution of decay electrons was measured. The life-time 
for u*, deduced from these measurements combining the top and bottom 
counters, was found to be: 

ti = (2.19 + 0.03)-10-* 8. 


U 


Combining the 3 values of O,,,./O.., Which are listed in Table I, one finds: | 


if = (25% = 2.8) es, | 


This corresponds to an upper limit for the EDM of e*(5.0 
is consistent with time reversal invariance. 
This value of «f» is to be compared with: 


5.0):10717 cm, and 


== (2.0 +1.5)-10- 


obtained (*) in ref. (8). 


{ ) It must be noted that by assuming that the incoming beam has zero polari- 
zation transverse in the orbit plane, a value: 


(AES el Po’, 


s obtained in ref. (3). However, no errors are quoted as to the possible magnitude 
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of the up-down asymmetry caused by misalignment of the cyclotron, ete. A partial 
discussion of these points is given by G. Gipaz: Thesis (Columbia, 1960) (Nevis 83), 
Appendix II. We believe that the focal point of the beam transport system in this 
experiment has not been demonstrated not to be above or below the target, which 
would contribute a vertical component of polarization sufficient to vitiate the second 


result quoted. 


RIASSUNTO 


Si descrive un metodo per la cattura di muoni in campi magnetici e la sua appli- 
cazione ad una più accurata determinazione del momento di dipolo elettrico del muone. 
Il valore trovato è: 


EDM <e(5+5):10-17 cm 


consistente con l’invarianza rispetto all’inversione del tempo. 
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Inner Bremsstrahlung-Nuclear Recoil 
Angular Correlation in Æ-Capture (*). 


S. D. BLoom and J. L. URETSKY (**) 


University of California, Lawrence Radiation Laboratory - Livermore, Cal. 


(ricevuto P11 Aprile 1960) 


Summary. — It is suggested that a measurement of the inner brems- 
strahlung (IB)-nuclear recoil angular correlation in K-capture leads to 
a method for determining the ratio |(,|?/|C,|?, which is substantially 
free from dependence on either lifetime or spectral-shape factors, in 
contrast to the estimates of |C,|?/|0y|? made from the ft values. So far 
a determination such as suggested here has been made only in the case 
of the @ decay of aligned neutrons. In the case of 7Be, the only known 
super-allowed K-capturer, a similar measurement might be of considerable 
interest since here some fairly dependable calculations of | Mer], the 
Gamow-Teller matrix element, might be made. (Knowledge of | Worl 
is essential in the evaluation of |C,|?/|0,|?). A general formula is derived 


for the IB-nuclear recoil angular correlation and applied to the particular 
i Æ-capture … . “ne STE 
case of the 7Be —°"5 7Li p decay. In addition, it is shown that the 


formula, which in its general form assumes an initial (arbitrary) polari- 
zation for the A-capturing nucleus, offers the possibility of doing an 
unambiguous time-reversal experiment by looking for a left-right asym- 
metry in the count rate when polarization, recoil-momentum direction, 
and IB-direction are all mutually perpendicular. 


1. — Proposed experiment. 


The discovery of parity breakdown in f-decay (') has given rise to a com- 
plete re-examination of the problem of the form of the @-interaction and, in 
fact, weak interactions in general. In particular, the work of FEYNMAN and 

(*) This work was performed under the auspices of the U.S. Atomic Energy 
Commission. 

('*) LRL, Berkeley; present address: Dept. of Physics, Purdue University, Lafa- 
yette, Ind. 

(1) T. D. LEE and C. N. YANG: Phys. Rev., 104, 254 (1956); C. S. Wu, E. AMBLER, 
R. W. Haywarp, D. D. Hopper and R. P. Hupson: Phys. Rev., 105, 1413 (1957); 106, 


1361 (1957); R. L. Garwry, L. M. LEDERMAN and M. WEINRICH: Phys. Rev., 105, 
1415 (1957). 
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| GELL-MANN (2), SUDARSHAN and MARSHAK (?), and others strongly suggests 
that the weak interaction is of the form (V— A). Moreover the experimental 
evidence, for }-decay at least, very clearly supports this hypothesis, including, 
for instance, experimental redeterminations of the *He recoil spectrum (#) and 
the recent discovery of the (e+ v) mode (4) of decay for the pion. Calculations 
of the ratio of the coupling constants, |C,|?/|C,|2, have also been made; to 
date the most reliable estimates of this ratio are based largely on the ft va- 
lues (5) for super-allowed transitions, using the method proposed and utilized 
by WINTER and KOFOED-HANSEN (°). Although the theories mentioned above (?) 
require that the ratio be exactly one, it is to be expected that in the case of 
complex nuclei corrections due to the virtual pion field will in fact alter this 
situation. In consonance with this the values so far calculated (5) have varied 
rather widely, nor merely from the value one, but also from each other, being 
highest for the neutron decay, 1.42 (5), and going down to a value of around 
1.15 for super-allowed £ emitters of atomic weight greater than about 15. It 
might therefore be interesting to make a measurement of this ratio which is 
independent of either the lifetime or spectral-shape determinations. Actually, 
such a measurement has been done on one case, the decay of the neutron, 
where the correlation coefficient between the B-particle momentum and the 
neutron spin direction leads to a value of 1.56 -- 0.14 for the ratio in ques- 
tion (7). It is interesting to note that this agrees, within statistics, with the 
value obtained from the newest ft value for the neutron decay (?), but the 
possibility of there being a difference between the two numbers is not beyond 
consideration. It is clear that another measurement, not relying on measured 
ft values, might well be a useful addendum to present data. A measurement 
like this would be one such as we propose here, in which one observes an angular 


(2) R. P. Feynman and M. GELL-MANN: Phys. Rev., 109, 193 (1958); E. C. G. Fu- 
DARSHAN and R. E. MARSHAK: Phys. Rev., 109, 1860 (1958); a more recent reference 
containing a reasonably up-to-date bibliography of work in this area is Y. TANIKAWA 
and S. WATANABE: Phys. Rev., 118, 1344 (1959). 

(3) The most recent experimental results for 8 decay may be found in the Minutes 
of the Conference on Weak Interactions Held at Gatlinburg, Tenn., Bull. Am. Phys. Soc., 
4, 1, 77 (1959) and Rev. Mod. Phys., 31, 782 (1959). 

(4) T. Fazzini G. Fipecaro, A. W. MERRISCH, H. PAUL and A. V. TOLLESTRUP: 
Phys. Rev. Lett., 1, 247 (1958); G. ImPEDUGLIA, R. PLANO, A. PRODELL, N. SAMIOS, M. 
SCHWARTZ and J. STEINBERGER: Phys Rev. Lett., 1, 249 (1958); H. L. ANDERSON T. 
Fos, R. H. MiLLER, and L. Tau: Phys. Rev. Dett., 2, 53 (1959). 

(5) O. C. Kistner and B. M. Rusrap: Phys. Rev., 114, 1329 (1959). 

(6) A. Wrntuer and 0. KoroeD-HANSEN: Kgl. Dan. Videnskab. Selskab-Fys. Medd., 
27, no 14 (1953). 

(7) MoT. Burey, V. E. Kroun, IT. B. Novey and C. R. Rixco: reported by 
M. GoLDHABER at the Annual International Conference on High-Energy Physics at CERN, 
July 1958, edited by B. FERRETTI (Geneva, 1958). 
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correlation between the nuclear recoil and some other detectable particle ac- 
companying the @-decay. The electron-neutrino correlation, or, what is the 
same thing, the electron-nuclear recoil correlation for super-allowed transitions 
of the (J+J, T >T) type, wherein both Fermi and Gamow-Teller matrix 
elements may contribute to the decay, would fall fairly well into this category. 
This is true in spite of the fact that the correlation is still somewhat dependent 
on the 6 spectral shape due to the finite mass of the electron (8). However, 
radio-isotopes of this class that are easily fabricated and also amenable to 
nuclear-recoil detection, and, morevoer, that have lifetimes long enough to 
accumulate sufficient data, considering the reduced count rate imposed by the 
coincidence requirement implied above, are not only few and far between, but 
also possibly even non-existent (*), so far. At any rate, K-capture offers the 
possibility of doing the equivalent of an electron-neutrino correlation exper- 
iment, but with the electron replaced by the inner bremsstrahlung (IB) radia- 
tion, which radiation accompanies the K-capture during 10-' to 10-* of the 
time. Furthermore, the zero mass of the IB photons offers the advantage of 
eliminating an energy-spectrum dependence present in the case of the ordi- 
nary electron-neutrino correlation. 

The formula for the IB-neutrino correlation (cf. the following section) is 


(1) W(g) = Constant-v?(k,, — »)[1— 0, cos y], 


where: p = the angle between the neutrino and the IB photon, 
y = the neutrino energy, 
len 2 the maximum photon energy, 
dQ, = the photon solid-angle differential, 
dQ, = the neutrino solid-angle differential, 
Fai rl Lai se) Or|*) 3 | Men |*(| Cx |*—| 041?) 
| My |?(| Os |? +] Or|?) + | Mow |?(| Or |? + | Cu |?) ’ 


M,,= the Fermi matrix element, 


M, = the Gamow-Teller matrix element, 


C, = the coupling constant as per subscript X. 


(*) J. D. Jackson, 8. B. TREIMAN and H. H. WyLp, Jr.: Phys. Rev., 106, 517 (1957). 
(*) Nevertheless it is worth while mentioning 19Ne, on which recoil spectrum measu- 
rements have been recently made by HERRMANNSFELDT et al. (cfr. ref; (3), op: 77). 
Its high 6+ end-point energy (2.18 MeV) minimizes the spectral dependence factor, 


but its 18s half-life requires the extended use of an accelerator for manufacturing 
the radioactivity. 
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To go from this formula to one for the IB-nuclear recoil angular correlation 
is a straightforward problem in kinematical and solid-angle transformations. 
The final result is 


n P2ky | K+P cos O 
2 W(0) = Constant: enna Pre KO KO 
(2) (0) onstan era ae 5 || L 00 ik fue d2,d2>, 
where: P = the nuclear recoil momentum, 


© = the photonuclear recoil angle, 


k =the photon momentum. 


It is to be noted, just as in any other type of 8 experimentation, one does 
not escape here the inevitable linking together of the coupling constants and 
the matrix elements in the form 0, M,. It is for this reason that one is most 
inclined to work with a super-allowed decay rather than one that is merely 
favored, since in the former case one might be able to calculate values for 
the Gamow-Teller matrix elements, either by the magnetic moment method 
or by direct computation (51%), Such a super-allowed K-capture is not merely 
rare but unique, namely, “Be. The radio-isotope "Be is unique also in that 
a good deal of recoil detection work has been done with it (11), so that a con- 
siderable knowledge has already been built up with respect to that aspect of 
the problem. Furthermore, the half-life of 53 days, plus the known methods (1?) 
for producing copious amounts of the isotope, make it experimentally attractive 
from the handling point of view, since the accelerator operation produces the 
source and the source measurement itself is quite separate in the length of time 
available. 

In Fig. 1(a)-(d) we have plotted the 7Be IB-nuclear recoil correlation in 
the case of a pure Fermi transition, pure Gamow-Teller transition, and also 
for exact equality of |0,|?|M@,,|? and |C,|?|2f,|?, assuming the (V— A) form 
for the 8 interaction. As can be seen the differences between the three cases 
are quite pronounced. Therefore it is to be hoped that with the advanced 


(10) A. M. LANE: Proc. Phys. Soc. (London), 68, 189 (1955); in the Appendix of 
this paper the work of LANE on the 7Be +"Li @ decay (which is of particular interest 
here) has been brought up to date and compared with more recent calculations 
(see ref. (5)). 

(11) P. B. Suirx and J. $. ALLEN: Phys. Rev., 81, 381 (1951); R. Davis: Phys. 
Rev., 86, 976 (1952). 

(13) F. AJZENBERG and T. Lauritsen: Rev. Mod. Phys., 27, 77 (1955). 
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techniques available today a measurement of the sort proposed here is not very 
far off. i. ; 
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Fig. 1. — Relative IB-nuclear recoil angular correl 


ation in K-capture for four photon 
energies, assuming the V—A interaction form. 


It is to be noted for photon energy 

greater than %/km = 0.5 there is a minimum angle, Omin = cos-1[(#2 —v)2/k], beyond 

whick no correlations can be observed, due to simple kinematic relations. Also note 

all curves have been arbitrarily normalized to the same value for 0,, (angle between 

IB photon and neutrino) equal to 90°. The three cases plotted for each energy are 

labeled with regard to the assumptions made, « intermediate » meaning absolute equality 
of |C,|?|Mr|? and | 0, |? | Mor]? 
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2. — Theory. 


The theory of internal bremsstrahlung in K-electron capture has been dis- 
cussed in much detail by GLAUBER and MARTIN (1%). Although our calculation 
follows along the same general lines, there are two points requiring comment: 
1) We must be a little more careful with the nuclear kinematics (because we 
are interested in observing the nuclear recoil), and 2) we can take advantage 
of GM’s results which are based on the assumption that Coulomb effects in 
intermediate states are negligible. (This treatment is appropriate for high- 
energy y-rays to which our interest is confined.) 

For the convenience of the reader we repeat some definitions of the GM 
paper. The radiation and 8-decay Hamiltonians are, respectively, 


(3) H'(x) = eyp(w) A(x) y(2) 
and 
(4) MEN — >; CO, Pl?) Ty ,Pp(@)* x(a) Triy(a) = x(a) t(x) p(x) . 


In these expressions we have introduced the field operators for the electron (y), 
nucleon (gx), neutrino (x), and electromagnetic (A,) fields, as well as the 
6-coupling operators T,, with the concomitant coupling constants O,. The 
amplitude for photon emission accompanied by electron capture comes from 
the second-order S-matrix element, 


— 
OL 
— 


Sa — Java (H"(2), H°(Y))+ ’ 


where we have used a conventional notation for the time-ordered product. It 
is to be understood that we are interested in the matrix element of eq. (5) 
between an initial state with one bound electron of total energy H, and a final 
state containing a photon of 4-momentum k, and a neutrino of 4-momentum »,. 

The space-time integrals in eq. (5) are evaluated by introducing the field 
operators for the initial electron state: 


| pr (x) = p(x — Y) exp [—i#,a]a, = 


(6a) | = mIa fap) exp [i[q: (x — Y) — iE,%]], 


(13) R. J. GLAUBER and P. C. MARTIN: Phys. Rev., 104, 158 (1958); hereinafter 
this reference will be quoted as (GM); also efr. P. C. MARTIN and R. J. GLAUBER: 
Phys. Rev., 109, 1307 (1958). 
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which involves the co-ordinate Y of the nuclear center of mass, and the final 
neutrino state: 


(6b) a (27) (È) fs(v) exp [— 4 vi x]as(y) . 
0 
The result is 
e [m,\? 
193 - k Ie . 
” ASI a 


Jar d'gò(po—v— Pour T(p — v) Sr(p, ¢ — *) ¢ y-(Q) = 
= O(k + Vv DE P) 0(E,— ko — vo— P)M; 5 


in which we have introduced the causal propagator S,(p, q) for an electron 
in the presence of the nuclear Coulomb field. In order to reach eq. (7) it was 
necessary to assume that the wave function for a bound electron-nucleus 
system could be separated into electronic and nuclear co-ordinates. This is 
equivalent to assuming an infinitely heavy nucleus, and it is expected that 
the correction terms should be proportional to the ratio of the electronic to 
nuclear masses. A second approximation is made in performing the sum over 
intermediate electron states; this is necessary to obtain the propagation func- 
tion, S,, since states are included in the summation that are occupied by other 
electrons. In the GM paper it is argued that this procedure is correct because 
of the possibility that the intermediate-state electron is captured before the 
initial-state electron makes a radiative transition. At any rate, we do not 
expect such considerations to be important in the present case where the y-ray 
energy is large compared with the electronic binding energy. 

It was shown in the GM paper that the effect of the Coulomb field upon 
the intermediate-state electron is not important for large y-ray energies 
(k > Zam,). Since this is the case in which we shall be interested, we may 
expand the propagator in powers of the external field and drop all but the 
first term. We have then 


(8) SR(D, 4 —k) = 0(p—q+k)(ip + m)-! + O(Za) . 


It is consistent with this approximation to use a non-relativistic form for the 
initial electron wave function and to evaluate the integrals in eq. (7) only to 
the lowest order in the parameter Zx. For allowed transitions T(p — v) is 


Dupin of its argument, and the integrals give, when evaluated for an 
initial (1s) state, 


: CPs 0 my E 
(9) M,=x Ss ( 7 | BV) Pp — v= 0): (ih + iy, E, + m,) ¢ Q, . 


Kym. J\2kavi 
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In the last expression q,,(0) is the non-relativistic electronic wave function 
evaluated at the origin and © is the normalized spinor eigenfunction of y, 
satisfying 


(10) 142, = Q,. 


We next square the expression in eq. (9) and sum over the initial electron 
spins (the (1s) shell is assumed filled) and the final neutrino spins, and average 
over final nuclear orientations. After multiplying by the energy-momentum 
ò-funetions and integrating over the neutrino variables and the magnitude of 
the nuclear recoil momentum, we are left with an expression that is appropriate 
to an experiment carried out with initially aligned nuclei. The angle between 
the recoil nucleus and the photon is to be observed along with the magnitude 
of the photon momentum. If the reference magnetic field is chosen as the 
Z axis, then the transition probability is 


LR K(k nm © = 
7’ CS Eee 
AD Wo Be cong UH dQ,dQply: 
{;_-fe+k-P È M? \[k:-3k+k-P— dala] 
LÉ of ee o ©: + 
MAT) 3 JJ+1 K(Kem — k) 
+ KP, 
9 i * VE. * Kin kz 
(2/1) Re M;( Mer) Pz. Me Ole Km fra k) Sie 
| vk ALP Es 2k) cr kP, | ER de (On (kXP), | 
di» 152 de Hi] 


The relevant definitions are 
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The factor | M,,|? is the square of the Gamow-Teller matrix element after 
averaging and summing over initial and final orientations. 


(Mor) = (I * Mr (M/J) ’ 
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where -%@,, is the Z component of the Gamow-Teller matrix element evaluated 
in the state where M — J. The maximum possible photon energy is denoted 
lp Wve 

When the initial nuclear orientation is not observed we obtain the simpler 
relation (given in slightly different form in eq. (2)): 


ke + k-P 
bk k) © 


W |P PACE, = k) 
| P+keos@ 


jar dQ JON TES 


The equivalent of expression (12) has been previously obtained by YANG. 
HUANG, and LEE (#) for the special case of u-meson capture by a nucleon. 
The relative dependence of W on the angle between k and P is shown in 
Fig. 1(a)-(d) for the various situations described above. 


3. — Time-reversal experiment. 


It is worth while to make some additional remarks on the form of eq. (0 
which is the correlation func- 
tion for an IB-recoil experi- 
ment where the decaying 
nuclei are polarized along 
the Z axis. If we assume 
% the complete validity of the 
sd 0 V— A theory and further 

detector | postulate an experimental 
configuration such as given 


direction of 
(Be Polarization 


Time reversal ih 
experiment 2 


ao 


( 
/ 


k 
SE Inner 


Bremsstrahlung in Fig. 2, where the spin po- 
detector ; 5 2 : 
larization direction, the IB 
Fig. 2. — A schematic representation of a time- 


ì : 3 k direction (k), and the recoil 
reversal experiment with polarized 7Be nuclei. The 


solid and dotted vectors for k correspond to the po- 
sitive and negative signs, respectively, in eqs. (13) 
and (14). 


(3) Wat 90°) oc PA(k,, — k) dkdQ, 49,1": 
n Ji k 0 1 AM x alle (—)7-4#-1 
es | en ye ty kJ 
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direction (P) are all mutually 
perpendicular, eq. (11) takes 
a very simple form: 


a (Cy: Car — Cyn Cut) | 


| ? 


Ly 


(14) K. Huane, C. N. Yane and T. D. Luz: Phys. Rev., 108, 1340 (1957). 
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where the I and R subscripts denote, respectively, the imaginary and real 
parts of the coupling constants. If we now concentrate our attention simply 
on the relative count-rate ratio for + 90° (plus and minus are as defined in 
Fig. 2), we find 


E of 1 M? 
4 F 90° | Lag SALO || VB 
( ) W paul DIE 90 )oc il = Km i, k (20 ai 8 IHR DÌ sui 


k 2PM My.A or fs 
et ie AU sai "(Cy Car — Cnt) , 
m e 0 


It is easy to see that roughly speaking the east-west difference in count rate 
will be of the order —2(0,+0,)/(10,1+10,]) Here we have taken advan- 
tage of the fact that in super-allowed @ emitters the matrix elements 
are necessarily diagonal, i.e., real. This is important because it means the 
existence of an effect is unambiguous proof of the breakdown of time invariance 
in $-decay, since such an effect would require non-zero values for C,, and/or 
C,,- Conversely, it also means the absence of an effect would be evidence, to 
the extent of the precision of the experiment, of the validity of time-reversal 
invarianee in -decay. 


The authors would like to express their thanks to Mr. DOUGLAS VERDERY 
for his invaluable help in preparing the various versions of the manuscript 
and to Dr. Ervin M. ScHwarcz for his interest and comments during the 
progress of the writing. 


APPENDIX 


In the case of the ground-state decay of 7Be, LANE (1°) has made a cal- 
culation of the Gamow-Teller matrix element (the Fermi matrix element being 
simply unity) as a function of the intermediate coupling parameter, ajK. 
In Fig. 3 we have replotted Lane’s results so that | Mgr|* may be directly 
read off as a function of y or (a/K)[y = (a/K)/(5 + a/K)]. Although Lane com- 
pared theoretical ft values (as computed by him) with experimental values, 
it was necessary for him to insert in his theory a value for the Fermi coupling 
constant, an improved number for which has become available since then (1955). 
However, in the case of 7Be or "Li the values over which a/K might range 
are restricted by many more consistent data than just the ft value. Accordingly, 
in this calculation we will take Lane’s estimates for this range: 2 < a/K < 4. 
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In Fig. 3 the range of values of | Mur]? corresponding to this range in alk 
can be seen to be 1.43 to 1.53. If we use the value of 1.48+0.05 for | Mer]? 
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Fig. 3. — A replot of Lane’s 
(see ref. (1°)) theoretical com- 
putation of |Mgr|® as a fune- 
tion of a/K, the intermediate 
coupling parameter. Actually 
the abscissa is laid out in units 
of y, where y =(a/K)/(5+a/K). 
From other pieces of data, re- 
lating not only to "Be and "Li 
but other p-shell nuclei as well, 
Lane concludes that a/X must 
lie somewhere between 2 and 4. 
As can be seen from our graph 
this leads to a rather restricted 
range of values of | Mgr|?, na- 
mely, 1:43 to 1053) 


in fitting to the Y-X plot of Kistner and RUSTAD (5), we obtain the results 
shown in Fig. 4, where it can be seen that the point for the 7Be 8-decays fall 


Fig. 4. — A duplication of the Y-X plot 
of Kistner and Rustad (see ref. (°), 
which describes in detail how this re- 
presentation is achieved). With the 
exception of the value of 0? obtained 
from the pure Fermi transition (0-0), 
the graph does not show the experimen- 
tal points they used in order to pick the 
three values for O%,/O} leading to the 
three different lines. It is clear that a 
value of C2,/C% of around 1.16 is ne- 
cessary to account for the ft value of 
"Be, just as in the other super-allowed 
cases of 150 and 17F, 
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essentially right on the Y-X line for the super-allowed B-emitters 150 and “FP. 
Cér/0% of 1.16. In view of the fact that 


This line corresponds to a ratio for 
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{Me 2 is so very insensitive to large changes in a/K this agreement must 
be considered as more than fortuitous, therefore lending good support to the 
validity of the calculated value for | Mew? 3 


RIASSUNTO (*) 


Si fa Vipotesi che la correlazione angolare del rinculo (IB)-nucleare della brems- 
strahlung interna nella cattura A porti ad un metodo per determinare il rapporto 
|0,|?/|C,|?, che è sostanzialmente privo di dipendenza sia dalla vita media sia dai 
fattori di forma spettrale, in contrasto con le valutazioni di |C,|?/|0y|® dedotto dai 
valori ft. Sinora una determinazione, quale quella qui suggerita, è stata eseguita solo 
nel caso del decadimento 8 di neutroni allineati. Nel caso del "Be, che è il solo ele- 
mento superdotato di capacità di cattura Æ, una misurazione simile potrebbe essere 
di considerevole interesse poichè in questo caso si potrebbero effettuare alcuni atten- 
dibili calcoli di | Mr], l'elemento della matrice Gamow-Teller (La conoscenza di | Mur 
è essenziale alla valutazione di |0,]|?/|0y|?). Si deduce una formula generale per la 
correlazione angolare del rinculo IB-nucleare e si applica al caso particolare del deca- 
dimento & "Be catturaX 77 i Inoltre si mostra che la formula, che nella sua forma 
generale postula una polarizzazione iniziale (arbitraria) per il nucleo che esegue la 
cattura K, offre la possibilità di effettuare un esperimento di inversione del tempo, 
privo di ambiguita, ricercando una asimmetria sinistra-destra nel conteggio quando 
la polarizzazione, la direzione dell’impulso di rinculo e la direzione IB sono tutte reci- 
procamente perpendicolari. 


(*) Traduzione a cura della Redazione. 
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Summary. — This paper is a preliminary report on the emission of heavy 
fragments in nuclear disintegration. It is assumed that the probability 
of emission of a heavy fragment depends on the nuclear temperature (7) 
and the potential barrier (V) of the excited nucleus. It is the purpose of 
this investigation to study the dependence of the differential energy 
spectrum on the excitation energy. The data collected were from four 
stacks of Ilford G-5 emulsions, exposed to (a) a stopped pion beam; 
(b) 400 MeV proton beam; (e) (4-5) GeV pion beam and (d) cosmic rays 
with an exposure in Texas, Magnetic latitude of 41° N, and on the Island 
of Guam at the magnetic equator. In this report emphasis will be placed 
on the emission of 8Li as the heavy fragment. The results of the expe- 
riment indicate that for high excitation energy there is a collimation of 
lithium and beryllium fragments in the forward direction which may 
be explained on the basis of evaporation theory by assuming the fragment 
emission to take place during the flight of a highly excited nucleus. In 
the case of low excitation energy the angular distribution of SLi fragments 
is isotropic and the energy distribution indicates that the majority result 
from the absorption of a 7--meson by a light nucleus. 


(*) Work supported by the National Science Foundation. 
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1. — Introduction. 


Whenever a high energy nucleon collides with a nucleus there occurs within 
a time interval of about 10-2% seconds a series of collisions between this nucleon 
and the individual nucleons which make up the nucleus. These collisions can 
be classified into two main groups, radiative and elastic. In all present day 
theories it is assumed that all collisions are of the nucleon-nucleon type since 
a nucleon cannot interact with an assemblage of nucleons as a whole. This 
argument appears reasonable in view of the low binding energy of the nucleons, 
and also since the wave length of the incident nucleon is small compared with 
internuclear distance and the range of nuclear forces. However, the exper- 
imental observation of energetic fragments of helium, lithium, beryllium and 
even boron is inconsistent with this simplified hypothesis. Also the rare 
observation of a «clean» jet indicates that collisions between the incoming 
nucleon and a nucleon in the nucleus must cause a great disturbance in the 
nucleus. 

Throughout the random exchange of energy by collision we can assume 
that after some time a definite amount of energy has been transferred to the 
nucleus leaving it in a highly excited state and that all the fast particles (gray 
tracks) have already left the nucleus. We can assume that the exchange of 
energy among the nucleons in the nucleus is completed before the expulsion 
of any further particles. The disintegration from this excited state may then 
be treated as quite independent of the process by which it is formed and the 
disintegration that follows is analogous to a process by evaporation. If this 
is the case, the energy spectrum of the emitted particles is of Maxwellian form, 
and their angular distribution should be isotropic. 

In this paper we will consider only those particles finally emitted with ener- 
gies below 30 MeV and they will be observed as « black » tracks and have 
grain density greater than 6.8 times minimum. 

The statistical treatment of nuclear disintegration was initiated by Bonr 
and KALCKAR (1), amplified by BETHE (?), WEISSKOPF (3) and LECOUTEUR (*), 
and leads to some very definite predictions. If 1) the excited nucleus is always 
in approximate thermodynamic equilibrium, 2) the emission time between 
successive particles is longer than the time required for the convection of 
energy across the nucleus, 3) the excitation energy is much less than the total 
binding energy of the nucleus, then the theory predicts that in the center of 


G) N. Bonr and F. Karckar: Kgl. Dan. Vid., Selsk., 24, 10 (1937). 
?) H. A. Berne: Rev. Mod. Phys., 9, 69 (1937). 

(3) V. WEISSKOPF: Phys. Rev., 52, 295 (1937). 

(4) K. J. LE COUTEUR: Proc. Phys. Soc. (London), A 63, 259 (1950). 
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mass system of the evaporating nucleus the probability of a particle being 
emitted with kinetic energy between Æ and £-+d£ is given by 


Te 
T2 


(1) pag == exp| To jar, 

where the parameters V and T are the average height of the effective poten- 
tial barrier and the mean nuclear temperature in MeV. 

| During the last years numerous reports have been given of the observation 
of heavy nuclear fragments emitted from stars in photographic emulsions (519). 
The study of nuclear disintegration with emission of particles with charge 
Z — 3 has been extensively studied in our laboratory. The program consisted 
in studying the emission for various excitation energies. The exposures ranged 
from high energy cosmic rays at the magnetic equator to 140 MeV energy 
realized in stopped =~ events. 


2. — Experimental procedure. 


21. Exposures. — Various stacks of Ilford G-5 emulsions were exposed to 
radiation of different primary energy and disintegrations, accompanied by heavy 
fragments, were analysed. These primary beam energies were, 1) stopped 
x beam, 2) 400 MeV beam from the University of Chicago cyclotron, 
3) 4.5 GeV pion beam (*) and 4) Cosmic Ray at the magnetic equator, Island 
of Guam. 


(9) A. BonETTI and C. Ditwortn: Phil. Mag., 40, 585 (1949). 
(9) P. HopGson and D. H. Perkins: Nature, 163, 439 (1949). 
(7) D. H. PERKINS: Proc. Roy. Soc. (London), A 203, 399 (1950). 
(8) S. O. SORENSEN: Phil. Mag., 42, 325 (1951). 

(°) J. P. LoncHaAMP: Ann. Phys., 10, 201 (1955). 

(9) B. A. Muir: Phil. Mag., 1, 355 (1956). 

(Ge) S 


5. J. GoLDsAcK, W. 0. Lock and B. A. MuNIR: Phil. Mag., 2, 149 (1957). 
(1) S. NAKAGAWA, E. Tamar, H. Huzita and K. Oxuparra: Journ. Phys. Soc. 
Japan, 11, 191 (1956); 12, 747 (1957). 

(13) 0. SKJEGGESTAD and S. O. SÜRENSEN: Phys. Rev., 118, 1115 (1959). 

(15) S. Karoorr: Phys. Rev., 114, 905 (1959). 

(5) R. HAGEDORN and W. MacKkE: Komische Strahlung (Berlin, 1953), p. 201. 

(19) R. H. Brown, U. CAMERINI, P. H. FowLER, H. Herter, D. T. Kine and CA 
F. POWELL: Phil. Mag., 40, 862 (1949). 

(5) M. TEUCHER: Komische Strahlung (Berlin, 1953), p. 69. 

(18) S. NAKAGAWA, E. Tamar and O. Nomoro: preprint (1958). 

(19) E. Tamar: Nuovo Cimento, 14, 1 (1959). 

(*) This stack of emulsions was kindly loaned to us by P 


dl i rofessors V. TELEGDI and 
R. Levi SETTI of the Enrico Fermi Institute of Nucle 


ar Study, University of Chicago. 
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The criteria for the selection and measurements of tracks were as follows: 
1) Tracks of all particles ended in the emulsion. 


2) The angle of dip was less than 30° before processing and range cor- 
rections were made for all angles. Width measurements on tracks with dif- 
ferent dip angles were made and found to be consistent up to 30°. 


3) The complete thickness of the emulsion was scanned. There was 
no significant change in track width from one emulsion depth to another for 
particles within these limits. 


4) All hammer tracks were assumed to be due to SLi nuclei stopping 
in the emulsion. The stacks for all exposures were gridded which permitted 
tracing back of all observed hammer tracks. Hence, the geometrical cor- 
rection to the energy spectra for tracks escaping from the emulsion is neg- 
ligible (?°). 


5) The energies of the fragments were derived from the measured lengths 
of tracks and range-energy relations for *Li in nuclear emulsions given by 
SKJEGGESTAD (2°) and BARKAS (°°). 


22. Determination of charge of the fragment. — The charge of heavy frag- 
ments was determined by the method employed by NAKAGAWA et al. (12). 
The procedure consists in measuring the average width of stopped particles 
as a function of residual range. By means of a bifilar the average width was 
obtained every 6 um of length and when plotted as a funetion of residual 
range produced a graph showing a plateau, the ordinate of which is taken as 
the width of the track. Our calibration consisted in measuring the track width 
of known particles. These were protons from stopped proton exposure, « tracks 
from thorium stars and hammer tracks, and Li from #Li hammer tracks. In 
the case of lithium fragments the flat portion of the width curve started at 
about 38 um, for alpha’s at about 30 um and for protons at about 18 um. 
In the identification of fragments it was only necessary to measure the width 
at every 10 um. The weighted average width of the *Li fragments was ob- 
tained by finding the centroid of the histogram or 


(1) en Dn; 


Assuming that the width of a track of charge Z varies as Z*, the average width 


(2) O. S. SKJEGGESTAD: Nuovo Cimento, 13, 927 (1958); Thesis (Oslo, no. 289). 
(21) W. H. BARKAS: Phys. Rev., 89, 1019 (1953). 
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Fig. 2. — Results of the study of track width of lithium fragments of four different 
range groups. 
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of any particle of charge Z was obtained and is indicated by arrows in Fig. 1 
and the corresponding width with standard deviations is given in Table T 


. 


TABLE I. — Average track width for fragments of charge Z=1 to Z=8. 
Z Average width (um) | Z Average width (um) 
eel 0.63 £0.05 5 1.39 0.03 
| 2 | 0.89 +0.03 6 1.52 | 
sae 1.07 +0.05 vert 1.65 
ike oho re} 1.24 0.04 lE 1.76 


It is immediately evident that for small Z the method employed becomes 
a fast and simple way for the determination of charge, but becomes less re- 
liable for Z greater than five. A study of the reliability of track width measu- 
rements as a function of the length of the track was made for "Li fragments. 
The tracks were classified into four groups of ranges up to 500 um. The results 
are given in Fig. 2 and show that the width of the track is independent of 
the range. 


3. — Results. 


31. Li, ‘Li and SL, fragments from cosmic ray stars. — A study was made 
of the energy distribution of the unstable SLi and the stable "Li fragments. 
Tf the lithium fragments are emitted during the evaporation of a highly excited 
parent nucleus, then their differential energy spectrum should be of the form (1) 


E—V 
pl 


jar, 


for E > V and where V is the potential barrier and 7 is the nuclear tempera- 
ture both expressed in MeV. 
The average £ for any fragment will be given by 


[E-N(E) dB 


| N(H) 48 
0 


where the value of V for heavy nuclei was taken to be (1) 


(4) V=5.5Z (MeV). 
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For each excitation energy £ (in MeV) the nuclear temperature 7 was calculated 
from the excitation energy H(exc) 


(5) E(exc) = KAT?, 


where A is the atomic weight of the excited nucleus, and A = 7?/4G@, where 
( is the kinetic energy of the highest occupied state in the unexcited nucleus 
and is taken as 25 MeV. Equation (5) sets an upper limit on the nuclear tem- 
perature T of 10 MeV, since the application of the evaporation theory to 
heavy nuclei requires that the excitation energy be less than the binding 
energy of the nucleus or about 800 MeV. 

The distribution for *Li and $’Li is given in Fig. 3. In all our work it was 
assumed that the evaporating nucleus came to rest before the fragment was 
emitted. In this figure we give the distribution of the lithium fragments ob- 
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Results of Barkow ETAL, 


Kinehc energy in MeV 
Results of Sorensen ETAL 
Fig. 3. — The energy distribution of ®7Li and *Li from cosmic ray stars. The angle 
of dip is less than 30°. The 6”Li distribution is for 340 events giving T—15 MeV, 
V=12.5 MeV and £ (av)=44 MeV. The SLi distribution is for 467 events giving 
T=9.5 MeV, V=5 MeV and E (av)=24.5 MeV. In both distributions T and V were 
obtained by the method of hest fit. 


served by us and those found by SKJEGGESTAD and SORENSEN by the method 
described in reference (1). By means of a projection microscope they draw 
contours of the fragment tracks and determine the charge by measuring the 
total area of the last 106 um of the track with a planimeter. SKJEGGESTAD 
and SORENSEN concluded that the marked difference in the form of the two 
spectra, position of the maxima, is not real, They believe that the lack of 
any strong increase in number of stable lithium fragments in the low energy 
region is due to a considerable loss of short tracks during scanning and the 
bias imposed by the method used. 
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Our results bear out these conclusions. Since SKJEGGESTAD and SORENSEN 
required tracks of length 106 um or greater, their energy cut-off was in the 
neighborhood of 29 MeV. Employing the method of measuring track width 
we are sure of identifying particles of charge three, having a range of 38 um 
or more which corresponds to a cut-off energy of 16 MeV. The inability of 
definitely identifying lithium particles of mass 6 and 7 of range less than 38 um 
introduces an error in the low energy portion of the distribution. Any increase 
in number in these two intervals would result in a shift of the theoretical curve 
to the lower energy region and should give complete agreement between the 
different isotopes. At the present time profile measurements are being made 
on short tracks in an effort to actually identify lithium fragments of range 
shorter than 38 um. 

From the curves the experimental results can be explained qualitatively 
by assuming that the fragments of lithium of kinetic energy less than 80 MeV 
are ejected as part of the evaporation process. For energies greater than 
80 MeV there is a surplus of fragments over the number predicted by the 
evaporation theory. 

The values of 7 and V were obtained by getting the best theoretical fit of 
Eq. (1) to the experimental histogram. This was done by making a plot on 
transparent paper of (E —V) vs. N(£)d# normalized to the total number of 
observed particles, for various values of 7. By superimposing the theoretical 
graphs on the histogram the best values of 7 and V can be obtained. For 
6.7Li the values obtained by the method described above were T= 15 MeV 
and V=12.5 MeV, while employing Eq. (2) and (3) we get T— 14 MeV and 
V =16.5 MeV which gives 
reasonably close agreement 257 2 in 
and an Æ,, — 44 MeV. The 8-16" 
E,, is biased toward the high 20] pr e 
energies. However for SLi the 
values obtained were 7 = 9.5 15 | 
MeV, V=5 compared with 
T=4 MeV and V= 16.5 and 107 
an Æ,,— 24.5 MeV. 

A study was made of the 5) 
energy distribution of *Li as 
a function of the size of the 
parent star. As a measure of 


the size of the star we used ; Peer see 
: : Fig. 4. — Energy distribution of SLi from cosmic ray 
n heavily ion- 8 à 
ie OL GI ù A stars as a function of the number of black tracks, 
LELE DECODES N, define do N,, of the parent star. The classification was: 
the number of tracks with 4) (1-7) prongs; b) (8--16) prongs and c) greater 
grain density greater than than 17 prongs. 
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1.4 times minimum. Further we have divided the *Li tracks into three groups 
with 1<N,<7, 8< N, <16 and N,>17. The three energy spectra are 
plotted in Fig. 4. As seen from the distribution the form of the three spectra 
does not depend on the energy released in the parent star. 


32. SLi fragments from 4.5 GeV pion stars. — A similar study to that in 
Section 31 was made on “Li fragments found in 4.5 GeV pion stars. The 


energy distribution of 363 lithium fragments is plotted in Fig. 5 and the energy 
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Fig. 5. — The energy distribution of 8Li Fig. 6. — Energy distribution of SJ1 from 
fragments found in 4.5 GeV pion stars. 4.5 GeV pion stars as a function of N,;. 
The distribution is for 363 events, The classification was: a) (1--7) an 
giving E, =22.5 MeV, T=11 MeV and b) (8+16) prongs and c) greater than 
V=8 MeV. 17 prongs. 


distribution of the fragments as a function of N, is given in Fig. 6. The energy 
distribution gives an H,, = 22.5 MeV which is in close agreement with the 
ne 24.5 MeV obtained from cosmic ray stars. The corresponding values 
T=9.5 MeV and V= 5 MeV for SLi fragments from cosmic ray stars must 
be compared with T=11 MeV and V=8 MeV for 4.5 GeV pion Stars. The 
energy distribution as in cosmic ray stars is independent of N,. 


È; D par nr y n n ù na ; 
33. Beryllium fragments from cosmic ray stars. — The energy distribution 
of °Be fraeme à ie pe A Oo lor 3 oe ly | 
a Een Ho Cosmic ray stars is given in Fig. 7. The angular distri- 
bution is given in Fig. 8. The observed values of 7=19 MeV and V=25 MeV 


appear to be in good agreement with evaporation theory values obtained from 
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eq. (2) and (3). The theory gives values of T—21 MeV and V= 22 MeV. 
However it is generally assumed that the evaporation theory breaks down 
when the energy of excitation approaches the total binding energy of the 
nucleus. If we assume an average A of 100 for AgBr and a binding energy 


NT & 
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Fig. 7. - The energy distribution for 98 *Be fragments from cosmic ray stars, with 
the following results: E,, —64 MeV, 7=19 MeV and V=25 MeV. 


of 8 MeV per nucleon, the upper limit will be of the order of 800 MeV cor- 
responding to a nuclear temperature of about 10 MeV. The value of 21 MeV 
for the mean temperature is much too high to be acceptable. This conclusion 
is even further supported by the results of nuclear ‘ascade calculations which 
vive the average excitation energies of the residual nuclei of only a few hundred 
MeV. We are at a loss 

to explain the relatively 20 
large number of ‘Be | 
fragments with kinetic 

energies less than 22 15 
MeV unless we assume 
the fragment to be the 
result of the evapora- 


N (6.46) 


10 


Fig. 8. — Angular distri- 

bution of *Be fragments. 

A collimation in the for- L 

ward direction is indi- lo 
cated. 
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tion of a light’ nucleus. The minimum energy of an evaporated particle from 
a heavy nucleus is determined by eq. (4) and for Be should be 21 MeV. 


84. SLi fragments from stopped 7 stars. — In connection with the study 
of heavy fragments originating from high energy nucleon interaction it was 
felt advisable to analyse the occurrence of *Li fragments originating in stopped 
n- stars, the thought being that since the excitation energy is known, namely 
140 MeV, a study could be made as to the evaporation of SLi from heavy 
nuclei contrasted to the evaporation of *Li from light nuclei. In the case of 
heavy nuclei, namely silver and bromine, the excitation energy is well below 
the total binding energy of the respective nuclei and the evaporation theory 
should apply. However, in the case of light nuclei (COHN) the excitation 
energy is greater than the average binding energy of the nuclei. 

A stack of 36, (2x 4) inch, 600 um G-5 emulsion were exposed to the x 
beam of the cyclotron at the University of Chicago. The beam entered the 
stack perpendicular to the 4 in. dimension and parallel to the emulsion. 
After travelling 1.5 cm of emulsion 95 percent of the 7- mesons stopped within 
a strip of length 2 cm. The pion flux was 4:10‘ 7’s per cm? per minute and the 
stopped u- contamination was less than one percent in the portion of the 
emulsion where the 7 mesons stopped, or less than one percent of the zero 
prong stars can be attributed to = absorption. Over 22308 o stars were ana- 
lysed and the following problems emphasized. | 


84.1. Prong distribution of o stars. — Allo stars found in the stack 
were catalogued and classified according to prong number. The results are 


TABLE II. — Sigma stars from stopped = beam ecitation energy 140 MeV. Total number 
of stars: 22308. 
| | 
No ADELMAN CHESTON | MENON, | Group from | Barkow, Barkow, 
DE and | and | MUIRHEAD | MARSHAK | ALUMKAL | ALUMKAL 
| 33 JONES GOLDFORB and ROCHETI P—182 | and CHA- and CHa- 
| Prongs| (9%) (%) | (%) (aver.)|  % © | many (no.) |. many % 
0 26.8 35.2 19.4 | 28 +2 6 684 | 29.9 
1 21.5 1975 33.4 125 SE 0 7389 SO 
2 27.0 21.9 23.4 23.5+0.7 4444 19.9 
2: | 15.2 | 18.9 20.1 | 1612203 | 2 624 | LS 
4 Te A e enc sy mA 
5 1.8 1.9 == (mn 45020 175000 Wee | 
6 _ | — — a One On il == 
7 = | = | = | CA | 3 = 
8 | == —- —— | = 1 | — 
| No. of stars | 
| studied 1500 | 500 | 2500 | 4000 22 308 


. SBLì FRAGMENTS ASSOCIATED WITH NUCLEAR DISINTEGRATION 327 


tabulated in Table IT and compared with results obtained by ADELMAN and 
JONES (22), CHESTON and GOLDFORB (**), MENON, MUIRHEAD and ROCHET (ea) 
and the average given by MARSHAK (2°), The possibility of including the recoil 
nucleus as a prong is very small. For a recoil nucleus of charge Z=5 and 
mass M =10 to have a range of 2 um, which is barely visible, required one 
proton to have energy of nearly 45 MeV. The probability of observing protons 
with energy greater than 45 MeV is less than two percent and since this is 
divided between light and heavy nuclei, the actual observing of a light nucleus 
recoil is less than one percent. 


84.2. Prong distribution of stars having a ‘Li fragment. 
The frequency distribution and energy distribution of ‘Li fragments observed 
in o stars were studied. The criteria for selection of events were, 1) all SLi 
hammer tracks were assumed to have stopped in the emulsion and their co- 
linearity checked within 5°; 2) energy of measured ‘Li fragments had angle 
of dip 30° or less in the undeveloped emulsion and 3) all questionable hammer 
tracks were assumed to be due to *Li. Three strongly suspicious *B fragments 
were observed and are not included in the results. The reason for selecting 
the 30° dip criterion was to simultaneously check the isotropy of ‘Li emis- 
sion, since the solid angle subtended by the surface area of a sphere from 0° 
to 60° is equal to that subtended from 60° to 90°. A typical four prong o star 
caused by an absorbed stopped z~ is shown in the mosaic photograph in Fig. 9. 

The prong distribution of the parent star from which the SLi fragment ori- 
ginated is syown in Table IIT, and contains both the data for all *Li fragments 


Tage III. — Prong distribution of parent star of *Li fragments 
(data from over 235 000 o stars). 


| . | Frequency of SLi | 

No. of | Frequency of *hi (dip angle less Percent of 

prongs | (total) than 30°) | total 8Li 
| | : pe 
| | | | 
| 1 | 5 2 eal | 
| 9 | 51 20 | 11.4 
| 3 | 266 139 | 59.5 | 
| 4 | 113 59 | 25.3 | 
| 5 | 12 5 | DT 
| - 
Ron 447 225 | 100 


(22) F. L. ADELMAN and JONES: Science, 3, 226 (1960); F. L. ADELMAN: Phys- 
Rev., 85, 249 (1952). 

(28) W. B. Cnesron and L. J. B. GOLDFARB: Phys. Rev., 78, 683 (1950). 

(24) M. G. K. Menon, H. MUIRHEAD and 0. RocHar: Phil. Mag., 41, 583 (1950). 

(25) R. E. MARSHAK: Meson Theory (New York, 1952), p. 182. 
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Fig. 9. — Mosaic photograph of the capture of a r--meson by a nucleus resulting in 
c=) 


a four prong star, one prong of which is a 8Li fragment. 


Lal 
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observed and those that satisfy the 30° dip criterion. The fact that the emis- 
sion is isotropic for both cases is apparent and this isotropy was verified by 
an angular distribution plot. From Table IIT 
the *Li fragment favors the environment of 60] 43 
two other fragments as 60 percent of the 
hammer tracks were observed in 3 prong stars. 

The energy distribution of *Li from 3 
and 4 prong stars and from all stopped 77 30 
stars is given in Fig. 10 and appears to be 
independent of the size of the parent star. 
Since 90 percent of the *Li fragments have 
observable energy of 16 MeV or less it is On 
evident that they must have originated from 
light nuclei. This follows from the fact that 
16 MeV is the minimum energy of an evapo- 
rated particle from a heavy nucleus. 

In an attempt to check this conclusion 
further we carefully analysed the fragments 
accompanying the *Li in three prong stars. 


50 


40 


20 


10 


I 
' 
' 
' 
' 


Era ork Sees 


5 
ON se 


Fig. 10. — Energy distribution of *Li from stop- 
ped = stars. The upper curve is the distribution 
from all stars regardless of the prong number. 
The middle curve is the distribution of $Li from 
4 prong stars while the bottom curve is the distri- 
bution from 3 prong stars. Since the potential 
barrier for *Li fragments from heavy nuclei is 
about 16 MeV, it is evident that 90% of the ag © 6% 21 274 

fragments originated from light nuclei. Energy in MeV 


Taste IV. — Analysis of n interactions from three prong stars. 
ones | ou Reaction scheme Percent 
fragments | events 

D prsli 9 ma +420 > 21H+SLi+2n 7.5 
Da li 48 NC Ho an 39.0 
a, a, SLi 15 n-+0 + 2:5He4 sai 12.5 
D, ali; SLi 17 NOs Hi lion =| | re 
DA Lili DE; noie tH hai SLI I ‘ 
heavies, 5Li 32 r-—heavy 27.0 
121 100 — 


oe 22 - Il Nuovo Cimento. 
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Of the 139 events that satisfied the dip criterion, 121 could be classified 
as clean stars giving no evidence of a short recoil. The fragments were 
identified by the method discussed previously with no width corrections made 
when a fragment had a dip angle greater than 30°. The analysis is given in 
Table IV and indicates that about 73 percent of the three prong stars could 
owe their origin to the capture of a 7- meson by a light nucleus. 


4. — Conclusion. 


At the present time nuclear evaporation is the only theory from which 
relative frequency, energy spectra and angular distribution of fragments can 
be calculated. It is generally assumed that evaporation breaks down when 
the excitation energy of the nucleus exceeds the total binding energy. The 
upper limit for application of the theory to heavy nuclei (Ag and Br in emul- 
sion) is generally taken as about 800 MeV resulting in a nuclear temperature 
of about 10 MeV. Hence, any temperature greater than this, like the value 
of 19 MeV for beryllium fragments is somewhat too high to be accepted. Also 
the presence of the high energy tail accompanying the distribution and the 
contribution of the fast nucleon-nucleon knock-on stage which precedes the 
evaporation process cannot be neglected. The observation of high energy heavy 
fragments which cannot be explained on the basis of elastic collision (18) strongly 
suggests the presence of hot spots in the nucleus. It is generally accepted 
that the fair agreement between calculation and experiment in the low energy 
region is a good indication that the evaporation mechanism is dominant, but 
that at high energies another mechanism takes over. 

For low excitation energy as in the case of the 77 exposure 90 percent of 
the ‘Li were emitted from supposedly light nuclei and the application of evapo- 
ration theory for light nuclei becomes questionable. The light nuclear case 
together with the presence of a small Coulomb barrier is complicated by the 
additional processes which lead to the formation of SLi as spallation residues. 
For a nitrogen target, an average light nucleus in the emulsion of 4 equal 
to 14 and Z equal to 7, the concept of evaporation does not apply since the 
‘Li has about half the mass of the target. In other words, what is the meaning 
of the potential barrier for a nucleus that absorbs a 77 meson resulting in a 
3 prong star whose fragments are ‘Li, an « particle and a proton? Results 
predicted by the evaporation theory especially those relating average energy 
of the fragment with nuclear temperature and the Coulomb barrier are in good 
agreement with experimental results (Fig. 10). The frequency of observing a 
‘Li fragment originating from a heavy nucleus is in agreement with the 
theoretical consideration of LECOUTEUR (25). 


(26) K. J. Le CourEUR: Nuclear Reactions (Amsterdam, 1959). 
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An additional controversy arises in the study of the angular distribution 
of Li from cosmic ray stars and heavy fragments Z > 4. In Fig. 11 we have 
plotted the combined angular distribution of over 600 heavy fragments of 
Z > 4, observed in cosmic 
ray stars, 4.5 GeV pion 


stars and 400 MeV proton °° 
stars. All observed frag- I 
ments greater than lithi- 
um, regardless of whether si 
they were identified or 
not, are included in the a 
plot. The distribution also eA 
includes the 98 fragments 
identified as ‘Be from fe 
Fig. 8. At high excitation 
energy there appears to 55 
be a definite collimation of 


fragments in the forward 


direction, while for re- O 2 40 60 80 100 120 140 160 180 
Unit solid angles 


latively low excitation 

(140 MeV) the angular Fig. 11. - Angular distribution of heavy fragments for 

Z>4. All observed fragments for Z greater than 3 were 

plotted regardless of whether they were definitely iden- 
tified or not. 


distribution indicates an 
isotropic emission of ‘Li 
fragments. The anisotropy 
observed in the case of 
high energy fragments can only be explained on the basis of nuclear evapora- 
tion by assuming the fragment emission to take place during the flight of the 
highly excited rest nucleus formed after the initial nucleon cascade is com- 
pleted. The theoretical discussion of SKJEGGESTAD and SORENSEN (1), calculated 
on the assumption that the evaporating nucleus is in motion predicts a colli- 
mation of fragments in the forward direction for high energies and becomes 
nearly isotropic at energies less than 30 MeV. 

In communication with LEVI-SETTI and SKJEGGESTAD a very similar distri- 
bution of *Li fragments from stopped K stars is being obtained. From pre- 
liminary results 90 percent of the *Li fragments have energy less than 16 MeV 
and range less than 50 pm. 


The authors are very grateful to Professor MARCEL SCHEIN, of the Univer- 
sity of Chicago, for his generous assistance in regard to this problem. It is 
also a pleasure to thank Drs. ERICK LOHRMANN and MARTIN TEUCHER, of the 
University of Chicago, for many stimulating discussions. 
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APPENDIX 


The Le Couteur formula can be obtained if we consider the nucleons within 
the excited nucleus to follow a Fermi distribution of the form 


constant g(E) exp [— E/kt} 


= dE , 
exp [x] exp [— E/kt] 


(1) n(E) dEj= 


where n(E) is the number of nucleons per energy interval within the nucleus 
and 4(E) the number of levels per energy interval and is proportional to E*. 
If exp « > exp [— E/kt] and when compared to exp [— E /kt] is a slowly varying 
function of kt and can be considered constant, eq. (1) becomes 


(2) n(H) dE = C,g(E) exp |. al dH 
or 

E 
(3) n(£) dE = C, Ht exp | — A dH. 


The probability of a nucleon escaping from the nucleus is proportional to 
n(E), provided the energy of the particle Æ is greater than the height of the 
potential barrier V, and also proportional to v the velocity of the particle or 
to E?. The probability that a particle of energy between H and H+d# will 
escape is given by 


P(E)jdE = 0, Et: Ht exp |— A dH 
or 
4 P(E)d Tepe | 
(4) P(E)dE = C, E expj|— Sd, 


where 7 is in units of kt and usually expressed in MeV. If we do not permit 
leakage through the barrier, the smallest value of £ for a nucleon that escapes 
will be V, the height of the Coulomb potential barrier. Normalizing eq. (4) 
to obtain 0, we have 


fe exp 
Vv 


since the observed energy of a nucleon outside the nucleus is (E —V) we have 


E 
— pe = 4 


(co) 


C; i E exp 


Ve 


Tulcea lye 
TRS 
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rom which OC, = 1/22 cor 


Ve | 
5 E) dE = - ex rs 
(5) P(E)dE qa °XP | T |a 


which is Le Couteur’s formula. 


RIASSUNTO (*) 


Questo è un rapporto preliminare sulla emissione di frammenti pesanti nella disin- 
tegrazione nucleare. Si presuppone che la probabilità di emissione di frammenti pesanti 
dipenda dalla temperatura nucleare (7) e dalla barriera di potenziale (V) del nucleo 
eccitato. Lo scopo di questa ricerca è di studiare la dipendenza dello spettro differen- 
ziale dell’energia dall’energia di eccitazione. I dati raccolti provenivano da quattro 
pile di emulsioni Ilford G-5, esposte a: (a) un fascio di pioni frenato; (è) un fascio di 
protoni di 400 MeV; (c) un fascio di pioni di (4--5) GeV e (d) raggi cosmici con espo- 
sizione nel Texas, latitudine magnetica 41° N, e sull’isola di Guam all’equatore magne- 
tico. In questo scritto si sottolinea l’emissione di 8Li come frammento pesante. I risul- 
tati dell’esperimento indicano che per elevata energia di eccitazione si ha una colli- 
mazione in avanti dei frammenti di litio e berillio, che può spiegarsi sulla base della 
teoria dell’evaporazione, supponendo che l’emissione dei frammenti abbia luogo durante 
il volo di un nucleo fortemente eccitato. Nel caso di bassa energia di eccitazione la 
distribuzione angolare dei frammenti di 8Li è isotropica e la distribuzione dell’energia 
indica che la maggior parte proviene dall’assorbimento di un mesone zx in un nucleo 
leggero. 


(*) Traduzione a cura della Redazione. 
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On the Anisotropy of the Cosmic Radiation. 


E. W. KELLERMANN and M. S. ISLAM 


Department of Physics, The University of Leeds - Leeds 


(ricevuto il 14 Aprile 1960) 


Summary. —— An attempt has been made to reproduce some results of 
McCusker and his coworkers who had reported diurnal effects of cosmic 
ray showers. The design of the apparatus follows closely that of these 
workers, but allows more detailed information about each recorded 
event to be obtained. The main categories of observed events, including 
those examined by McCusker et al., show no significant diurnal variation. 
A possible sidereal diurnal variation of a particular group of events is 
critically examined, and it is concluded that, in spite of strong statistical 
weight, the implications of accepting the variation as physically real do 
not at present allow it to be regarded as such. 


Many attempts have been made to determine whether an anisotropy of 
the cosmic radiation, apart from that imposed in the immediate neighbourhood 
of the earth, can be detected. Although for obvious reasons one would regard 
as most promising the study of the products of very energetic primaries, a 
number of workers who report some indication of anisotropy have been con- 
cerned with effects due to primaries of much lower energy. Notably, McCuskER 
and his co-workers (1) have presented evidence for the detection of a sidereal 
time variation of large amplitude, and hence a directional anisotropy of inci- 
dent primaries, with a relatively simple counter arrangement sensitive to 
penetrating showers. Other workers, too, have reported similar, if smaller 


. B. A. McCusker, J. G. Darpis and B. G. Witson: Proc. Phys. Soc., A 68, 
i. B. A. MCCUsKER: Nuovo Cimento, 2, 1340 (055) 

(5) C. B. A. McCuskEer and R. J. Rein: The Oxford Conference on Extensive Showers 
(Harwell, 1956), p. 58. 
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effects, and the great difficulty in seeing any reasonable interpretation of these 
observations must be set against the fact that considerable agreement is to 
be found in such work regarding the phase of the possible sidereal variations. 


The work reported here was prompted by 


the need for investigation of the large side-: 


real variation reported by MCCUSKER and of 
the very large variation in solar time repor- 
ted by him (3) by similarly simple counter 
arrangements. Data for 14 months are now 
available, and so a stage has been reached 
when significant statistical improvement will 
be slow; they are perhaps now of sufficient 
interest to justify this report. 

The experiment is being carried out at 
Leeds. Following MCCUSKER closely we use 
a penetrating shower detector of three trays 
of Geiger tubes (Fig. 1a). Each tray contains 
9 counters of area 240 cm?, and they are re- 
spectively shielded vertically by 30, 40 and 
50 cm lead while there is 15 cm lead shield- 


50 cm 


Fig. la. — The shielded trays S. 


ing at the sides. This unit is placed at the centre of a triangle (Fig. 1b) formed 
by three unshielded trays of counters each containing two of area 240 cm? 
and two of 35 cm? area. The discharge of every counter of the array is separa- 


Fig. 1b. — The geometry of the shielded 
trays S, and the unshielded trays A, B 
and C. 


tely displayed, and the form of a pene- 
trating event together with an estimate 
of the accompanying electron shower 
density may be obtained. The selec- 
tion was made for penetrating showers, 
involving the discharge of at least two 
counters in each of the three shielded 
trays. 

The hodoscope arrangement, enabling 
the discharge pattern to be examined, 
allowed different kinds of events to 
be distinguished. Some of these could 
be compared with observations of other 
workers, affording a useful check on 
overall operation; it was felt, too, that 
any anisotropy found to occur in events 


of any but extreme energy might belong not to a broad grouping of events, 
but rather to some small, but perhaps recognisable component. Accordingly, 
the rates of the separated components making up the total observed group of 


+ 
1 
à] 
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( penetrating showers » were analysed both in solar and in sidereal time, tested 
for fit to isotropic incidence, and only where this failed subjected to Fourier 
analysis over a solar or sidereal day. 

The period of operation reported here was of about 10000 hours (between 
February 1958 and June 1959); in it the following events were distinguished: 


1) Knock-on events, 4047 in number. 

These, arising from the successive elastic collisions of muons near to each 
shielded tray were selected from the hodoscope record by inspection. The 
selected number was later compared with a calculated rate based on the 
established muon flux. Agreement was reasonable, and the separation of these 
events is probably good. 


2) P,-showers, 18 482. 
Penetrating «local» showers, discharging at least two counters in each of 
the three shielded layers, by their geometry probably not of type 1) above, 
and not accompanied by any discharges in the unshielded trays. 


3) P,-showers, 2 455. 
Penetrating showers accompanied by at least one counter discharge in the 
unshielded array. 


The pattern of discharges of shielded counters within the central lead shield, 
represented in Figs. 2(a), 2(b) by the 
total number, 6 < n < 27 of these coun- 
ters discharged, shows a significant dif- 
ference between P, and P, events; the 
latter yield a very much higher average 
number of discharged counters, and 
therefore probably refer to more ener- 
getic events. 

In the P,-events (Fig. 2(b)) there is 
a tendency for the frequency-number 
relation to become constant for a high 
number of discharged counters. This 
may represent saturation, but might also 
be the result of superposition of two 
effects: the P,-events have according- 
ly been divided into the categories 


Relative frequency 


Fig. 2. — The frequency of numbers of 


| | shielded counter discharges for: a) penetra- 
10 ting local showers P,: b) penetrating ex- 
tensive showers P,. 


iG. 2 
Counter Discharges 
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P,(> 19) and P,(< 18) depending on whether the count does or does not lead 
to the discharge of at least 19 of the shielded counters. They are the events 


(3a) ee 19)5 471 
and 
(35) PIS) 40934. 


A final subdivision of the P,(>19) events was made on the basis of the 
density of accompanying electronic shower particles as estimated by the de- 
tailed discharge pattern in the unshielded counter trays surrounding the central 
unit (Fig. 3). The figure leads one to make a separation at an electron density 
of about 50 particles per square metre, and the final subdivision there leads 
to the categories 


(3a, a) Ria SS 0 ne) ANAL 
and 
(3a, db) Pee Lo D0 Mia) eS. 


In order to acquire some information comparable to the results published 
by McCusker (*) referring to showers of high particle density, we selected, for 
some six months, also 
events triggered only by 
the small counters in the 190 
unshielded trays, where 
at least one small counter 
in each of the trays was 
discharged. We estimated 
the average density of 
these events as larger than 
200 particles m-?. About 
730 of these high density, 
H.D., events were regis- 
tered. 

All of the groups listed 
above were, as has been ® 
indicated, tested for com- Ro Sur VE 
patibility with isotropic Density 
incidence in both solar and Fig. 3. - The density distributions of the soft component 
sidereal time and were accompanying the P,(> 19) events. 


T 


oO 
(æ) 
el 


(©) 
T 


Relative frequency 


(4) C. B. A. McCusker, D. E. PAGE and R. A. REID: Phys. Rev., 89, 1172 (1959) 
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subjected to harmonic analysis where the level of compatibility was low. The 
raw data were not corrected for pressure or temperature since such correction 
would only be possible for recognizable components, and such identification 
is not possible as yet: moreover, our statistics are not yet sufficient to derive 
a pressure or temperature coefficient from our data. 

The results of the analysis are given in Tables I and II; the compatibility 
of the data with isotropy is measured by a y?-test for constant intensity through 
a solar and sidereal day respectively, and the probability of fit, P, is given in 
column 3 of each table; the usual convention of examining the significance of 
the value of P at the 5%—and at the 1%-—level has been adopted. 


Tape I. — Results of solar time analysis. Total running time 9803 h solar time. 


Ist harmonic | 224 harmonie 


Type of shower Ce ea ES Po atene] È Reano 4 
: events with st. line (0%)! Ampli-| Phase | (%) ‘Ampli- Phase) (%) 
| | tude ee!) | tude | (h) 
| | | 
1) Knock-on 4047 | ? 3064222.7/2.80- Ls NOTES 
2) P, 18482 | Yess| 390 | — Cul SA 
SP, 2455-| Yes | 36 | — = | 2 
Baie PL (S19) 471 | Yes | 63 = = aa IESPA = o 
SEP (LS) LOST p45 NS Sa ii e 
Bd GP, (19, = 50 m=*)) 191) Yes | B00) e = =| COLE 
3a, b) | Py (> 19, < 50 m—) 280 ves EL A 3 a = | a 
4) | II, (six months) 729 Yes | 70 == Le == Sa DI 


l'ABLE II. — Results of sidereal time analysis. Total running time: 9830 h sidereal time. 


RER Ist harmonic 2™¢ harmonie 

Type of shower | No. of | Straight P sn —| P| È P 
events, line | (0, | Ampli- | Phase h (%)| Ampli- |Phaseh (%) 

tude” Gi) Gatti ITU Mir 
(1) Knock-on | 4047 | Yes |87 = | = — I a 3 
(2) TE 118482 | Yes |33 3 al econo di EA 
Cie wes | 2455 | Yes |50 nn ee fra RU IR SERRES 
(34) Pal ee 471) No | 0.05)20.7 4 7.4/21.6+1.4| 7 |20.7+ 6.5|!29t28| 5 
(3b) Pr (<18) 1984 | Yes 10 : | == ro = 

x D TOR ~2 7 no le ‘ = 5 | 
CE RIO 50 De 191 No 0.2 |30.9-+11.3}21.5 41.3] 2 |35.1411.9| 556507 | 16 
(3a, b) | Pz (219, <50 m >?) 280 | Yes | 8.6 = es SS = = = 

| (4) \ H,, (six months) 729 | Yes |40 x | | | 
| | = = = 


À & x ‘ > ‘ tn bd D D © Q A È 

The solar and sidereal analyses may be commented upon separately; the 

lr . . Ù ‘ 
only indication of a solar effect concerns knock-on events, and even here the 
Significance of the variation is not particularly high. Neither amplitude nor 
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phase of the possible effect agrees with that established for the general meson 
component, and there is no reason to postulate a real effect. It should be 
mentioned here that the combined groups of knock-on and P, events roughly 
correspond to the « penetrating 


local showers » examined by Rates 
McCusker and WILSON (1). 

The P, showers as a whole cor- der 
respond to the penetrating exten- ae ° > 6 
sive showers for which McCuskER SS |° Py n : ° o 
and Rem (*) have reported a 2 œL* ‘ ; 
sidereal variation. Since our re- DI è Me Cusker 
sults have statistics comparable Ol }- » Leeds | 
to those of MeCusker and Reid, 0 6 = a cia 
both their results and ours have (sidereal time) 
been plotted together in Fig. 4. Fig. 4. — The rates of penetrating extensive 
The figure illustrates the findings showers versus sidereal time, as measured by 


: ’ ; crea SN à 3 Ta 
in Table I, namely that, as for the McCusker (*) and in Leeds. 


local showers, our results confirm 
neither the phases nor the amplitudes found by these workers. 

The possible sidereal effect concerns the results grouped as P,(> 19) and, 
in particular, those related to high density electron showers, P,(> 19, > 50m?) 
rather than those related to lower shower densities, P,(>19, < 50m?) 
although one must bear in mind that the subdivisions leading to the isolation 
of the P,(>19) events are arbitrary. The rate of P,(> 19, > 50 m?) events 
as a function of sidereal time is plotted in four-hourly and in two-hourly inter- 
vals in Fig. 5. The continuous graph is the calculated second harmonic, the 
dotted line is the first harmonic. For comparison we have plotted, too, the 
rates in solar time and «antisidereal » time. 

The «antisidereal» time scale, first used by FARLEY and STOREY (5), has 
no physical meaning and is arrived at in a purely formal way. FARLEY and 
STOREY have pointed out that as the sidereal time can be regarded formally 
as a modulation of the solar time, so an «antisidereal» time can be obtained 
by a modulation of solar time opposite in sign. Since this time scale is devoid 
of physical meaning no variation should be expected and none, as in solar time 
representation, is in fact found. 

The evidence presented in favour of the reality of sidereal variation must 
be weighed carefully in conjunction with the following arguments: 

Firstly: it is important to establish whether the events in question, the 
P (219, >50m?) are a selection from extensive air showers, (EAS). If 


(6) F. J. M. FarLEY and J. R. STOREY: Proc. Phys. Soc., A 67, 996 (1954). 
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: a x 
they are, then in view of all the accumulated knowledge on EAS, the time 
variation reported here is most likely a statistical fluctuation. Certainly, the 
rate of the events, — 0.5 per day, is about 5% of the rate of showers of den- 


ESSE | | 1 | L 
4 ROBE AUOT du 


(antisidereal time) 
0.02! : ri [eo $e vena 
l | | | ] Il = = 
2 16 20 0 £ hours 
es : ; (solar time) 
20.02) $ 
=) 
È r] 
S [SS | ] ; je cs Al 
S 4 8 12 16 20 0 EN 
È (sidereal time) 
"50. 04h 
a ® 
0.021 + } ; 
| o) Eee e [ i L a 
4 8 12 1 20) 
0.04+- _ hours 


(sidereal time) 


2 | 
4 8 12 16 20 0 4 hours 
(sidereal time) 


Fig. 5. - Pz (> 19, > 50 m-?) events plotted in antisidereal and solar time, and in two- 
hourly and four-hourly intervals in sidereal time. 


sity of (100200) particles m-* expected, according to the known density spec- 
trum of HAS, to fall on the pile. Such a percentage is not unreasonably high, 
and the identification of these particular events with a selection from EAS 
is therefore not ruled out. The energy one would have to associate with the 
primary giving rise to the events would judged by the rate, be of the order 
of (10151016) eV. 

A further argument against the reality of a sidereal variation is, of course, 
the amplitude (30%) of the effect which is unplausibly large. 


The arguments against the existence of the effect are at least modified by 
the following: 
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First, a systematic error is less likely to occur in a sidereal than in a solar 
analysis. Here, the analysis in solar and in «antisidereal» time shows the 
expected absence of the effect. There is agreement too, with the results of 
other workers, but not with McCuskEr, on the absence of a large effect in 
the group comprising all penetrating extensive showers. 

Secondly, the main maximum of the variation reported here (21.5 h) co- 
incides with the phase for which other workers (see for example GREISEN (°), 
GALBRAITH (7)) have found—albeit not statistically significant—indications. 

The size of the variation reported here relates to a separated component 
and so does not correspond to the total cosmic ray flux investigated by pre- 
vious workers. In particular, normal observation of extensive air showers 
would, even if this effect were real, only show it in an extremely diluted form. 

However, even if one does assume that one is dealing here not with EAS 
but with a correctly separated and independent component of the cosmic 
radiation, there are still great difficulties in accepting the effect as real: such 
an independent component would have to be due to a neutral particle in order 
to account for its arrival in a preferred direction and thus to have escaped the 
diffusion process. At the latitude of Leeds (53.8° N) and at 21.5 h the appa- 
ratus faces along a direction near that of the spiral arm of the galaxy. One 
can see that a neutron of 101° eV, though not any known particle of shorter 
lifetime, could have a range comparable with galactic distances, but an attempt 
to connect the effect with the existance of neutrons leads to further impli- 
cations which are difficult to reconcile. 

A strong directional effect implies that if neutrons arose by a charge- 
exchange of protons after acceleration, this must take place in the neigh- 
bourhood of sources rather than after complete diffusion of the protons. Since 
some concentration of matter is essential for an exchange process, this is in 
a sense satisfactory, but the implicationt that the maximum intensity should 
be observed in the direction of maximum concentration of sources does not 
fit with the observed direction. A further difficulty which still remains is how 
such a neutron could give rise to effects which are experimentally distinguishable 
from the main EAS. 

All these considerations, therefore, lead us to the conclusion that the reality 
of the effect remains doubtful in spite of its statistical weight. Certainly, our 
investigation does not confirm the existence of large time variations as re- 
ported by McCusker in the case of the total flux of cosmic ray air showers. 
We are continuing our experiments to examine the small component of cosmic 
ray showers characterized by high multiplicity under lead. 


(5) K. GREISEN: Progr. in Cosmic Ray Phys. (ed. J. G. Witson), vol. 3 (Amsterdam 
1956), p. 124-125. 
(7) W. GALBRAITH: Exlensive Air Showers (London, 1958), p. 162-163. 
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RIASSUNTO (*) 


Si è fatto un tentativo di riprodurre alcuni risultati ottenuti da MeCusker e colla- 
boratori che hanno descritto effetti diurni degli sciami di raggi cosmici. La struttura 
dell’apparecchio segue da vicino quello di questi autori, ma permette di ottenere infor- 
mazioni più dettagliate su ciascun evento registrato. Le principali categorie di eventi 
osservati, comprese quelle esaminate da McCusker et al., non presentano significative 
variazioni diurne. Una possibile variazione diurna siderea di un particolare gruppo 
di eventi viene esaminata criticamente, e si conclude che, malgrado il forte peso sta- 
tistico, le implicazioni derivanti dall’accettare la variazione come fisicamente reale 
non permettono che essa sia considerata tale. 


(*) Traduzione a cura della Redazione. 
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Radio-Frequency Discharges in Magnetic Fields. 


F. Bist and B. DE MICHELIS 


Istituto di Fisica del Politecnico - Milano 


(ricevuto il 14 Aprile 1960) 


Summary. — Measurements on radio-frequency discharges in argon have 
been performed by using a magnetic field normal to the r.f. electric field. 
The power absorbed by the discharge has been determined as a function 
of the pressure inside the chamber and the magnetic fields: different 
curves have been obtained by increasing or decreasing the pressure. 
Different geometrical arrangements have been used. The qualitative 
behaviour of the discharge has been interpreted by using an average 
electron picture. 


1. — Introduction. 


Many experiments have been made on electrical discharges in a.c. fields. 
In the last few years some results on r.f. discharges in magnetic fields have 
been reported and discussed (1%). The essential features of the processes in- 
volved are sufficiently known and numerical data can be obtained in some 
simple configurations by applying the Boltzmann theory that takes into ac- 
count the electron energy distribution inside the discharge. However due to 
the large number of discharge parameters involved, some processes have not 
been adequately studied. 

In our measurements we used a magnetic field normal to the r.f. electric 
field. This arrangement is responsible for some of the observed phenomena, 
namely the appearance of a high intensity discharge at low pressure, because 
the electron oscillation amplitude in the direction of the electric field is strongly 
constricted. We have used different geometrical arrangements of the electrodes. 
(1) S. C. Brown: Basie Data of Plasma Physics (London, 1959). 

(2) L. GoLpsrEIN in: Electronies and Electron Physics, 7, 399 (1955). 
(3) S. C. Brown in: Rend. del III Congresso Internazionale sui fenomeni di toniz- 
1 


zazione nei gas (Venezia, 1957), p. 169. 
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Some preliminary measurements were made on a discharge between parallel 
plates, so as to have the advantage of a simpler configuration and a better 
analytical elaboration of the obtained data, but the large drift velocity of the 
electrons in the direction normal to both electric and magnetic fields prevents 
the high intensity discharge at low pressure. The phenomenon is encountered 
in discharge between coaxial cylinders: we have used two different sizes of 
internal electrode. The results obtained with the smaller electrode are the 
most significant and most measurements were made in this configuration. We 
have also studied the discharge with a larger internal electrode, so to approach 
the condition for a constant electric field, as postulated in the development 
of the theory. The results are in agreement with the theory, but due to re- 
duced dimensions, only a few discharge states are present. With radial con- 
figuration and smaller internal electrode, the discharge shows a different be- 
havior for curves obtained by increasing or decreasing the pressure: in the 
latter case the experiments do not only show the maintenance of discharge 
at lower pressure than that of threshold, but also the appearance of new 
States of discharge not observed when pressure was increased. Using an average 
electron picture a qualitative interpretation of phenomena is given (14). 


2. — Apparatus. 


Our measurements have been performed in a discharge chamber as in Fig. 1. 
The discharge was excited inside the cylindrical tube of pyrex by the r.f. field 


q ae 


3 |E —D Pirani gauge 


—-To Argon reservoir 


0 9 10cm 


To vacuum pumps 


Fig. 1. — Discharge chamber: a) plexiglass window; b) coils; c) axial graphite electrode; 
d) steel mesh; e) pirex cylinder. 


(*) B. Lass, W. P. Artis and S. C. Brown: Journ. Appl. Phys., 21, 1297 (1950). 
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applied between the axial internal electrode of graphite and a steel mesh around 
the cylinder: the direction of the magnetic field was the same as the axis of the 
cylinder. The gas used (argon at 99.9%) was introduced by a needle valve 
and the pressure was read on a Pirani gauge: measurements were made in the 
pressure range between 1 and 500 um Hg. 

The principal geometric and electric data are reported in Table I. 


TABLE I. 


Chamber length 


Oscillator power . 
Output voltage. . 


3. — Results. 


Chamber diameter . IS (SONO ES 
Diameter of internal electrode: small. . 


Electric field frequency . . 


Maximum magnetic induction . . 


large. . 


120 mm 
50 mm 
2 mm 
34 mm 
1 MHz 
100 W 
500 V 
600 G 


In Figs. 2 and 3 are reported the results of our experiments. The power 
absorbed by the discharge was determined by measuring the variation of the 


oscillator plate current as a func- 
tion of the pressure inside the 
chamber. The voltage between 
the electrodes, measured by an 
electronic voltmeter, did not 
change greatly during the meas- 
urements: the small varia- 
tions were compensated by 
means of a 1000 pF capaci- 
tance inserted in the discharge 
circuit. In every position the 
oscillator was syntonized. The 


Fig. 2. — Plate current increase as 
a function of the pressure for dif- 
ferent magnetic fields. The curves 
were obtained by increasing the 
pressure. In curves a), b), c) and d, 
the small electrode was used; in 
curve e) the large electrode was 
used. 
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curves in Figs. 2 and 3 are reported for different values of the magnetic field 
and were obtained respectively by increasing (Fig. 2)or by decreasing (Fig. 3) 
the pressure; a small internal electrode was used. 


a) Small electrode-increasing pressure (Fig. 2). — At magnetic fields under 
300 gauss, the curves show two different states of the discharges; a low inten- 
sity discharge excited at low pres- 


ad sure is suddenly changed, at higher 
307 pressure, into à high intensity dis- 


Oy D FT 8=250gauss charge, the absorbed power increas- 
al | ing to a higher level. The varia- 


20 | tion is also visible due to à sharp 
15- J 8=380 gauss increase in brigthness. 

04 == = = 

304 At higher magnetic field the 


154 EE B=450 gauss threshold pressure is lower and 
07 


the initial discharge appears as a 


40-4 high intensity process; however 
the intensity decreases sharply and 
a low intensity discharge is main- 
tained at higher pressure. By 


207 8=500 gauss 


405 increasing the pressure, a new thre- 
204 8=580 gauss shold is found for a high intensity 
al discharge. The threshold pressure 


of the first state seems to be 

indipendent from the value of the 

magnetic field. 

, cr = 5 oe on i the Di Magnetic ee 
curves the initial discharge, dis- 

playing a peak, was always pre- 

ceded by a very low intensity 


8=650 gauss 


Fig. 3. — Plate current increase as a function 

of the pressure for different magnetic fields. 

The small electrode was used. The curves : : 

were obtained by decreasing the pressure. discharge excited at pressures well 

below the peak threshold. This dis- 

charge state, not shown in Fig. 2, is 

visible as a feeble violet discharge and the correspondent plate current increase 

of the oscillator is quite negligible. At low magnetic fields the very low inten- 
sity discharge disappears. 

The data are not well reproducible; the general behavior of the curves 
does not change, but the values of the threshold pressures can vary with dis- 
charge conditions. In particular when the pressure of the residual gas exceeds 
10~*mm Hg, it seems to modify the value of minimum pressures required to 
excite the different states of discharge and to favour the appearance of the 
very low intensity discharge. 
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b) Small electrode-decreasing pressure (Fig. 3). — In the curve of Fig. 3 
the discharge is studied by slowly decreasing pressure. At high magnetic fields 
a new state of discharge appears and the high pressure, high intensity discharge 
decreases to the low intensity region by two steps. The high intensity, low 
pressure region is still visible and the quenching of the discharge is displaced 
towards lower pressures. At lower magnetic fields the high intensity, low pres- 
sure-discharge is observed over a more restricted region and only three states 
of discharge are present throughout the pressure range. 


¢) Large electrode (Fig. 2e). — The results are reported for a magnetic 
induction of 500 gauss: the curves obtained by increasing or decreasing pres- 
sure are quite similar. A net increase of the plate current at low pressure is 
well visible: then the discharge decreases and no other states are visible in the 
studied pressure range. Different values of magnetic induction do not change 
the process: under 300 gauss the discharge is no longer observed. 


d) parallel plates. — Some measurements were made for a discharge be- 
tween parallel plates in a gap of 12mm. A discharge was observed only well 
above 500 um Hg: the discharge is not characterized by a high intensity ini- 
tial state. 


4. — Discussion. 


A detailed description of the diffusion controlled breakdown and the phe- 
nomena accompanying low pressure discharges can be given only by taking 
into account the electron energy distribution. A physical interpretation of the 
controlling mechanism can be given more simply by referring to a simple re- 
presentative electron: however the results are not valuable for the determi- 
nation of the threshold pressures, because the radial electric field causes a 
strong asymmetry in the electron energy distribution. 


a) Diffusion processes. — To describe the breakdown some fundamental 
variables must be considered, the electric field intensity £, the ionization po- 
tential V;, the diffusion length A describing the discharge chamber, the angular 
frequency © of the exciting field and the pressure p (°). 

The simplest breakdown condition is the balance between the loss of electrons 
by diffusion to the walls and the ionization rate. If the average energy an 
electron gains in a mean free path 4 is w., the total energy gain before reaching, 
the wall, in absence of anelastic collisions, is w,(3/2/A?). 


(5) S. C. Brown and A. D. MacDonaLp: Phys. Rev., 76, 1629 (1949). 
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The balance condition requires that before the electron strikes the wall 
it must gain the energy required for ionization: 


(1) Vue = = 6V;. 


The factor y, less than one, takes into account the fact that electrons can lose 
energy in collisions by excitation or processes other than ionization. The 
energy u,, depending on collision frequency », can be obtained through the- 
kinetic theory and is given by: 


e2E? 
Ur or 
m(v? + w;) 
when a magnetic field B is present and © € @ = eB/m. 
A can be calculated from the diffusion equation 


D VN(r, 2) +v;,N(r, 2) = 0, 


where D is the diffusion coefficient for electrons, N is the electron density de- 
pending on the radial and the axial variables, r and 2, and », is the net pro- 
duction rate per electron. By imposing the condition that N=0 on the 
walls of the discharge chamber, the expression for A results the sum of two 
terms, corresponding to diffusion lengths A, and A, in radial and axial di- 
rection: 


In our calculations we used a mean value of the field Z: we think that the 
error introduced by this approximation does not modify the validity of the 
qualitative conclusions. 


Under these conditions the radial solution of the diffusion equation is 


N= oe | Si | LA, (1) 
Ve, SU 


where J and Y, are the zero Bessel functions of first and second kind. The 


application of boundary conditions allows a numerical evaluation of A, which 
for our vessel results to be: 


A,~0.6 em small internal electrode , 


A,~0.2em large internal electrode . 
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The axial solution of equation is 
N (2) = € cos (i) : 


The buoundary clndition gives A, = 3.8 em. 

The presence of a magnetic field in the axial direction strongly reduces the 
radial diffusion of electrons: the effect can be evaluated by using a diffusion 
length A, corrected by the presence of the magnetic field so that 


1 v? al il 
ren on: 
At vo A, A 
The validity of the diffusion theory is assured when the mean free path 
is shorter than the diffusion length: therefore equation (1) gives the breakdown 
equilibrium only wen A< A. 
In Fig. 4 and 5 1= VeV,2/3yu,, A» and 4,, the ionization mean free path, 
are reported as functions of pressure p for two different values of the magnetic 


cm 
104 


1044 


P jm Hg 
+ 


| 
1 10 LC | 10 10° 
(a) (b) 


Fig. 4. —1, À;, Ay, L/2 and d as a function of the pressure for B—500 G. Small electrode. 
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field B= 250 gauss and B = 500 gauss, assuming a mean diffusion velocity of 
electrons w=3-108 cm/s, corresponding to an energy of about 25 eV, and 
y= 0.02. 


P pm Hg 


T To 
; 10 ioe 10° 10° 


Fig. 5. — 1, 2;, A,, L/2 and d as a function of the pressure for B=250 G. Small electrode. 


The values of y and w are only representative and have been chosen in a 
range of reasonable values so as to have a rough agreement with experimental 
data: however the values of both w and y can vary in a fairly large range of 
values without compromising the validity of the qualitative interpretation of 
the phenomena. In particular, the mean electron energy used in calculations 
is largely superior to the effective mean energy of electrons in the discharge; 
such an approximation is justified by the fact that, due to the presence of a 
radial electric field, the discharge starts along the axis where the electric field 
is Maximum and is mainly maintained by the contribution of the more ener- 
getic axial electrons. 

In Fig. 6 are reported the same quantities as in Fig. 4 and 5 for the dis- 
charge chamber with the large internal electrode, at B= 500 gauss. 
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cm 


P pm Hg 


102} | n T = 
1 10 10? 10° 10° 
Fig. 6. —1, A;, A, L/2anddasa function of the pressure for B=500 G. Large electrode. 


b) Oscillation amplitude limit. - When the amplitude of the electron 
oscillation in the electric field is sufficiently high the electrons can travel across 
the chamber and be lost at the walls. As a limiting case, all the electrons will 
hit the walls when the oscillation amplitude is equal to one half of the electrode 
separation length: at this limit the electron loss to the walls is greatly in- 


creased. 


The motion equation of an electron in a plane normal to the 2 axis is: 


dv 


m— (r, 0) + vmv(r, 0) = — eE(r) — e[v(r, 0)\B]- 


In the approximation E = const and @, <vy/r well verified for magnetic fields 
higher than 200 gauss, the radial solution of motion equations, for an electron 
of initial co-ordinates r=r,, 0 = 0, is: 


er : ’ , 
g = ——,, [ C1 (a Sin @t — y cos ot) + (0,(a, cos ot + » sin œt) | + 
Tao 
Aw-+ By) . Bw — Av) Ay— Bo , Civ — Co 
be ) sin wt + Lies cos wt + sto 
OW, OW, WOW, DL Wy 
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where A and B are constants depending on discharge parameters and 0, and C; 
are evaluated by the initial conditions 0,=0 and v, = 0 for t= 0. 

Omitting the transitory solution, unimportant for o<«@, the oscillation 
amplitude is: 


(2) eae Vee ee 
, mo | Avo? + (1° + wm; — w?)? 


and has a maximum for v= ©. 

In Fig. 4, 5 and 6 d is reported as a function of p in the same conditions 
as above; half of the interelectrode length, L/2, is also reported. 

The physical significance of the oscillation amplitude curves is clear by 
calculating the path of electrons in crossed magnetic and electric fields. In a 
plane normal to the 2 axis the trajectories at low pressure are circles with center 
on the + axis. When the pressure is increased, the collisions create a radial 
component of velocity and decrease the tangential drift velocity; the trajectories 
then become spirals. As pressure increases, the electrons cannot do a com- 
plete circuitation at high pressure, the motion is almost rectified; the length 
of oscillation decreases and the effect of the magnetic field is negligible. 


c) Interpretation of the experimental curves. — From the curves in Fig. 4, 
5 and 6 a qualitative interpretation of the experimental curves can be tried. 
It must be noted that only a rough interpretation can be given for three reasons, 
at least: 1) the negligibility of the true electron energy distribution; 2) the 
approximation used in the expression of the radial electric field; 3) the uncer- 
tainty about the form of the electric field in the plasma that affects the cal- 
culation of the thresholds of the different states when the discharge is present. 
A net difference can be seen in the experimental curves obtained with the 
small internal electrode by increasing or decreasing the pressure: in the latter 
case new states of the discharge can be found. In general a low intensity dis- 
charge seems to be maintained in states normally unstable, when it follows 
a high intensity discharge: the phenomenon is probably due to the contribu- 
tion of secondary electrons emitted from hot electrodes. 


x) Small electrode-increasing pressure. — For magnetic fields 
strong enough to bend the path of the electrons bringing them to hit the axial 
electrode, a discharge, probably maintained by a secondary electron resonance 
process (°), can be excited at low pressures, as shown in Fig. 2. The discharge 
intensity is high because the oscillation amplitude is shorter than the chamber 
dimensions and the loss of electrons to the walls is Strongly reduced. 


By increasing the pressure the discharge intensity decreases, because the 


(9) E. W. Gitt and A. von ENGEL: Proc. Roy. Soc., A169, 446 (1948). 
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breakdown resonance conditions are no more satisfied, the loss of electrons to 
the walls is higher and is not balanced by ionization inside the gas. However, 
the discharge is maintained probably by electrons extracted from the walls; 
the internal graphite electrode is heated to a bright red, during the discharge. 
When the loss of electrons is balanced by ionization in the gas, a high inten- 
sity discharge is again excited at pressures corresponding to breakdown. dif- 
fusion threshold (at 1 = A,). 


B) Small electrode-increasing pressure. — In the curves of 
Fig. 3 the high intensity discharged at high pressure decreases abruptly at b) 
when the diffusion equilibrium is no more maintained. The discharge remains 
in an excited state until the ionization mean free path becomes longer than 
diffusion length: under this threshold a) another sharp decrease in intensity 
allows the maintenance only of a low intensity discharge. A new high intensity 
discharge is present at lower pressure when the oscillation amplitude becomes 
shorter than chamber dimensions and a resonance breakdown mechanism is 
possible; the loss of electrons is also strongly reduced. 

An interesting feature of the experimental curves is the behaviour of the 
two thresholds, indicated as (a) and (b) in Fig. 3, when the magnetic field is 
varied: the values of such thresholds cannot be obtain by the average electron 
theory developed above, but the behaviour can be qualitatively explained. 

From the curves in Fig. 5 and 6 it can be seen that the threshold (a) varies 
between two values, given by the intersection of A; with the asymptotic values 
of A, and reaches a constant value for high magnetic fields. The threshold (b) 
represented by the intersection of / with A, increases with the magnetic field, 
at least in the studied region: by lowering the magnetic field the two thresholds 
become identical and one state of discharge is observed. These conclusions 
are in agreement with the results. 

The first state, in the resonance breakdown region, shows a constant 
quenching pressure and covers a large pressure range at increasing magnetic 
field. This last result can be justified, as is shown from eq. (2), by the fact 
that, due to reduced electron diffusion at high magnetic fields, the intersection 
between oscillation amplitude and interelectrode distance L/2 moves towards 
higher pressures. 


y) Large electrode. — A high intensity, low pressure discharge can 
be observed, the oscillation amplitude being shorter than the interelectrode 
distance: a second state of discharge with a sharp threshold is not visible for 
pressures below 1mm Hg, because, as shown in Fig. 6, there is no inter- 
section between A, and J and no diffusion equilibrium can be reached. No 
metastable state of the discharge in the decreasing pressure curves can be 
observed, because these are possible only when the metastable state follows 
a higher intensity stable state. 
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6) Parallel plates. — No high intensity-low pressure states were 
observed. In this geometrical arrangement the electron drift velocity, normal 
to both the electric and the magnetic field, is directed towards the walls of the 
discharge camber and causes a high electron loss. 

The drift velocity v, is much greater than the velocity in the direction of 
the electric field, at least in high magnetic fields: v,/v = @,/o. 

So the electron motion in the drift direction becomes important in discharge 
processes that are related mainly to its characteristics. The oscillation ampli- 
tude of the drift motion is much greater than the chamber dimensions at low 
pressure: as a function of the pressure it does not shown a maximum but 
decrease continuosly. These conclusions explain the absence of the pheno- 
mena observed in other geometrical arrangements. 


Thanks are due to Prof. G. BoLLA and Prof. E. GATTI for their kind inte- 
rest and helpful criticism. 


RIASS UNTO 


Sono state studiate scariche a radio-frequenza in argon, usando campi elettrici e 
magnetici perpendicolari fra loro. La potenza assorbita dalla scarica è stata deter- 
minata in funzione della pressione e del campo magnetico. Le curve ottenute aumen- 
tando o diminuendo la pressione presentano un diverso andamento: la scarica è stata 
studiata in diverse disposizioni degli elettrodi. I diversi stati di scarica osservati sono 
stati qualitativamente interpretati riferendosi al comportamento medio di un elet- 
trone nella scarica. 
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An Analysis of Infrared Divergences. 


E. R. CAIANIELLO and S. OKUBO (*) 


Istituto di Fisica Teorica dell Universita - Napoli 
Istituto Nazionale di Fisica Nucleare - Sottosezione di Napoli 


(ricevuto il 15 Aprile 1960) 


Summary. — A concise and rigorous treatment is given ot the infrared 
divergences of quantum electrodynamics, by using a technique pre- 
viously developed in connection with ultraviolet renormalization. The 
result is essentially equivalent to an extension of the classical Bloch- 
Nordsieck formula, and is independent of perturbative methods. 


1. — Introduction. 


The present investigation is intended to provide a concise and rigorous 
proof, with methods recently developed by one ot us (1), of very plausible state- 
ments made on the problem of infrared divergences by several authors (2°) 
on the basis of the traditional analysis of graphs. It may serve, therefore, 
also as a test of the effectiveness of different techniques in questions of this 


nature. 
As is well known, unlike the ultraviolet divergences which are caused by 


short-range effects, the infrared divergences originate from the long-range 


(*) This work was done in part at the Ecole Normale Supérieure, Paris. The authors 
express their warm thanks to Prof. M. Livy for the kind hospitality they enjoyed there. 

« The research reported in this document has been sponsored in part by the Office 
Chief of Research and Development, U.S. Department of the Army, through its 
European Office under contract no. DA-91-591-EUC-1096 ». 

(1) E. R. CAIANIELLO: Nuovo Cimento, 13, 637 (1959); 14, 185 (1959). 

(2) J. M. JaucH and F. ROHRLICH: Helv. Phys. Acta, 27, 613 (1954) and Theory 
of Photons and Electrons (Cambridge, Mass. 1955), pp. 390-405. 

(3) D. R. YenwIe and H. Suura: Phys. Rev. 105 1575211057) 
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interactions among fields; one may say, in the usual phraseology, that the first 
are related to the innermost, the latter to the outermost part of the structure 
of the electron. It was shown recently (*) that ultraviolet divergences may be 
conveniently removed if the ordinary notion of integral is modified in a suitable 
way (provided the theory is «renormalizable »). We propose to show here 
that a similar approach permits to take care also of the infrared divergences, 
although the problem for them is not to remove, but rather to factor out of 
a matrix element the divergent contributions, so as to achieve an appropriate 
generalization of the classical results of BLocH and NORDSIECK (4); our work 
is thus similar in scope to that by YENNIE and SUURA (°). 

We restrict our considerations to the case of electron-electron scattering, 
with the explicit exclusion of forward scattering; this will suffice, since the 
generalization to all possible cases is immediate and does not involve, with 
our method, any complication of the calculations presented here. Reference 
is made, for notation and results, to the works quoted in Ref. (1), which were 
primarily concerned with the removal of ultraviolet divergences; since, 
however, the present problem is very much simpler than the ones treated 
there, we have tried to render the present work as independent as possible 
from those referred to, by limiting our quotations to a few preliminary for- 
mulae and repeating for the rest the considerations which are relevant for 
our present purposes. For the same reason, we consider only perturbative 
expansions and use equations among matrix elements only as short-hand for 
relations among their expanded forms; this is irrelevant for the final result 
and expedites our proofs. 

The reader will remark that there is in principle no interference between 
ultraviolet and infrared divergences, except for photon self-energies; only for 
the latter, indeed, can a continuous chain of short-range vacuum polarizations 
have an overall long-range effect in our problem. As we deal with a theory 
which is already freed from ultraviolet infinities by renormalization, any such 
interference is taken automatically into account by our treatment. 


2. — Infrared divergences. 


We consider the matrix element for non-forward scattering of two electrons 
with initial momenta q, and q, and final momenta pi and p,. 
It is given by (°) 


Pi Pa co 72 à Pi Do È or 
Gb) AA > o l= > Jas [as » RAD VE : [epcrcaailts 


We (Qn)! o o 
di VE ù | qh VE di UO Su 


(4) F. BLocH and A. NORDSIECK: Phys. Rev., 52, 54 (1937). 
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with: 


(2) (BE) =U); (EN = ul); (pa) =0, 


where u,(£) denotes the wave function of the electron of momentum 4. 

Ultraviolet renormalization amounts to changing the integrals i in (1) with 
a new operation, denoted with the symbol | in Ref. (*); this was proved in 
all generality in that work and will be the object of a specific study for electro- 
dynamics in a forthcoming paper. In order to avoid any possible confusion, 
however, we shall retain here the old integral symbol in (1) and in its conse- 
quences; the considerations that follow must be understood, though, as per- 
formed with the operations [ throughout, and with renormalized values of 
the parameters À (charge), ete. 

The 1-derivative equation for kernels (°) yields immediately, because of (2): 


d Pi Pa d d n Pi Ds È, fs 
(3) d(22/2) + a] PNA [€,, €2] NI È 


The integration DG gives infrared divergences whenever &,, é, are in 


positions such that: 


a) a pair Pi, Ps; OF qi, G2, is directly connected with &,, é; 

b) a pair p, q is directly connected with &, &; 

c) one momentum state p or q is connected to itself by & or È, through 
an ultraviolet self-energy part. 


Typical instances of these situations are shown in Fig. 1, to which all those 
those resulting from the appropriate symmetrizations must be added. ¢) shows 
the interference of ultraviolet 
with infrared divergences, which 
is however taken care of auto- 
matically by the use of | 
= our J) 

In analogy with the nomen- 


clature established in (1) we NES 

say that there is a confluence 4, 921 gi qo q,4 "A 
of the pair é,, é., with 1, Po: a) b) c) 
etc., whenever we find factors Fig. 1. 


in the integrand of (3) of the 

form (p,€1)- (P22) [462], etc. The preceding considerations have merely the 
purpose of illustrating in a familiar language the instances in which such 
confluences lead to infinite results; in all other cases one obtains finite results. 


(5) E. R. CAIANIELLO: Nuovo Cimento, 11, 220 (1954). 
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We neglect all other integrations which are contained in the perturbative 
expansion of (3) and decompose Il dé, Î dé, as follows: 


(4) fas fasie. G25 oe) = | as, | 46. (6, GODO 


fag, fac. f in (4) will be such that it coincides with Î dé, fas f whenever the 
latter is convergent (and then D“f=0); it will stay finite, however, if 
| dé, | fAË, f is divergent, the divergent part being given by D. Since we work 
SA a renormalized theory, there are no other infinities to worry about; it 
is interesting to remark that for non renormalized, and therefore for all non- 
renormalizable theories, this treatment of infrared divergences may break down, 
as will be seen from the fact that all divergent contributions in which 4 is the 
renormalized coupling constant are factored out by it as an exponential tactor 
exp [— (A?/2)1]. 

Tag | dé, is still far from being uniquely defined at this stage, nor do we 
wish to discuss here the optimum way to define it for practical purposes. We 
prefer to take advantage of the fact that our considerations apply as well 
to perturbative expansions, which permits to use in (4) a definition which 
agrees entirely with the standard procedure. The comparison of the results 
obtained by our method with those of the other authors will then be immediate. 

With | dé, LE we find in (3) always the integrand factor: 


1 
=z, exp Lik(E — &)] : 


(5) [661 =a, dre [a 1 k 


we define 


(6) ae, | a, viv [EE.lg(é1, È, ...) = lim | a, dé,- 


1,2 a 4 1 T(E 
WANA - (22)! Opus | fa k (UR exp [1k(é, — &)] + 


KR >x 


1 by 
i | dik =, (exp [ik(E,— &)] — 0° exp [io h(E: — Fan) Luce RCA 


— 
=] 
— 


DP yy [E, E, ]g(É,, Sr E = tim a. (ae, 
o—0 


eayl | i n 
DE la (22)! Ope 0? fase pesi exp [io k(É, SÒ) g(é.; + .….) ? 


|k|<K 
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where K is an arbitrary positive cut-off which is taken to be sufficiently small 
in the next section. 

Elementary computations prove that the def. (6), (7) satisfies all the wanted 
requirements (5): the new symbol i} differs from the old one | = if only when 
there are infrared divergences, which are set apart in (8). 


We find: 
E, d FES fas ae a ( au 
= see] =| oes EST TEL x + 
SERATA Perey 
Da Da lab 
2yiy2(Ee| Ml _ 
di VE Ei Si 


where the divergent term is computed explicitly from (1), (3) and (7), since 
the terms in the expansion of 


Le] [e] 
woe ey 
(o) [e] 
di LE È, DC ITS 
which give rise to infrared infinities such as 4), b), ¢) of Fig. 1, can be easily 


studied one by one. 
A straightforward calculation gives then (°): 


7 Pi Da bs bo Pi Do 
(9) inayat e SM |; 
di da Sa 3 di qo 
with 


(10) I=J(p1, Pa) +T(413 da) — IF (Pry da) — I (Po, la) — J (Pry de) — I (P25 La) + 


+4[F(p., pi) +I (po, Pa) +I (hh, LD) + IF (4, La) 
where 


Sit 1 Il 
a 


|RI<X 


(5) For example, in the case of the confluence of pair & and & with p, and p,, we have. 


DiPs E È, Ter, 
DE 717216185] M (' x La or 


il 
TOTP N e fat 
\1142 E È, = oo k 


|k|<x 


(ps — ok) — m (Pat ok) —m = c —ok, PB “) SER (Fs) 


(p, — ok)? + m?  (pa+ ok)? + m? eee do da do 


where jf has he same structure as M except for the absence of the Dirac matrices. 
“(p,) and (po). 


2567 


360 E. R. CAIANIELLO and s. OKUBO 


I is therefore given by 


(12) La r d*k agri ); 
with 

2 NA Pi Pa di do i 
4) ME Tio DO 


Exactly with the same arguments one would find, more generally: 


d Pa De Hy... Ur 
(14) 31/9510 N M Bes = 
HB) \e ps 


Pi Ds EE, Mie. Uy Pi Ds Dies Ur 
.M pr 


di qs Ses Wk ob OF di ds Ya Yr 


[ae dé, DEV Ree 


which is true also if y, =. 


8. — Factorization of divergent contributions. 


We can regard the r.h.s. of (3) as the result of an operation © performed 


D, De 
on M 1 a and write 
Mi Ye 
d Di Ds : Di Ds 
(15) d(22/2) M Si = O-M > o 5 
FU) NG Ge qd: 


© could be easily written down explicitly by means of functional derivatives, 
but this is unnecessary here. 
The decomposition (8) of 0, which we may symbolically express with 


(16) Vv tf À 6 V i" A Px c ") 
M — 07 an DM 
72 È o 0 >) a o 0 ? 
LA hh ds dn d 
leads immediately, if we define 
Pi Ds A. DD 
(17) M, cei oso = $ M o 
hh % dh q» 
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to 
d Di D APE PEUX] Di Pa 
(18) ae) À ai CD À 516 exp IE Ak 
de CCE CES Gi Qi do 
or 
a aay i | Dida bb 
(19) d(2/2) M, a = [as | tirarsi M, a , 
di do, si di UE E & 
if 
(20) 6D-M=D-6-M. 


The relation (20) is certainly verified by our prescription (6), (7); it suffices 
to use (14) to calculate first (d/a(22/2)) ©- M and then @- (djd(A2/2)) M: (20) 
follows by equating the two expressions thus obtained, due to the symmetry 
properties of the integrands. 

We conclude that the perturbative solution of (19), by iteration of 
fac. fag, is entirely free from infrared divergences, and that the original matrix 
element is related to the one thus calculated by: 


as a consequence of (9) and (17), with I given by (12). This proves that in- 
frared divergences, computed with the prescription used here (which is the 
close counter-part of the one usually followed), can indeed by factored out in 
an exponential factor. This result is furthermore independent of perturbative 
methods, because entirely similar considerations can be made directly on the 


branching equations (1). 


4. — Soft-photon emission. 


It follows from (21) that the cross-section for the process is: 
(22) o = exp [— /?-I]-o.. 
We must consider now the cross-section for the emission of » photons with 


94 - Il Nuovo Cimento. 
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momenta k,, ..., k, respectively. The corresponding matrix element is given by 


ky dos 1 = ze fa [a > ee ga y" Zalé) Biel, AIN 
Pi Pa DUO Con 
MUR A 
di VE gy proie Ean 
We are interested in the limit for vanishing k,...k,. Using the same termi- 
nology as before divergences with respect to & when k, >0 may appear in 


the cross-section when &, has confluences with either of the external p,, Pa, 
g, and q,. We define, similarly, 


(24) fas, = | dé, + DI, 


and neglect the finite contributions from ICT The remaining divergent part 
is easily calculated from (23) and (24), with considerations identical to those 
made previously, as 

ko A : 


where e, is the polarization vector of the photon k,: By induction we find 
that the divergent part of (24) is given by: 


(25) Z L a (pare) (81) ni 7 È Ds 


(nt (248)? (pk) (pet) (dk) (dh) 


ie} 


(26) An i L cal i (Pa 83) (41° 83) put 


min LE oui (Paks) (pk) (&k) (dk) 


je} 


di de 


so that the cross-section for » soft-photon emission has a divergent part 
given by 


(27) 6. = exp [— A2-T]-o., 
Ne 
where 
43 1 
28 = == 
(28) = RS foste, 
1k\<X 


with @ given by eq. (13). It is then easy to see that 
(29) HA 


because, since Q-k=0, Y(Q-e)? =Q?. 


pol. 
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Eq. (27) reduces thus to 


(30) On = ie exp [— f]:0., 


ne 


which is the desired Poisson distribution originally given by BLock and Norp- 
SIECK (4). The total cross-section is obviously: 


ao 
(31) DR N asia 
n=0 


which exhibits the cancellation of the infrared divergences against the di- 
vergent contributions due to soft-photon emission. 

Energy conservation is neglected in eq. (31). The condition Æ + 0 is un- 
realistic; if we impose, instead, for all n, 


(32) DL AE, 
i=1 
we have 
AH 
(33) Otot (AE) = F(B)-exp Blog GO 
where 
NA: 
Ze 2 
(34) 227) sap (4 RRO. 
and 
+0 1 
n di exp [it] _ of exp ata 
(35) F(B) 2 fas Ge XP È foe : = 
ies E 
Que SU 
= 1-5, f+ 0189) 


so that the divergent term 


(36) tes ake. 


disappears. 
The derivation of eq. (33) is almost the same as that given by YENNIE 
and SUURA (3) and JAUCH and ROHRLICH (7), and need not be repeated here. 
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In conclusion, we wish to emphasize that our aim has been here solely to 
give a rigorous and relatively simple derivation of eq. (21). We have limited 
our considerations to perturbative expansions and to the prescription (6), (7) 
in order to simplify the discussion as much as possible; neither restriction 
appears however to be necessary, and it may well be that other subtraction 
rules prove more advantageous. The prescription (6), as it is, could be taken, 
for computational purposes, as a formal re-definition of the free photon pro- 
pagator (which includes the additional limiting process 0 — 0): it is easy to see 
that, if everything is calculated with (6), performing the limit 9 +0 before 
the integrations gives M, performing it afterwards gives M.. 

One of us (S. O.) presents separately (“), for purposes of comparison, a more 
conventional treatment of this subject based upon use of (6) in momentum 
space, yelding identical results, although with less generality. 


RIAS SUNT O 


Le divergenze infrarosse della elettrodinamica quantistica vengono trattate in modo 
rigoroso e conciso, mediante tecniche elaborate a scopo rinormalizzativo. Si ottiene, in 
modo indipendente da metodi perturbativi, la naturale estensione della formula clas- 
sica di Bloch-Nordsieck. 
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(ricevuto il 16 Aprile 1960) 


Summary. — The energy spectrum of photoneutrons from bismuth irra- 
diated with 31 MeV bremsstrahlung is studied at 9 angles of 30°, 60°, 
90°, 120° and 150°, by means of nuclear emulsion. A comparison between 
the experimental energy spectra and an evaporative spectrum calculated 
according the expression F(E,) = KE, exp [— E,/T], shows a surprisingly 
high contribution of high energy neutrons around ÿ—90°. The angular 
distribution is isotropic in the region (2 < EH, <4) MeV, is approximatively 
of the form 1+0.7 sin? 9 for (4 < E,<5) MeV, and shows back and front 
asimmetry for E,>5 MeV. The contribution of the photoneutrons 
arising from the Bi(y, 2n) process is evaluated. 


1. — Introduction. 


The energy spectra of photoneutrons from Bi irradiated with bremsstrahlung 
have been studied by PRICE (?) at E.,ma = 22 MeV, by ZATSEPINA et al. (2) 
at Emax = 18.9 MeV, and by Bertozzi et al. (*) at.,ma = 14.3 and 15.8 MeV. 
These spectra show a greater contribution of high energy neutrons 


(1) G. A. Price: Phys. Rev., 98, 1279 (1954). 

(2) G. N. ZarsePINA, L. E. LAZAREVA and A. N. PospeLov: Sov. Phys. JETP, 
5, 21 (1957). 

(3) W. Berrozzi, F. R. PAoLINI and ©. P. SARGENT: Phys. Rev., 110, 790 (1958). 
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(E, > (3--4) MeV) which cannot be accounted for by an evaporative process 
than that found for instance in the photoneutron spectrum from Au (*). A 
remarkable contribution of high energy neutrons has been found also in the 
photoneutron spectrum from Ta (°), but this contribution is explained as being 
due to neutrons emitted by direct interaction and coming in part from a second 
resonance in the (y, n) cross-section. In the case of Bi, the known cross-sec- 
tions (5) give no evidence of a second resonance as in the case of Cr, As, I, 
Ta (59). So the high percentage of high energy neutrons from Bi seemes to 
be a peculiarity of this element. 

In the present experiment the energy spectra and the angular distribution 
of the photoneutrons from bismuth are studied by irradiation with a 30 MeV 
bremsstrahlung beam. At this irradiation energy, the study of the photo- 
neutrons includes a large region beyond the giant resonance in the o(y, n) 
process so that, far from the giant resonance, the contribution of other pro- 
cesses may also be investigated. 


2. — Experimental procedure. 


A cylindrical target of Bi was irradiated by a collimated y-ray beam from 
the Betatron of the University of Turin, working at the maximum energy of 
30 MeV. The photoneutrons emitted at angles of 30°, 60°, 90°, 120° and 150° 
with the photon beam where recorded by means of proton recoil tracks in L4 
plates 400 um thick. 

The experimental arrangement was analogous to that used in a previous 
work (1). The only difference being the substitution of the water wall—screening 
the plates against the neutrons from the betatron—with a paraffin moderator 
as shown in Fig. 1. 

The scanning of the plates exposed with no target in the beam shows that 
the background was negligible (<1%). 


HAE A SI E x 

(*) S. CAVALLARO, V. Emma, ©. MILONE and A. RuBBINO: Nuovo Cimento, 9, 
736 (1958). 

5 N Nr ATT i MAT 

(°) G. CORTINI, C. MILONE, A. RUBBINO and F. FERRERO: Nuovo Cimento, 9, 
85 (1958). 

6 SRN J 

J. HALPER? AT ‘Kee ! Ç IN 7 5 
(9) J. HALPERN, R. NATHANS and A. K. MANN: Phys. Rev., 88, 679 (1952). 


) 
(7) R. MONTALBETTI, L. Karz and J. GoOLDENBERC: Phys. Rev., 91, 659 (1953). 
) F. FERRERO, A. O. Hanson, R. Matvano and C. TRIBUNo: Nuovo Cimento 
4, 418 (1956). 
(9) F. FERRERO, R. Marvano and C. TRIBUNO: Nuovo Cimento, 6, 385 (1957). 
(!°) C. MILONE and A. RuBBINO: Nuovo Cimento, 13, 1035 (1959). 
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The experimental data are summarized in Table I. 


TABLE I. 


367 


Element 


209Bi (100%) 


(y, n) threshold (MeV) 


7.44 


Energy of the cross-section peak (MeV) 


Meg lev) 30 
Dose (Roentgen) 2 960 

d 30° 60° 90° 120° 150° 
Scanned volume (mm?) 50 54 88 O1 54 
Slow scanning. Number of experimen- Ee 5 ni 
tally accepted tracks E, <5 MeV PES DE 24) di 208 
Slow and fast scanning. Number of ex- ie DI ras 55 36 
perimentally accepted tracks E, > 5 MeV i 


\ 
\(Bi scatterer) 


0.5m 


Fig. 1. — Experimental arrangement; Pb: lead collimator; Paraffin (moderator) screen 
against neutrons from the betatron; L: photoplates; Bi: target of Bismuth. 


The experimental data on the observed recoil protons have been treated 
following the standard method (1) in order to derive the neutron energy and 


angular distributions. 
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3. — Results and discussion. 


31. Energy spectra. — The experimental spectra obtained at different @ 
angles are shown in Fig. 2. The curves reported in the same figure are calculated 
according the following expression, for evaporated neutrons (141°) 


F(E), = const E, exp [— ET], 


assuming for the nuclear temperature the value 7 = 0.8 corresponding to the 
T values found by many authors (1315). 


À = — aoe a = 
Bi(y,n) 
al 60° | 1200 | 1500 | 
£ 
= 
3 
asi | I ] 
g 
ë | | 
È 
N li | + 7 
4 ¥ 
0 SER iS ARTE : ue tite AE 
0 CO e YO ea ET ON 2 AN & € 


E (MeV) 


Fig. 2. — Energy spectra of the photoneutrons from Bi at angles #=30°, 60°, 90°, 
120° and 150°. The continuous curve is the calculated evaporative spectrum (see text). 


The evaporative curve (Fig. 2) is normalized to the average of the exper- 
rimental figures around H, = 2 MeV relative to all five experimental spectra. 

A remarkable contribution of high energy neutrons outside the evaporative 
curve is shown in Fig. 2. The amount of these neutrons divided by the total 
neutron number over H, = 2 MeV is reported in Table II as «neutron frac- 
Clone) (ES) 


(11) V. F. WrissHoPr and D. H. Ewing: Phys. Rev., 57, 472 (1940). 

(12) D. L. Livesey: Can. Journ. Phys., 33, 391 (1955). 

(8) G. CoRTINI, C. Mitonr, T. Papa and R. RinziviLLo: Nuovo Cimento, 14, 
54 (1959). 


(4) PoC. Gueruorey Phys. mev., 81, 51 (1951), 

(15) L. Corti, U. Faccuini, I. Iori, M. G. Marcazzan and A. M. Sona: Nuovo 
Cimento, 18, 730 (1959). 

(16) V. Emma, C. MILONE and R. RINZIVILLO: Nuovo Cimento, 14, 1149 (1959). 
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TABLE II. — « Neutron fraction F » at E, > 2 MeV for Bi irradiated with 30 MeV and 
18.9 MeV bremsstrahlung. 


Ÿ angle | 30° and 150° | 60° and 120° | 90° | References 


0.32 © 


=| | 


«Fy (30 MeV) | 0.22 0.38 Present work 


0.58 (2) 


«Fy (18.9 MeV) 0.44 


In Table II is also pelle ee 
reported the neutron frac- 
tion F relative to the + Ey =30 MeV (Pr. work) 
experiment at En, = Ms 
— 18.9 MeV (). The eva- 
porative curve and the 
point of normalization are 
the same as in the exper- 
iment at 30 MeV. This 
fraction F, is surprisingly 
high in a nucleus as heavy 
as bismuth. 

Yet, the contribution 
of high energy neutrons in 
our experimental  spec- 
trum is lower than that 
found by ZATSEPINA et 
al. (2) at ÿ = 90° and 270° 
with 18.9 MeV bremsstrah- 
lung. This may be seen in 
Fig. 3 where the two exper- 
imental spectra obtained tg ll A EEE ; 
at E,...=20 and 18.9MeV Al ee TI NIE tr 
are compared after nor- ù Sa LO 15 
malization in the ener- Fig, 3. - Comparison between the experimental spectra. 
gy range 5<H,<8 MeV obtained at E,max=18.9 (?), 22 (°) and 30 MeV after 


(see below). In Fig. 3 normalization in the energy range (5 <H,< 8) MeV 
(see Sect. 3°3). 


Mu 


"(Ref 


Arbitrary units 


dio) 


is also shown the energy 
spectrum of photoneu- 
trons from Bi obtained by PRICE (1) with irradiation at E.,na = 22 MeV. 

We will show in Section 3°3 that this different behaviour is due to the 
Bi(y, 2n) process that for Ema: = 30 MeV gives a considerable number of 
neutrons with energy E,< 5 MeV. 
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We consider the (y, 2n) process as the only process competing with the 
(y, n) in the neutron emission. Indeed at the excitation energy that we have 
used, the (y, np) process is negligible respect to the (y, 2n) because the emis- 
sion of a charged particles is strongly inhibited by the Coulomb barrier, and 
the Bi(y, 3n) process is not appreciable at this energy. 


32. Angular distribution. — The angular distributions are reported separa- 
tely in Fig. 4 for the following neutron groups: 2<H,< 4 MeV, 4<E,< 5 MeV 
and H,>5 MeV. The experimental results are normalized to unity for ÿ = 90°. 

Neutron group 2<H,<4 MeV: 


sl 1 i dre The angular distribution is isotropic 
| ari È ai | inside the experimental errors, as ex- 
OBF > i | pected, the emission of the neutrons 
| . su k 3 i being due principally to evaporation 
0.4} i + | in this interval of energy. 


+ Present work 
+ Reference (8) 


8 roc cle | FI 
& 087 Ae os 1 
% piste IL De 
hg rt 7 Sassi È 
04 | 
(4<E 5) MeV il 
0 - —— 
Tori =) ] 
mort set sar Parce +---- ES 
0.8} | 5 
| = 
- 1 a 
| =e IT 7 
04° I 
| (2<E,< 4) MeV | oy, 2n) 
| 3 
o a SERI RR, e | 
30° 60° 90° 120° 150° 10 5 20 Mey 
È 
ay 
Fig. 4. — Experimental angular Fig. 5. — Bi(y, n) and Bi(y, 2n) cross sec- 
distribution of photoneutrons from tions deduced by the total photoneutron 
bismuth for the various neutron cross sections assuming the statistical 
groups. theory of nuclear reactions. 


Neutron group 4<E,<5 MeV: These photoneutrons may be attri- 
buted prevalently to y absorption in the region of the giant resonance; see 
the behaviour of the o(y,n) in Fig. 5. 

An angular distribution of the form A+B sin? 9 should be expected for 
electric dipole transition, The experimental figures may be fitted by a curve 
of the form A+B sin? 4 with (B/A) ~ 0.7. 
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Neutron group #H,>5 MeV: A back and front asymmetry appears 
evident in agreement with a similar result found previously (8) with Si(n, p) 
detectors (Fig. 4). 


3°3. Contribution of the Bi(y, 2n) process.,— The competition between the 
(y, n) and (y, 2n) process can be predicted by the statistical theory of nuclear 
reactions (1417) and we use as a good approximation the following expression (17) 


o(Y, 20) 
oY, Da) 


(1) 


exp 


= 1 L | Ee: E,) 


e i 

Al 5 BB], 

where a is a constant for a given A, E is the energy of excitation, E, is the 
threshold energy for the (y, 2n) process and o(y, n,) =o(y, n)+o{y, 2n) ac- 
cording to our previous discussion. 

This expression compares the (y, 2n) contribution to the (y, n) cross-section 
in heavy nuclei; it may be observed that the calculated ratios o(y, 2n)/o(y, n,) 
are in fair agreement with the experimental results (71819), 

The experimental total neutron cross-section curves, OygxperltYs 0,), for the 
photoneutron emission from bismuth are given by several authors: up to 
20 MeV (7:8) and up to 25 MeV (*); so the unknown o(y, 2n) cross-section may 
be calculated according to expression (1). Let 


o(y,2n) _ o(y, 20) 
Glyn o Nha oly, 2n) 


= Î(E y); 


and being 


Opxper(Y; n) =. oly, Hae 20(Y, 2n) 
o(y, 20) o(y, 20) À 


We then have 


0(Y, 20) = Orxper(Y, De) ga i 
1-f(E,) 
ott) cn nd (14 155) 


We have introduced for this calculation the values of the o,,,..(y, n,) found 
by FERRERO et al. (3) up to E. = 20 MeV, that are in very good agreement 


(7) B. T. FeLp, H. Fersupacu, M. L. GOLDBERGER, H. GOLDSTEIN and V. HE 
Wersskorr: Final report on the fast neutron project, USAEC NYO 636 (1958). 

(18) J. GOoLpENBERG and L. Karz: Phys. Rev., 90, 308 (1953). 

(19) E. SiLva, J. GOLDENBERG, P. B. SMITH and L. Marquez: Nuovo Cimento, 
9, 17 (1958). 
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with that obtained by MONTALBETTI et al. (7). For E, between 20 and 25 MeV 
we have adopted the values given by HALPERN et al. (°). 
The ratio 


EYmax EYmax 


fon n) ab] [ot 2n) dE, 


so calculated is reported in Table III. In the same Table, the experimentally 
known ratios for 181Ta (2°), 141Pr (21), ND (19) and Cu (21) are also reported. 
As Table III shows, the calculated ratio for Bi is very near to the experimental 
ratio for Ta as theoretically predicted. This result gives assurance of the relia- 
bility of the calculated Bi(y, 2n) cross-section taking into account that the 
behaviour of the o(y, n) cross-sections are very similar for Bi and for Ta (?°). 


TaBLe III. 
Element IS Gu ES DEN Ea 181Pa 200 By 
#,(y, n) (MeV) | 10.6 8.7 10.0 7.6 7.44 
E, foe sa max (MeV) ve (An 2 15.2 14.0 133 
| E,(y, 2n) : (MeV) | 2 OTO ies 18.3 14.0 13.3 
E. for 0.,,2n Pa (MeV) | 25 D || es} 20 15 15.5 
E,(y, 2n) —E,(y,n) (MeV) | 9.4 | 8.6 8.3 6.4 5.9 
| VIE (MeV) 35 22 31 31 | 25 25 
RE nee alal ino Experimental Calculated 
foro n) ar) [at 2n) dE, = - 
li : 6.80| 6.25| 5.5 | 3 | 2.7.0.3 |'2.6540.3 2.9 
= o me) ; (9) a) (22) (29) Pres. work 
= |- 


The deduced curves of the B i(y, n) and Bi(y, 2n) cross-sections are shown 
in Fig. 5 together with the o(y,n,) obtained by FERRERO ef al. (5) up to 
E,=20 MeV and by HALPERN et al. (*) for E, > 20 MeV. 


. H. Carver and W. TURCHINETZ: Proc. Phys. Soc., 71, 613 (1958). 
. H. CARVER and W. TuoRCHINETZ: Proc. Phys. Soc., 73, 110 (1959). 
J. BERMAN and K. L. Brown: Phys. Rev., 96, 83 (1954). 


A ranean: R. MaLvano, E. SiLvA, J. GOLDEMBERG and G. Moscati: Nucl. 
Phys. 10, 423 (1959). 
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Now in order to show and confirm that the different behaviour of the 

spectra reported in Fig. 3 is due to the contribution of photoneutrons arising 

from the (y, 2n) process, our spectra and Soviet group’s spectra observed at 
9 — 90° are connected by the following functions: 


‘ 


dg, = OY, n)-I(H,, 30) ; tig = Oly, n)-I(E,, 19) ; 
bs, = 2o(y, 2n)-T(H,, 30); Br» = 20(y, 2n)-I(E,,19). 


Normalizing these func- 
tions so that a, and dy, 
assume the same max- L 
imum value we observe 
that the curve a, nearly 
coincides with the norma- 
lized curve a® (Fig. 6). 
Taking into account this 
result, we have normal- 
ized the two experimental 
spectra in the energy 
‘ange 5< HE, <8 MeV; the 10 15 Ey 20 MeV 
lower limit (E, >5 MeV) pos Du te In A ai ii 

0 5 10NÉ ES CV DIN IS 
selects the neutrons aris- 
ing almost exclusively Fig. 6. — Functions 44, and 4%, 28 defined in 
from the (y, n) reaction, Sect. 3°3. 
while the higher lmit 
(E,<8 MeV) excludes the low contribution of neutrons coming certainly 
from the region in which aj, becomes different from a, (Fig. 6). 

Now, if F(E,) is the flux of neutrons of energy É,, we have 


T 


units 


Arbitrary 


30-E;(Y; 2) 30 
F,(E)dE, [tes + bs) d E, 
0 Es. 0) 
19 — E;(Y, n) et 19 
(a) | pt) 
(En) OE, fa 19) AE, 
0 Es(Y: D) 


Ignoring the negligibles values of the integrals 


E;(Y, 2) +1.5 30 30 


adE,, fond, and far. (Fig. 5); 


E;(Y: 0) 25 25 
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taking no account of the 
Ey(Y: 2n) +3 


bdk,, 
E(y, 2n) 
values that give prevalently neutrons of energy less than 1.5 MeV; and ob- 
serving that the F(E,) values for E, <1.5 MeV are connected to the F(E,) 
values in the energy range (1.5< E, <5 MeV), we have with good approx- 


imation: 
30—Z.(Y; n) 25 25 

| P;,(H,) dk, i 3) dH, + | bs) dE, 
R 1.5 | _ __ Esty.m)+1.5 E(Y: 2n)+3 er 
1 ~~ 19-E;(Y; n) 19 19 2 

Ù (n) ; 
F9 (Ln) dEn i a dE, + | bo dE, 
1.5 Es(Y. n) +1.5 Es(Y. 20) +3 


In our case the ratios R, and R, obtained from the curves reported in Fig. 3 
and Fig. 6 assume respectively the values 1.4 and 1.3. This confirms that the 
different behaviour between our and the Soviet group’s spectra, in the energy 
range E, <5 MeV, is due to the different contributions of the (y, 2n) process. 


4. — Conclusions. 


a) Photoneutron spectrum from Bismuth shows around 90° a surpri- 
singly high contribution of neutrons outside the calculated evaporative spectrum. 


b) The angular distribution is isotropic for Æ,<4 MeV, is approx- 
imately of the form 10.7 sin? 9 for (4< E, <5) MeV while a back and front 
asymmetry appears for E, > 5 MeV. 


c) Evidence and evaluation are given of the photoneutrons from (+, 2n) 
processes. 


RIASSUNTO 


Si studia lo spettro a 30°, 60°, $0°, 120° e 150° dei fotoneutroni emessi dal bismuto 
irradiato con bremsstrahlung di 30 MeV. Un confronto degli spettri sperimentali con 
uno spettro evaporativo calcolato secondo l’espressione F(E,) = cost-E,-exp [— E,/T] 
mette in evidenza intorno a 90° un contributo di fotoneutroni, al di sopra della curva 
evaporativa, sorprendentemente alto per un nucleo cosi pesante. La distribuzione 
angolare risulta praticamente isotropa per neutroni con energia inferiore a 4 MeV. 
Nell'intervallo (4 < E,<5) MeV si ha F(3)=14-0.7 sin? è mentre sopra 5 MeV si nota 
una evidente asimmetria avanti indietro. Si 


DRS valuta il contributo dei neutroni prove- 
nienti della reazione Bi(y, 2n). 
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Quantum Electrodynamics in an Expanding Universe (*). 


F. L. SCARF 


University of Washinglon - Seattle, Wash. 


(ricevuto il 19 Aprile 1960) 


Summary. — The equations of quantum electrodynamics are investigated 
in an expanding universe characterized by a de Sitter metric. It is 
assumed that the expansion occurs on a microscopie level and that 4, 
the bare electron mass and the bare charge are all constant. For a par- 
ticular coordinate system the operators may be expressed in terms of 
coupled flat-space fields which have m replaced by mR(t). The equivalent 
quantized field problem is soluble using standard perturbation theory 
with modified commutators and Green’s functions. As a first approx- 
mation, adiabatic changes caused by the expansion are evaluated by 
imaking the replacement m +mR in all final formulae. Renormalization 
effects are examined using asymptotic expressions which contain ultra- 
violet cut-offs. Assuming that the cutoff procedures are meaningful and 
are related to the structure of the physical electron, it is noted that the 
renormalized mass and charge could vary with time for a particular 
kind of cut-off process. 


1. — Introduction. 


It is customary to interpret the observed red shift of stellar spectra as a 
Doppler effect associated with a large-scale expansion of the universe. This 
expansion has been described as a Newtonian phenomenon, and it also has 
been discussed using a curved space-time metric such as 


(1) ds? = di? — R(t) (da? + dy? + de?], RO) 210. 


(*) Supported in part by the U.S. Atomic Energy Commission. 

(1) The velocity of light is defined by this metric; it is then equal to unity and 
independent of cosmic time in the conformally flat T-system introduced below. Fur- 
thermore, we assume that # is the same universal constant of unit magnitude in 
all systems. 


2583 


376 FORTS SCARF 


For macroscopic motions the two points of view predict the same effects (2). 
This note is concerned with a comparison of the two theories on a microscopie 
level; we shall examine the consequences of an assumption that the relativistic 
expansion occurs for atomic and subatomic phenomena. 

In order to have any observable red shift for electromagnetic radiation, some 
unit of length (that is, a third scale factor to go with c, #) must be independent 
of the cosmic time {. In the past it has usually been implicitly assumed that 
the bare electron mass is the meter stick which defines (dx) of eq. (1) and we 
shall accordingly take m, and e (bare charge) as parameters which are constant 
in all systems. 

Ordinary quantum mechanics then predicts that only the standard red 
shift will be seen, even if the Bohr orbits expand; when the co-ordinates in 
the Schrédinger equation, | 


(2) Ey = [p°]2m, — Zei/rly, 


undergo an adiabatic variation which is consistent with eq. (1), r is repiaced 
by E(t)r(0) [(dr?) > R?(t)(dx? + dy? + dz?)| and p, by x(0)/R(t) (the latter fol- 
lows from invariance of the phase integrals). This leads to 


(3) E,(t, m) = E,(0, m) R(t))/R(t) + O(dR/dt) , 


but for a Coulomb force, 


(4) E,(0, m') = m'f(é, h, n), 


so that E,(t, m,) is the same constant at each star, if dR/dt terms are negli- 
gible. The only appreciable modifications of the observed spectrum stem from 
the relative velocities of source and observer. The same result follows when 
the adiabatic variation is applied in the classical relativistic hamiltonian. In 
fact, one finds that Z, is independent of R(t) to first order using only the virial 
theorem for the Coulomb potential. 

Although these simple arguments might seem reliable when applied to 
determinations of atomic energy levels, there are several uncertainties con- 
nected with the above discussion. The results, eqs. (3), (4), depend in a de- 
tailed way on the form of the potential. It would be desirable to have a for- 
malism which can treat any problem in electrodynamics, including those with 
magnetic forces, spin-orbit forces, etc. Furthermore, the Schrédinger and 
Dirac equations must be regarded as approximate relations which follow 


(*) H. BonpI: Cosmology (Cambridge, 1952), chap. 9; W. H. McCrea: Rep. Progr. 
Phys., 16, 331 (1953). 
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from the complete quantized theory of interacting electron and photon fields. 
Accordingly, in the next section the general equations of quantum electro- 
dynamics are examined for a geometry characterized by eq. (1). It is found 
that the original problem may be mapped onto one Having effective flat-space 
fields with mg = mo R(t). 

In the third section the analysis is completed with a discussion of renor- 
malization effects. It is noted that another fundamental length (the cut-off, A) 
enters at this point and the renormalized charges and masses are functions of 
(A/m)). For some cut-off models, À has a value which depends on Rit). We 
conclude that the quantized field theory may lead to results which are not 
contained in the single particle picture, and that frequency variations other 
than Doppler shifts could be observable even if the parameters for bare fields 
are constants. 


2. — Field quantization. 


In a Riemannian space with metric g,,, the equations of electrodynamics are 


(5) iy" Vp + mp = 0A, , 

(6) {y", y} = 29", 

(7) (—g)*, [(— 9) F4] = ew", 
(8) F,,=0,4,—04,, 


(g = detg,,). The covariant spinor derivative is defined so that V, acting on 
yy or yy'y gives the ordinary gradient of a scalar or vector, respectively. 
One finds (*) V,y = 0,y— Ty, 


(9) Ir y, VO 


where 71%, is an ordinary Christoffel symbol. 


For the present it is convenient to consider these equations in a co-ordinate 
system (T-frame) with 


(10) ds? = h(T)[AT: — de — dy? — de®], 


TL 
t—[h(T) AT, R(T) = R(t). Then g,, = h(LG,, (I, is flat) and eq. (6) may be 


(3) D. R. Britt and J. A. WHEELER: Rev. Mod. Phys., 29, 465 (1957). 
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satisfied by y" = AT) (9 is flat). In this T-frame the spinor derivative 
also has a simple form, y"V,y = hapa (hey), and eq. (65) becomes 


(11) "0,9 + h(Lymop = PA, 9 ; 
with g=h*y. In the same system eq. (7) reduces to 
QTD) 2 0, (hE) = elP7"9) - 


The tensor F” must be expressed in terms of /’,,, the quantity which de- 
fines A,. The transition is quite straightforward: 0,— In (Gu. is Con- 
stant) so that the left-hand side of eq. (12) is 0*(ÿ,,hk2F#). If the entire equa- 
tion is multiplied through by ÿ,,, the new source term, ¢(~Y,;@), is equal to 
002,0, FR) However =, 4009-20 h, so that 


(13) 00,4, = e( PY, ) 


if we require that c”A, be zero in this particular co-ordinate system (*). 
The last condition is not covariant and therefore it is not the Lorentz 
condition (4,” = 0). However, although eqs. (11), (13) are equations of mo- 
tion for field operators, the gauge condition is only a subsidiary condition on 
the state vectors, (0’A,)‘*?|m> = 0. Moreover, the definitions of all state vee- 
tors, positive and negative frequencies, etc., shall be made with respect to 
a particular co-ordinate system, and our subsidiary condition does not intro- 
duce a new lack of general covariance. (We shall not speculate here on the 
fundamental problem of reconciling the non-covariance of quantum field theory 
with the general covariance of the relativistic scheme; we merely note that 
this difficulty is common to all theories [see ref. ()].) In the T-frame it insures 
that the things which are ordinarily called longitudinal photons are eliminated, 
but the «energy » of the free photon field is not positive definite (5). Since. 
the concept of energy is degraded in a curved space we believe that the effect 
is not necessarily significant. At any rate the quantization of this field theory 
is straightforward and consistent although it may not correspond to the gene- 


(*) SCHRÔDINGER has mapped the source-free electromagnetic field strengths into 
flat space quantities for conformally flat spaces (E. SCHRÔDINGER: Proc. Irish Acad., 
46, 25 (1940); see also M. von LAUE: Helv. Phys. Acta, Suppl., 4, 42 (1956)). 

(9) GUTZWILLER has examined these equations in terms of a potential 4 defined 
by Pur— Qu A”— Ar (M. GurzwILLER: Helv. Phys. Acta, 29, 313 (1956)). His 4 is 
not related to our A by any simple equation such as 44— Ju A, Furthermore, in this 


work the condition A”,,=0 was used so that the final equation for 4 involves 


I 
. . . (ae 
in an explicit and .complex manner. 
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ralization of ordinary electrodynamies. It should be noted, however, that in 
the relevant adiabatic limit of h—0, the two conditions are actually 
identical. 

The coupled equations for p and A have flat-space derivatives and matrices 
and are identical to the flat-space equations for a g-electron field with the 
single modification m, = mh* (TL). In the T-frame it is possible to quantize 
eqs. (11) and (13) by imposing the commutation rules 


(14) LA), At (e)| = ig.,, D(æ— x") = Gy Doe ve 
and 
(15) {pa (2), Gs (@')} = STE 


~ 


Here D, which satisfies DD — 0, is the ordinary flat-space photon commutator. 
The anticommutator S is not the conventional S function, but is a solution 
of [1ÿ,0,+m(T)]S = 0 with the usual boundary values. (If we had required 
A,” =0, the D function in eq. (14) would be a solution of L(0,, h)D = 0.) 

The quantization in the T-frame is completed by setting up Schrodinger 
equations for the operator representatives of the classical constants of motion. 
The action principle which gives the equations of motion leads to 


+ Ga Na AV 
(16) OT = +37 (09, / 0%"), 
where 7% = v—gT! has the usual form with covariant derivatives replacing 


the flat-space quantities. Since h = A(T), the spatial component S of the total 
momentum, 


(ere) | P,=[aeey= gr" («=1, 2, 3), 


are constants of motion (in both the t and 7 system) although P, and P” are 
not. The state vectors in the T-frame are thus the simultaneous eigenstates 
of P, and 


(18) Q= | da j, (ej, = 0), 


and they may be constructed in the customary manner. For instance, the re- 
normalized in-state with a single physical electron is (°) 


(19) ive lim | D, (x, T)|OHHo(x, T), 


(6) W. THIRRING: Principles of Quantum Hlectrodynamics (New York, 1958), p. 177. 
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where v is a positive frequency ¢-number solution of [1p"0,,+ m(T)]v, = 0. 
. È 2 A (+) (et) PA) a 
The vacuum is defined so that the above is zero when v,” > v,’, and v,- are 


given by 


9 Aber hy eta AN 
(20) Ve (4, %) = Dari (Tp E Ge) c 


If the gauge condition 0 À, = 0 is used, the photon states are defined in the 
conventional manner. In the above scheme all state vectors are constructed 
in terms of entities which are not spinors or vectors under general co-ordinate 
transformations. This means that a special kind of time keeping is required 
for the decomposition into positive and negative frequencies (7). 

The general perturbation treatment of quantum electrodynamics is again 
valid for this theory with modified S function and suitably modified propa- 
gators. Since energy is not conserved, the out-state has additional pairs and 
photons and the matter production is proportional to (dh/d7) (eq. (17) with 
y—0). We are primarily interested in the adiabatic limit because È is extre- 
mely small for the actual universe. Therefore, it is of interest to compare the 
predictions of the theory for two values of h, neglecting terms of order h. 
(Again note that in this adiabatic limit 0’A,—= 0, does not differ from the con- 
ventional Lorentz condition.) 

If è is very small, P, depends on T only through the adiabatic variation 
of m,, and therefore P,(T, ms) = Po(0, mh*)+0(h). However, P° =P, varies 
as h-!; in order to have a meaningful definition of energy we must transform 
to a system such as the ¢frame with g°=1, g°*=0. In the t-frame the 
energy is P, = P= P,/R(t) and thus the adiabatic limiting process applied 
to the field theoretical problem again gives eq. (3), but it is now meaningful 
for any electromagnetic phenomenon. The result is more general than was 
indicated by the original deduction and it is valid for hyperfine splittings with 
anomalous g-factors, radiative corrections and so forth. 


3. — Renormalization. 


All properties of the physical particles (y-field) are given in terms of the 
observed (renormalized) charge and mass with e — Z'e, and m=m,+5m. We 
assume that m, and è, are finite and find, as usual, that 3m and Z are repre- 
sented by divergent integrals. Several authors have suggested that there may 


(*) For many forms of R(t) the out-state has unusual properties in at least one 
of the two frames. A discussion of this point shall appear in the Proceedings of the 
Royaumont Conference on Gravitation (C.N.R.S., Paris, 1960) 
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be some degree of meaning associated with the introduction of an ultra-violet 
cut-off although there is no unique way to perform the operation. We shall 
assume that observations (|¢|> 0, m< oo) imply that the true theory does 
contain an additional length in some unknown manner; we also assume (with 
much less foundation) that a cut-off theory has significance as a phenomeno- 
logical way of accounting for this new parameter with dimension of length. 

These ideas may be used to speculate on the possible effects of a micro- 
scopic relativistic expansion on Z, dm. LANDAU (8) has argued that formulae 
of a certain type are meaningful so that, for instance, in the T-frame (è small) 
the renormalization constant for the (effective) g-field is 


(21) Li (=) = exp [— 2e In (A,/m,)], 


where m, is the slowly varying effective mass for the y” field and A, is an 
effective cut-off energy for the virtual quanta of the bare y and À fields. 

Some physical meaning must be attached to di in order to complete the 
analysis for an expanding universe. A plausible suggestion is that A, is à 
constant multiple of the effective bare electron rest mass and that it corres- 
ponds to the maximum energy of the effective virtual electrons. Then A, km 
is constant and A, = km, = h}(T)km,, so that A,/m, =A,/m, =k and Z is 
independent of 7. If this were true, there would be no real difference between 
the predictions of the single particle and second quantized theories. (One 
argument in favor of this choice is based on the observation that all charged 
particles have the same value for e. This suggests that Z should indeed be a 
universal constant. However bosons have the same e, but eq. (21) does not 
describe their renormalization.) 

On the other hand, A, might represent a maximum energy independent 
of T (that is, À, might change). For instance, A, could be the energy asso- 
ciated with a critical photon frequency or field strength, hr (the photon 
properties are independent of T for h—0 but the y—A coupling term 
contains A(T); on the other hand, the g — A coupling is independent of h). 
It is not unreasonable to associate the cut-off for the electron charge with some 
property of the photon field since the divergence is certainly an aspect of the 
interaction between the two fields. In this case Z(T)= Z'(t) is a scalar func- 
tion of x, which depends only on time in the two special frames under con- 
sideration. 

While the second suggestion may be highly improbable, it is of some in- 


(8) L. D. Lanpau: Niels Bohs and the Development of Physics (New York, 1955), 
pp. 52-69. A somewhat different form for Z is given in this reference. 
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terest to note that the variation 


Di 4 
(22) Z(t) = éxpi= EW [Z (present) |-1 In |4/R(6)] 


will give an increase in e with time (the present value of (137)-1 has been 
inserted for x and & is defined as fia,,,/m)). This would mean a decrease with 
retardation (distance) for the atomic frequencies in stellar spectra. Such a 
decrease could be separated from the ordinary Doppler shift by examining 
the ratios of fine structure frequencies to principal frequencies for sources of 
discrete spectra such as planetary nebulae. If one inserts E(t)= exp [é/r], 
eq. (22) becomes 


(23) Z(t) = Z(0) exp [2t/137 rZ(z)] 


and the maximum frequency variation over the observable universe corres- 
ponds to bitin 0, tn t (*). The self energy also varies with time if 24 
is constant, but these changes cannot be isolated from the standard Doppler 


shift. 


4. — Conclusions. 


If all lengths and times, except for the bare electron Compton wavelength 
(A,), expand uniformly (e.g., Ar=h?Ay,, —AT=htAT,) one might expect 
that electrodynamics could be connected with a theory in which only the 
effective electron Compton wave (A,) decreases while other lengths remain 
fixed (Ar = Ar,, AT’ =AT,). We have shown that this reasoning is indeed 
correct if a special non-covariant generalization of the Lorentz condition is 
imposed, but that the connection does not exist in general (unless h is neg- 
ligible). 

It has also been noted that the renormalization « constants » of quantum 
electrodynamics might be functions of time, even in the adiabatic limit. In 


(°) A frivolous estimate of the greatest permissible variation in e? can be made 
using these formulae if one assumes that kR-1(t) had its maximum real value (no ghosts), 
at t-0. Then Z(0)=exp[— 2 (In k/137) Z(r)] and Z(t) =exp [2(1 — In k)/137Z(x)]. Set- 
ting Z(t) =COZ(0) we get Z(0) In Z(0) = —21nk/1370. The smallest real value of Z 1n Z 
(largest real k) is —e-1 = — 0.368. The equation for Z(t) then yields 0 In O=2¢/137 
or C=1.34. This magnitude for O gives a maximum Ay comparable with ordinary 
Doppler shifts, but any numerical result is, of course, extremely sensitive to the func- 
tional form of the cut-off factor. 
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order for this to cecur the ultraviolet cut-off must represent some maximum 
energy for the virtual photons; if the cut-off is correlated with a maximum 
bare electron energy the renormalization constants are time-independent in 
this limit. 


I am grateful to Professors E. HENLEY and L. WiLets for several helpful 
discussions. 


RIASSUNTO (*) 


Si esaminano le equazioni dell’elettrodinamica quantistica in un universo in espan- 
sione caratterizzato da una metrica di de Sitter. Si suppone che l’espansione avvenga 
a livello microscopico e che %, la massa nuda dell’elettrone e la carica nuda siano tutte 
costanti. Per un particolare sistema di coordinate gli operatori possono essere espressi 
in termini di campi accoppiati dello spazio piano che hanno m sostituito da mA(t). 
Il problema equivalente del campo quantizzato si risolve usando la teoria standard 
della perturbazione con commutatori e funzioni di Green modificati. In prima 
approssimazione i cambiamenti adiabatici prodotti dall’espansione vengono valutati 
facendo la sostituzione m >mR in tutte le formule finali. Gli effetti di rinormalizza- 
zione vengono esaminati usando espressioni asintotiche che contengono tagli ultra- 
violetti. Supponendo che i procedimenti di taglio siano significativi e siano in rap- 
porto alla struttura dell’elettrone fisico, si nota che la massa e la carica rinormaliz- 
zate possono variare col tempo per un tipo particolare di procedimenti di taglio. 


(*) Traduzione a cura della Redazione. 
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Tripole Ghosts in Field Theory. 


Keele INVAGERT 


Istituto di Fisica e Scuola di Perfezionamento in Fisica Nucleare dell Universita - Roma (*) 


(ricevuto il 22 Aprile 1960) 


Summary. — An appropriate modification of the Lee model is discussed. 
For a definite set of parameters beside the dipole ghost state also a 
tripole ghost is to be found. Using the Källén-Lehmann method contrib- 
utions of tripole ghost states, supposed to exist in a relativistie quantum 
field theory, are calculated. 


1. — Introduction. 


According to the general properties of a state vector space with an indefi- 
nite metric even in a finite dimensional space, it is not sure that the eigen- 
vectors of a hermitian operator (in Heisenberg’s nomenclature pseudo-hermi- 
tian) form a complete system. The necessary and sufficient condition of the 
completeness is that the minimum polynomial m(P) (i.e. the polynom having 
the least possible degree for which m(P) = 0) of the operator P to be con- 
sidered should have the form (1) 


MP) = (P—p,)... (P— px), DPix~P, Lit 


Since the characteristic polynomial M(P) of P is, of course, equal to zero 
and has the form 


EM (PY =" (Pay et p2E Da ee 


(*) On leave of absence from Institute for Theoretical Physics of the R. Eòtvòs 
University, Budapest. 


©) L. K. Panpir: Suppl. Nuovo Cimento, 11, 157 (1959). 
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where «,+...+ x, is the dimension of the vector space, the minimum poly- 
nomial can always be written as 


m(P) = (P — p,)*... (P — pt, Bi <a. 


But if at least one of the Ps, let us say B;>1, the eigenvector |p,) cor- 
responding to the eigenvalue p; has necessarily zero norm. If /; = 2, a state 
vector different from 0, satisfying 


(1) (Pepi Dp, 


also exist and has to be added to the eigenvectors, for completing the vector 
space. The state |D> was called by Heisenberg dipole ghost. If 6; =3 be- 
side |D> another state vector |T) defined by 


(2) (PERIOD 


has also to be added to the set, and so on, depending on the degree of the 
minimum polynomial. |7> may be called « tripole » ghost. 

Choosing the parameters of the model appropriately, a dipole ghost state 
can be found in the Lee model (2). Ghost states with higher pole order 
emerge (*) in a proposed model theory of FROISSART (*), but in the sector 
containing higher number of particles. 

Here a modification of the Lee model will be treated resulting in dipole 
and tripole ghost states. The modification was originally made in TEODORO- 
vicu’s work (5), but there the consequences were analysed in a completely 
wrong way. Next we calculate the contribution to the propagators of tripole 
ghost states supposed to exist in a relativistic field theory. 


2. — The model. 


The modified Lee model starts with the Hamiltonian 


3) f=— ms | veyed + n [va + Jo(k) a*(k) a(k) dk — 


Jo 1 ee. : 
Tei k) {wx (q)vv(p) a*(k) — y(p) px(g) a(k)} dp dg dk + 
$e. | (p—q— k) {yx(q)yv(p) a*(k) — y(p)yx(g) a(k); dp dg 


+ 00] {ps — pops Op , o=VE+ mì, 
(2) W. HEeIsENBERG: Nucl. Phys., 4, 532 (1957). 
(3) K. L. NAGY: preprint. 
(4) M. Froissart: Suppl. Nuovo Cimento, 14, 197 (1959). 
(5) E. V. Troporovicu: Dokl. Akad. Nauk SSSR, 126, 1236 (1959). 
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where p,, py and a stand for the field operators of the V, N and 0 particles. 
© describes another fermion, with bare mass m. 4 and e are pure imaginary 


numbers. The commutation relations are the usual ones, except for yy, for 


which 
(4) {yv(P); pp} = — Sp —p'), 
which gives rise to an indefinite metric: 
NN a Neg Nae NS NN lo 


The real vacuum and N particle states are the corresponding bare states, 
the dressed V, g states have to be determined from 


(5) HA HER | He 


where | E) has the form 


(6) 1 —( by — Bg" 4 vx [otre (ak) 10) 
considering the Hamiltonian (3). (5) and (6) give 


Yo a 
(My — B)a — e,8 = —— — D(k) dk 
i À VAir) V20 ce) 


? 


— 
“1 
— 


Gin il 
(@ <= E) D(k) = € ——= 
Via V20 


x + (h—m)p 0. 


If H=w one obtains the scattering States, all of them have positive norm. 
We are interested in bound real particle states for which Z = ©. In this case 


{ Coa 
B = eae sao 
(3) m— H 
| oi) TRES 
[ Via V20 (© — E) 


‘and the eigenvalues have to be determined from 


(9) my— E 4 €0 = a] k? dk 
m 20 


(@ — E) 
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Using the identity 


1 it & 2? 
== eee ee oe 
© — Zz Oo oO | Do —2) : 
and defining the function H(z) i 
2 — My k? dk > 0 
(Oye EI IÈ i | = 
Jo 20(w —2) m—zZ 
È 1 1 
=a+be+ G(ez)e?+ 6 = h(z) + ¢e——_.,, 
m—e m— sz 
where 
my  [k?dk di "kh? dk i k? dk 
EB egg to a ol 
Jo 20? Jo 203 (Of: 2w*(@ — 2) 


the energy eigenvalues are given from (9) by the zeros of H(z). In this 
model a, b, ¢ and m are finite parameters, therefore 


CPP d 
Co oe MM, 
log @ log @ 


is the cut off energy. h(z) is the corresponding function in the Lee model. 
Considering (10) there are three different cases: 


H(z) H(z) 


+ +— Z + 75 
m Mg m Me 


Fig. La, b. — Discrete ghost case. 


H(z)- H (2) 
+ Z 
ml | pe i m Ma 
Fig. 2. — Tripole ghost case. Fig. 3. - Complex ghost case 
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For H; real, from (6), (8) and (10) one obtains 


Sa o dH 
(12) (E;|E> = Ô:lagol —| , 
de 
(A 
therefore in the discrete ghost case two of the states have positive, one has 
negative norm. This is the ghost state. In the tripole ghost case there is only 
one eigenvalue £, and one eigenvector. It has zero norm, since there dH/dz=0. 
| H,> is orthogonal to the scattering states having real eigenvalues and since. 
the metric is not degenerate (i.e. a state orthogonal to all state vectors of the. 
space, including the vector itself does not exist; cf. (1)) at least one state | D> 
satisfying (1) 


(13) (H — E,)|D> =|) 


/ 


necessarily exists. Indeed substituting (6) and 


(14) [D> = | vyy — dy* + vx fora") ax) 10) 
into (13) one finds 


A 0 1 
(My — Eyjy — 600 = ff o(k)dk + a, 
TT 20 


(15) EA ie D(k) 
Via vV20 i 
CV + (Ey a m)ù A, p E 


From (8) and (15) 


5 CA x 
m— E DE m= E)" 
x 


Yo 
o(k) = - —— | | | ; 
V 800 (wo — E) a m— E 


| 
(16) | 
[ 


Substituting (16) into the first equation of (15) and collecting the terms to 
one side, using (10) one obtains 


(17) go [Ey an 


which is equal to zero, showing that the dipole ghost |D) indeed exists with 
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an arbitrary x and y. From (6), (8), (14) and (16) an explicit calculation gives 


) < 2 AH | [ago |? d2H | 
(18) CH" 2719, a, et FAM 0 
mb 9 = i |Eo 


thus the dipole ghost |D> and |E,) are orthogonal. Similarly 


x 1 Ci: ay + ay” CPE ab. sp Ola | 
G9) dial, el qe, Glee! ae |, = 
SEE pe, 
6 de |, 


‘Therefore choosing |x|? equal to 
oD 


(20) Fee ay 


|D> can be normalized to +1. It is interesting to note that the normalization 
coefficient « of |H,> is determined by the normalization of |D>. Since |) 
and |D) are orthogonal, the supposition of a non-degenerate metric requires 
that a state |T> satisfying the equation 


(21) (H — E,)|T> =|D) 


exist. Substituting as before |7> and |D> 


(22) |T) = | epy — Àp* 4 vx [eu a(h) dk) |0> 


one obtains 


hg RE 
(my — Eye — A = al x(k)dk+y, 


Vin) v20 
(23) GIA 
o-E)x(k\=e— —+ o(k) 
eae VEL 
66 EE = mA =—O, 
from which 
2 Co i ia ia 
È = ale | m— E RA): 
«TESTE À eu 
(k) = —— | 1 + — le 
| te! V 8200 (0 — E) ( _o—#, 2a 
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Substituting these into the first equation of (23) and ordering it again as in 
(17) one obtains 
Eo 


so |T exists indeed with an arbitrary e and y. An explicit calculation gives: 


dH | iL Gel 
dz pe ' 9 de? 


(25) Jo (are 


/ | Ee * | * [go] d°H | Bay I yo) [Jo d'A | 
DTS = (yy° + 84° + xe") - 6 dss + (Nievo) 246 e: E 
logo |? d°H | 
1202 420 
By choosing ¢ appropriately one can therefore make 
(26) <i Oe 
Similarly 
= red |x|? d‘H 
r\n = |}! La | + PA A È 
| 6 de® [x 24 de* z] ? 


and now fixing y one can also prescribe 


(27) <TD 0 
Finally 
lun. |2 AS 
(28) HPP ee one ae 
6 — de} 


Here dH/dz|,, = d2H/de?| n, — 0 was taken into account. Thus fixing the para- 


meters x, y and e according to the equations (20), (26) and (27 ), respectively, 
we have 


EE} = <T|D = CBI Dy = LD) 0 


? 


These equations being results of an explicit calculation are, of € 


ourse, in: 
agreement with the general properties. 


Indeed, let us suppose that there exist: 
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three vectors, satisfying 


(H — E,)|Ey> = OF 
(H — E,)|D> =|H, , 


(H—E,)|T> =|D). 


In this case <H,|H,> — 0 since just this was the necessary condition for the 
system of eigenfunctions being incomplete. From the third equation 
CE, | (H— #,) | T> =H, | D>, from the first <7 | (H— H,)|H,>=0. Thus <#,|D>=0. 
But in this case <H,|7> 40 since the metric is non-degenerate. From the 
third equation <D|(H — E,)|T> = <D|D), from the second <T|[(H— H,)|D> = 
— (T|E,), from which <D|D) = <E,|T> #0. Now <T|D) and <T|T> can 
be chosen equal to zero without violating the requirement of the non-degen- 
erate character of the metric. 

Supposing that another state |Q> would be needed to complete the system, 
then 


(Ei Es). == |e >, 


and repeating the procedure as in (17) or (25) we obtain that such a state would 
exist only if 


d'A | 


dz* x, 


were also zero. In this model it is positive, showing that the system |M), 
|D>, |T> is complete, and the minimum polynomial in this subspace is 


MH) = HE 


Similarly it follows that |#,>, |D> and |T) are orthogonal to an eigen- 
state |Z if E E, (scattering eigenstates). So the projection operator to 
to this subspace has the form 


(80) P =|D><D| +[By><T|+| D> 


showing that one can meet negative norm only if a state contains both [o 
and |T). 

In the complex ghost case H(z) has two complex roots E, E' and a real 
root H,. E and E' are complex conjugated. The corresponding states have 
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the properties 
{E | E> == (H,| E" == {E|E) dl (E |E'y 2e 0, 
{E|E") #0 


as also follows from general principles. 

The probabilistic interpretation of this model is the same as that of the 
Lee model; in the discrete ghost case (°) and in the complex ghost case (7) 
they are not probabilistically interpretable since in the higher sectors negative 
scattering probabilities occur. The tripole ghost case is interpretable only in 
the absence of bound states with negative norms in the higher sectors. 


3. — Tripole ghost contribution to propagators. 


If, for the sake of simplicity, in the case of a simple relativistic scalar 
’ ’ 
(pseudoscalar) fiela there are discrete ghost states, the contribution of these 
states to the propagators is simply 


Joue) A(æ — x’; u?) du? 


but o is negative. In the dipole ghost case one finds the contribution (8) 


a 


fot) Aa—a'; up) dut + o(ut)(e —x'),0, Ae—-a';p°)dw, 


where 0, and 0, may be 20. 
Now we wish to find the tripole ghost contribution. The method applied 


here is the same as in (8), according to which we suppose: 


a) there exist à four momentum operator with the property 
(31) LESS p(x) | = i0,p(x) b} 


b) indefinite metric, with tripole ghost, 


c) and require relativistic invariance. 


The supposition b) in the spirit of the general principles and the preceeding 
model means: the eigenvectors of 1 do not form a complete system. To obtain 


(0) G. KALLÉN and W. PauLI: Mat. Fys. Medd. Dan. Vid. Selsk, 30, no. 7 (1955). 
(7) R. AscoLi and E. MINARDI: Nuovo Cimento, 14, 1254 (1959) 
(*) K. L. Nagy: Nuovo Cimento, 15, 993 (1960) 


. 
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such a set one has to consider all the states defined by 


2) (Pp) \o> = 9, <p|p>=090, pi—cpn, 9 p>, 


on 


aye) (Ep, )|Dp>—elp>, <Dp|Dp> = 0; 
(34) (P,.—p,)|Tp> = C..\Dp} lp Ep eas OS 


Here because of relativistic invariance 


(35) Oi), CUP O6, = 0,0. 
According to the indefinite metric 8 — +1. |p> and |Tp) are not orthogonal 
(36) (Tp |p’) = ayy, 


but |Dp> is orthogonal to |p> and | Tp). 
As in the general part, from (34) <Dp|(P,,— p,)|T'p> = 0, Dp | Dp>, from (33) 
<Tp|(P,,— P,)|Pp> = ¢,<Tp|p>, thus 


(37) g(p?)e = f*(p*)a* . 


The projection operator to this subspace is (cf. (30)) 


= HEC er Lee ul re 
(38) P= x{l2» a Poltlp7<Po]+]Tw) ae | 3 
and the contribution of these states to J’ 
/ /n | | x if / 5 ‘i / \ i / / \ 
(39) A'~ > fo |p(®) |Dp> — <Dp | p(x’) |0> + <O| p(x) |p> = Tp pl’) |0> + 


ab 
+ O|p(@)|2p> = Wo) es}. 
From (31) and (32) 

(40) <0| g(a) |p> = <O]@(0) [p> exp [ipa], 

from (31) and (33) 

(41) <0| g(a”) | Dp> = {<0| (0) | Dp> + ix,,¢,<0|p(0)|p>} exp [ipa] . 


and from (31), (84) and (35) 


(42) 0 y(«)|Tpy = Lo lp(0)| Lp) + 10 2, (<olpto |Dp> + 


+ 294% £0 | 90) ») exp [ipa 


26 - It Nuovo Cimento. 
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Substituting these matrix elements into (39), and requiring again relativistie 
invariance one obtains 
(43) Aa — x") = 0, A(w — e) + O,(@ — 2) 0, A(@ — #') + 


LOG) æ'),0,0, A(x — x), 


where (37) was also used. Each coefficient C; is real and may be negative. 
From (43) it follows that the Fourier transform of the tripole ghost contri- 
bution has the form 


! 1 À 112 2 Le Pa ll 2 ue il 2 è 
(44) 40) fot) pe + we Ry + fou (pi pa)? en fes (PERU) i 


DO 


each function 0, may be negative. 

The contributions from ghost of higher pole order can be calculated similarly. 
Here new and new derivates of A multiplied by the distances emerge. 

Since theories with discrete or complex ghosts do not seem to be inter- 
pretable, the general properties of a state vector space mentioned in the intro- 
duction and the probabilistic interpretability together seem to allow only 
pole-type propagators, i.e. the form (44) having as many terms as the « pole 
order » of the ghost. 

Nevertheless one has to be careful in applying the general results of a state 
vector space to the vector spaces of quantum field theories, which are con- | 
structed from the direct products of an infinite number of vector spaces. 


The author is greatly indebted to the Scuola di Perfezionamento in Fisica. 
Nucleare dell’Università, Rome, for financial support which made possible for 
him to work in the stimulating atmosphere of the Institute. 


IR ILA SS! UW INO) (>) 


Si discute una modifica appropriata del modello di Lee. Per un definito gruppo 
di parametri si trova, oltre allo stato fantasma dipolare, anche un fantasma tripolare. 
Usando il metodo di Källén-Lehmann si calcolano i contributi degli stati fantasmi tri- 
polari, che si suppongono esistere in una teoria relativistica del campo quantistico. 


(*) Traduzione a cura della Redazione. 
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Polarization Phenomena for Positron Annihilation-in-Flight (*). 


W. H. McMASTER 


Lawrence Radiation Laboratory, University of California - Livermore, Cal. 


(ricevuto il 23 Aprile 1960) 


Summary. — A method is presented for solving the Dirac equation for 
the annihilation of positrons using the Pauli spinors in such a way that 
all the polarization features are retained. Several specific examples of 
the polarization features are discussed, and finally the annihilation equation 
is put into matrix form for use with the Stokes parameters. This later 
representation is particularly convenient in that it readily points out 
the general features of the annihilation process. 


1. — Introduction. 


The computation of the differential cross-section for positron annihilation 
can be accomplished without the usual averaging over the positron-electron 
spins and summing (1) over photon polarizations by the use of the two-com- 
ponent Pauli spinors. Thus the wave function can be written in a form that 
greatly simplifies the calculation and leads to results that are given in terms 
of the positron and electron spin vectors (?). The final results of this compu- 
tation are similar to those obtained by PAGE (*), but are presented in a more 
general form. 

The computation is carried out in the center-of-mass system for a two- 
quantum annihilation. In this system the energies of the two photons, the 
electron, and the positron are identical. The geometry being discussed is shown 


(*) This work was performed under the auspices of the U.S. Atomic Energy 
Commission. 

(1) See W. HeITLER: Quantum Theory of Radiation, 3rd ed. (Oxford, 1954), p. 269. 

(2) H. OLsEN and L. C. Maximon: Phys. Rev., 114, 887 (1959). The author is 
grateful for having received a pre-publication copy of this paper. See also the Appendix: 

(3) L. A. PAGE: Phys. Rev., 106, 394 (1957). 
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in Fig. 1. For a transformation to the laboratory system one utilizes the stan- 
dard formulae: 


9/52 née) AA SE 
VE bs fos) © ci a DI 222400 


where f is given in the e-m system; # and e are the total energies in units of me? 
in the lab. and c-m system, respectively. 


Fig. 1. — Center-of-mass system for the two-quanta annihilation of positrons. The 
angles 9, and g, give the angle between the polarization vectors of the two quanta 
k, and k, and the plane of interaction (x-2 plane). 


2. — Computation of the cross-section. 


21. Evaluation of matrix elements. — The two-quanta-annihilation diffe- 
rential cross-section is given by (1) 

dae earn 

1) — = "||? 


where the matrix element responsible for the annihilation process is 


(uz )(u'*x34) (utazu”)(u'*x u) 
H— > i + È (ai 


e e! )) 


— 
bo 


and the summation is to be carried out over all spin directions and energy 
signs in the intermediate states. The quantity «,=a-e,, with e, being a 
unit vector in the direction of polarization ot photon k;, and, similarly, 
a = ae, for the second photon k,: The quantities e’ and e" are the total 
energies in units of me? of the electron in the intermediate states, where k:, 
and k,, respectively, are given off first. 
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Since the denominators are sensitive to the sign of the intermediate state 
we can use the relations H'u'= e'u' and H"u"= e"u" to rewrite eq. (2) as 


(wo, H'u')(u'*oœu) , (u*x;,H'u)(u"*æxu) 


(3) H=—y a Der 


€ ‘ € 


Momentum is conserved in the intermediate states, so we have p'—p—k, 
and p"=p—k,, where p is the momentum of the initial electron. Since 
k,=— k, and kj =k, =, we find e'?— 2e°(1—f così) and e"?— 2e2(1 +R così). 
Thus the matrix element becomes 


(4) 2e?(1 — B? cos? 0)H = 
— X.[1+fcos0)(ufaH'u')(u'*a,u) + (1—B cos 0)(ufaH"u')(u"*2,u)] 


The sum over the intermediate states can now be carried out in the usual 
manner, resulting in 


(5)  2e2(1— B? cos? 0)H = — [(1+8 cos0)(ufa,H'ayu)+(1— cos 0) (wpe, H'ou)] . 


To evaluate the matrix elements we will follow the spin formalism used by 
OLSEN and MAXIMON (?), which gives the result directly in terms of the spin 
vectors of the positron and electron. The initial (electron)—and final (po- 
sitron)—state wave functions satisfying the Dirac equation are 


i 
1)\i 
Uy = (= ) vo(S) » 
26 c‘p 
‘ese J 
(6) and 
, 1 
NAT 
U == | De | —o-p v(— P), 
e—1 


where v,(S) and v(— P) are the Pauli spinors, with S and P being the electron 
and positron spin vectors, respectively. The positron’s momentum is repre- 
sented by — p since it is moving in the negative 2 direction. 

Using these wave functions, we have 


Posse: 1 
? (e2 — 1)? JE op I ) 
(7) (fa H'asu) = — È dl sti [o (a: p'+ b)x,] CR OS 
Re ES 
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and similarly for (usx,H'au) where we have used the Dirac Hamiltonians: 
H'=a:p'+ 6 and H’=a:p’+ 68. 
Multiplying out the terms in eq. (7), one finds 


2pe(ufaH'asu) = 0 [(e —1)o-po-e,o'p'o-e,— (e + 1)o-e,0-p' G-e,0°p — 
—p?o-e,0°'e,—O'po-'e,o'e,a'p|v, 
and similarly: 
2pelusH'au) = vy [(e—1)o:po-e,o:p'a-e,— (e+1)0'e,0:p'G-e,0°p — 
— p?c'e,c'e, — O° po'e,o'e,o'p|v. 
Combining, to obtain the complete expression for the matrix element, we find 
(8)  4pe?(1— PE cos? 0)H = v%{[(£+1)(c-e,c-p'o:e,o:p+0:e,c'p'6‘e,s-p)— 
—(e—1)(o-po-e,o'p'c-e,+ 6'pa-e,o'p'o-e;) + 


+ p(o'e, c'e + 6:e,6°e,) + 


+ (o:'po-e,o-e,0'p+oa'poa-e,o:e,c:p)| + 


+ B cos Ü[(s + 1)(6-e,6'p'c-e,0:p—G-e, 6: p" oe, cp) 
—(e—I1)(o-po-e,o'p'c-e,— o'po-e,o'p'a-e,) + 
+ p?(G-e,0°e,—G:e,0°e,) + 
+ (o-po-e,c-e,0:p—o'po-e,o:e,o'p)]}v. 
In multiplying out the terms in eq. (8) we make use of the relation 
(c- A)(c-B) = A-B+ioc-(Ax B) 


which follows from the properties of the Pauli spin matrices. 


The result can be put in a simple form if we define a set of functions de- 
pendent upon the photon polarizations as 


0 = exp [dò] sin y, cos p, — exp [id,] cos g, sin P15 
01 = exp [dd] sin y, cos y, — exp [id,] cos g, sin Pi; 
(9a) 

0, = COS p, COS Y, + exp [1(d, + 6,)] sin Pa SiN Gy , 


03 = COS Pz COS p, — exp [1(d, + d.)] sin g, sin Pi + 
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where we are representing the electric vector of the photon by e = cos per + 
+ sing exp [id]e,, with e, and e, being unit vectors in the plane of inter- 
action and perpendicular to it, respectively, and 6 is the phase angle. For 
linear polarization (6 = 0), we see that eqs. (9a) reduce to 


Qo = Sin (p2— 1) 

01 = SiN (pr + a) 

0, = COS (P2— gu) = (er es) 

_ [2(pre.)(prez) 


03 = COS (Pr + Pr) TOGO) (e. :e.)|. 
With these conventions, eq. (8) reduces to 
(10a) e(1 — f? cos? 0)H = (H+ iH-o)v, 
where 
= È (0. + 03 Sin? 6) 
=— 0, B sin 0 
(10b) 


I 


Hy 
H, 
H, = 0,f sin 6 così 
H, 


I 


92. The cross-section in general form. — Upon substitution of the matrxi 
elements of eq. (10) into eq. (1) we find the differential cross-section 


do re 1) | | 
dQ Sz 1 e2(1 — B2 cos? 0)? {o*(H, — iH* +o) vot (Ho + iH-6)v} ; 


(11) 
(12) 00, = 3(1+S:0), 


however, is the projection operator for the electron. 
Therefore, the term in brackets in eq. (11) becomes 


(13) 4{0*(H, — 1H*-o)(1 + S-o)(H, + Ho). 
Upon expansion we utilize the properties of the positron spinors, 7.e., 


(ei == il and © 
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where P has the components P,=—sinycosw, P,=—sinysinw, and 
P,= cos y (see Appendix, eqs. (A.14)). Thus the differential cross-section is 


given by 


Aa eee Se rr HI RI ae Ss) (Ss He 


+ (S-H*)(P-H)]— iH, [(S-H*) — (S-H)]— 1H, [(P-H) — (P-H°)] + 
+ i[S-(Hx H*) + P- (Hx H*))— H,[(Sx P)-(H + H*)]. 


For the special case of plane polarization, where H* = H and A, =e 
eq. (14) reduces to 


5 les i! | 
(15) FG et — B® costo) 
-{Hî | H:— (H?— H')(P-S) — 2(S-H)(P-H) — 2H (Sx P)-H}. 


2°3. Special case of longitudinal spin. — For the special case of longitudinal 
spin along the ¢ axis it is instructive to rewrite eq. (10a) in the form 


(16a) (1 — 82 cos? 0)H dee IL) 
a e(1 — Bb? cos? =% D 
i Hy, Hogs i 
where 
Le î 
Hi = (105 Bos Pos sin: 0) , 
1 3 : 
(160) xs: (— i00+ Bor + Bos Sin? 0). 


A, 
H,, =— B sin 0(i0, — 05 cos 8) , 
H,, =— f sin Ü(io, + 03 cos 8) . 
In this cage the spin spinors have the simple forms: 
= | 
o = 10: electron spin along the positive 2 axis 


FES 
o = 01: electron spin along the negative e axis . 


0 
= ( : positron spin along the positive 2 axis 
ss | a positron spin along the negative e axis . 


Using these spinors, one readily finds that the matrix elements of eq. (16a) 
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refer to the following alignment of electron-positron spins: 


Ay: +—, (14) 
Ha: <+, (4) 
He Py) (tH) 
Hy: ——, (44) 


This form of the cross-section, then, is similar to that presented by PAGE (5). 


2°4. Photon polarization. — The cross-section depends upon photon pola- 
rization through the o-functions defined by eqs. (9a). For a complete de- 
seription of the polarizations of the two photons k, and k, we can specify the: 
o’s in either a set of plane-polarizations tates or circular-polarization states. 
From the definition of the o0’s one sees that for circular polarization 0, and 
0, are pure imaginary. To avoid confusion it is convenient to rewrite the 
matrix elements for the circular-polarization de- scription as 


= a (02 + Q3 Sin? 0) 


€ 


; He 10,6 510 0 


(17) 
H, = 03} sin 0 così 
H, =— 1Q/e 
and 
il 5 
Ay, = E (004 Bo Bos sin? 0) tit) 
H Sn 00 + Pos + Pos sin? 0) : (44) 
(18) Log — E Yo T PO2 T PGs 5 : Ÿ 
Ho a = sin 6(0, — 03 COS 9) : (tt) 
H,, = — B sin 0(o, + 03 Cos Ü) : (4) 


For the plane-polarization description, let p and s refer to polarization in. 
and perpendicular to the plane of interaction, respectively. For the circular- 
polarization description, let R and Z indicate right and left cireular polari- 
zation, respectively (4). With these conventions, the values of the os for each 
set are given in Table I. 


(4) In this representation, e.g., right circular polarization of both photons is spe- 
cified by the vectors e,—3(e,—ie,) and e,=43(e,+/e,); thus each photon advances 
as a right-hand screw. . 
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Tage I. — Photon polarization states. 
Linear Circular 
State en = n = 
pp ss ps | sp RR LL RL LR 
| 

valo 0 0 1 =| 1 =a 0 
01 0 0 1 1 0 0 =| 1 

0s 1 1 0 0 1 1 0 

| ees 1 =i fire RIO 0 0 1 1 


3. — Discussion of the cross section. 


3°1. Unpolarized electrons and positrons. — The differential cross-section for 
unpolarized initial and final spins is obtained directly from eq. (14) as 
do~ (H?+H?) or from eq. (16b) and (18) by averaging over the spin combi- 
nations; i.e., do }(H?,+ H}:+ Hi,+H;,). The result is 


do gi ih 


a) dQ  8Be (1 — B? cos? DEA 


1 
da [oo + P*(0» + 03 sin? 0)?] + f? sin? 9(0° + 03 cos? 0} 


for both the linear—and circular—polarization descriptions. From this equa- 
tion it is seen that in the low-energy limit (6 — 0) the photons tend to be 
cross polarized (ps and sp). In the circular-polarization description the photons 
are RR and LL, i.e., there is no net angular momentum. 

If the polarizations of the two photons are not observed, then one obtains 
the differential cross-section by summing over all directions of photon pola- 
rization as given in Table I, and obtains 


do tè [1+ 2? sin? 9 — B4— B4 sint 6 
dQ  4fe (1— B? cos 0) i 


As a check on eq. (19) we can substitute the vector forms of 0, and 03 as 
given by eq. (9b) and make use of the fact that gg =1— 0? and po? = 1 — gf. 
Then eq. (19) becomes 


do o, ro ET Co. 9)(e,-e,) + 4(p-e;)(p: es)(€," e») 4(p-e,)?(p-e,)? 
dQ 8fe (e — p? cos? 0) : (e2— p? cos? 0)2| ? 


which is Heitler’s formula. 
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Let us how investigate the possibility of producing a polarized y-ray beam 
using an unpolarized positron source and an unpolarized electron target. The 
experiment would consist of measuring the polarization of the lower-energy y 
in the laboratory system, which would then determine the degree of polariz- 
ation of the high-energy y. We will define the polarization correlation through 
the equation 


OO; 
pra I 
O11 9 
where o ,—0,,+0 gives the cross-section for cross-polarized y’s, and o, = 


= Oy + 0, gives the cross-section for parallel polarization. Thus, using eq. (19) 
we have 


ZII 
1428? sin? 0 — f4— BA sint 0’ 


and for the two limiting cases: 


_1--2*1— 2) _1-P? 
Pur) 1+ 2p — fa) and POy= 14° 


Thus we see that at high energies (8 —1) the cross-polarization correlation 
increases from zero at 0=0 to complete cross polarization at 7/2. This does 
not seem to be a good source of polarized annihilation radiation since the cross 
section. at 7/2 is quite small and the value of P does not increase above ~ 0.1 
until 0 is greater than about 30° for energies of interest. 


3°2. Longitudinal spins. 


302.1. Plane-polarized photons. — The cross-section for antiparallel 
longitudinal spins is given directly by eq. (15) as do > 2(H;+H?), or can 
be obtained from eqs. (16b) by averaging over the two spin orientations; 7.¢., 
do }(Hi,+H;,). The result is 


do (e 1 Lia à E 
i Du , + 0, sin? 0}. 
AL (a). pi pe II 


Summing over photon polarizations, one has 


DI do ls (1 iù | 24 2 Q1 4 
eo) (zo), 2Be? (1 — B? cos? 0)? Cons LI, 
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Similarly for parallel spins, we have from eq. (15) that do 2(H° — H°) =: 
— 2(H?+H?), or from eqs. (168) that do > }(Hj,+H;,). The result is 


do #6 6? sin? 0 at, 
CE (3) = * dpe? (1 — B? cos? 0)? KE Eck De 


Summing over photon polarizations, one has 


(211) (7) i rè B? sin? 6(1 + cos? 0) 


dQ 2Bez (1 — f? cos? 6)? 


From these equations we see that at low energies the annihilation takes 


place with spins opposed, while at high energics annihilation with parallel 
spins predominates. In both cases the y-rays tend to be cross polarized. Note, 
however, that for 0 — 0, (do/dQ),,,=0. This feature has been used exper- 


imentally to measure the longitudinal polarization of positrons in 8-decay (°). 


3°2.2. Circular polarization. — For the case of antiparallel longitudinal 
spins, the cross-section as derived from eq. (14) is do ~(Hj+H?+2iH,H,) = 
= (1/e?) [+ @o+ (02+ 03 sin? 0)]?, where the upper sign refers to the case where 
the positron spin is in the + direction, and is just H°, of eqs. (18). The four 
cases are then 


do (1 IL a B)? 


do,,~ (LE B?)( se B)? i 
do,, ~ (1 — 6?) 6? sin‘ 0, 


do, (1 = B*) B* sinto 


where the upper and lower signs refer to positron spin in the positive and nega- 
tive 2 directions, respectively. For positrons polarized in the positive 2 di- 
rection the ratio of right circular polarization (RCP) to left circular polariz- 
ation (LCP) for photons k, and k, is given by 


RCP (1—f)?+ 2 sint6 


LCP (+8)? + B sin® 6’ 


from which we see that at low energies (8 +0) the photons are unpolarized, 


while at high energies these photons are mainly LCP. The reverse is true if 


(5) See for example: S. FRANKEL, P. G. HANSEN, O. NATHAN and G. M. TEMMER 
Phys. Rev., 108, 1099 (1957). 
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the positron is polarized in the negative 2 direction. These results are illus- 
trated in Fig. 2 where we see that photon k tends to pick up the spin of the 
particle whose momentum is nearly parallel to k. 


1 k, K, 
LCP LCP RCE RCP 
Où 
6 
S < ap p 
LCP LC RCP RCP 
r r 
Kr D Ko Ko 
(a) Positive positron spin, H,, (b) Negative positron spin, H,, 
Fig. 2. — Cireular-polarization correlations at high energy for antiparallel-longitudinal 


spin. Case (a) with positive positron spin corresponds to matrix element H,,, and 

Case (b) corresponds to H,,. Two values of 0 are illustrated: one less than 7/2 and 

one greater than 7/2. In practice k, and k, are undistinguishable. The solid arrows 
designate the spin directions of the electrons and positrons. 


For the case of parallel-longitudinal spins, the cross-section as derived from 
eq. (14) is do (H°+H° +2iH,H,)= f? sin? 0(0, + 03 cos 0)?, where the upper 
sign again refers to a positron with spin in the positive 2 direction and is just 
Hi, of eqs. (18). The four cases are then 


(24) do,, =? sin? O(1 + cos 0), 

do,, = 6? sin? 6(1 += cos 0)’. 
For positrons polarized in the positive z direction, the ratio of polarizations of 
photons k, is given by 


ROB “(1 
POP 


cos 0)? 
cos 0)? ? 


[al 
MI 


and for photon k,, by the inverse. Since k, for 0 < 2/2 and k, for 0> a/2 
are undistinguishable, if we look at the photon for angles 0<0< x/2, then 
the ratio at high energies becomes 


RCP (1+ così)? 


125) LCP (1—cos6)?’ 


and we see that these photons are predominately right circularly polarized. 
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These results are illustrated in Fig. 3, where again the photon tends to pick 
up the spin of the particle whose momentum is parallel with the photon. 


k, 
RCP 
ID 
LCP RCP RCP LCP 
,’ ! 
k, Ko k; ky, 
(a) Positive positron spin, H,» (b) Negative positron spin, H,, 


Fig. 3. — Circular-polarization correlations at high energy for parallel-longitudinal spim 


Many times an experiment will be conducted with just one of the particles 
polarized. We will discuss two such cases (5): a) A polarized positron beam 
incident on an unpolarized electron target. This type of experiment might 
be done to study positron polarization in £-decay when parity is not conserved 
b) An unpolarized positron beam incident upon a polarized electron target, 
where the electrons are aligned by a magnetic field either parallel or anti- 
parallel to the positron’s momentum. In both cases we will limit observation 
to the higher-energy photon in the direction of the positron’s momentum. 

For the case of the polarized positrons incident on unpolarized electrons, 
let the positrons be polarized in the negative e direction, which covers the 
case of polarized positrons from 8-decay. Then the cross-section is given by 
eq. (14) as do Hj+H?+2iH,H.— 2iH,H,, or from eq. (18) as do H?,+H},: 
that is, 


do ~ (1 — B?)(@) + Bos + Bos sin? 0)? + B? sin? 0(0, + 03 cos 0)? . 
The four cases are then 
dos (1 — 6?)(L-— Bes 


Goes lee ie ae 


doz, (1 — B?)B? sint 6 + 6? sin? 6(1 — cos 0}, 


do;,~ (1 — B?)B? sin‘ 0 + B? sin? 6(1 + cos 6)2 . 


Observing the high energy photon means we are restricting ourselves to ky 
when x/2<60<x, and k, when 0<0<%x/2. Let us then redefine the angle 0 


+ 
ri 
co 
N 
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as being the angle between the direction of motion of the positron and the 
observed photon (0<0< 7/2). The polarization ratio is then determined 
from eqs. (26) as 


Rit Pe = pi) Cb e) 6? sin! 0] + B? sin? 6(1 + cos 0)? 


27 3 
(27) LCP (1—-p)[(1— pb)? + 6? sin‘ 0] + B's sin? 6(1 — cos 0)?’ 


from which we see that at all energies the observed photons are predominantly” 
RCP. This process is shown in Fig. 4a. 


RCP Lee 
pera ex 
of SOS 
(a) Polarized positrons, (b) Unpolarized positrons, 
unpolarized target polarized target 


Fig. 4. — Circular polarization produced by high-energy positron annihilation-in-flight 
when one of the particles is unpolarized. 


When the target electrons are polarized in the positive ¢ direction and the 
incident positrons are unpolarized, then the cross section is given by eq. (14) 
as do Hi+H’+ 2iH,H,+-2iH.H,, or from eqs. (18) as do~ Hebe Soli. Taney 


do ~ (1 — B?)(0) + Bos + Bos sin 20)? + p? sin? 6(0, — 03 cos 0)? . 
The four cases are then 
do er (1 — Bat ap) 3 
do,,~ (1—f?)(1 —B}, 
(28) 
do,,— (1 — B?)B? sin! 0 + B? sin? 0(1 + cos 0}, 
do,,~ (1 — B?)B? sin! 0 + B? sin? 0(1 — cos 6)* 
Again we shall restrict ourselves to observing the high-energy photon in the 


lab system, and thus redefine the angle 0 as before. The polarization ratio 
is then determined from eqs. (28) as 


ROR (p51 (la BE + B? sin 0] + f? sin? 0(1 — cos 6)? 
LOP (1 —#?) [(1+ 6)? + 6 sim 0] + B? sin? 9(1 + cos 6)? i 


(29) 


from which we see that at f =0 the photons are unpolarized, at high ener- 
gies they are LCP, and for an intermediate range they are RCP. This process. 
is shown in Fig. 4b. 
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33. Transverse spins. 


3°3.1. General. — Let us now consider the case where the electron and 
positron are polarized transverse to the direction of motion at azimuthal 
angles ©, and ©, from the x axis. Then, recalling eqs. (A.14), we find that 
eq. (14) becomes: | 


do 
dQ 


we 


(30a) H,+H?+(H?— H*)(cos @, cos 0,-+sin ©, sin w,)+2H? COS ©, COS + 


+ 2H? sin w, sino, — iH;(H% — H,)(cos w, — cos ws), 


where we have utilized the values of o given in Table I to eliminate those terms 
that are alway zero. Thus, 


ABe?(1 — b? cos? 0)? do Di; , © + © 
(30b) ri 2e Jon TEN CORE (25> + Ecos? or | 


Da iH, : 
+ H? sin? | | = °(H*— H,)(cos w, + cos ws) : 


di 


where the last term is non-zero only for the circular-polarization description 
when it becomes (?/e)0,03 sin® 0(cos w, + cos w,). If one averages over photon 
spins, then only the first four terms remain and lead to Page’s (?) formula 


9 


did 


(ET MU (Oy eo 
+ Be? sin? 0 [cos (rn) + così 0 sin? (as) : 
dd 


d (0a 5 a Di — Mo 
(31) TS gin? ( a | + B2(1 + sin‘ 6) cos? — =) Da 


a 


For antiparallel spins, ©, —«, = x, we have 


(32a) do per. IRA 5 o 
2 a (1— f?) + A? sin? O[sin? w, + cos? 0 cos? w,], 


and for parallel spins, ©, — w,, we have 
(32b) da ~ B*(1— B?)(1 + sin’ 6 2 sin? 0 [cos? :08? 0 sin? 
TO ) + 8? sin? 0[cos? w, + cos? 0 sin? w,]. 


From egs. .(32) we see that in the low-energy limit, 6 =0, the annihilation 
process occurs only for antiparallel spins, which agrees with that found for 
longitudinal spins. 

A close look at eq. (30) shows that the forward y-rays are unpolarized. 
The only polarization that can occur arises from the H? term, and in this case 
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the polarization. is found to be in the plane of interaction with a maximum 
at 9 =a/2. This polarization is a maximum at low energies (6 40) and de- 
creases with increasing positron energy. 


83.2. Polarization correlation when one particle is polar- 
ized. — Let us consider the case where one particle is polarized in the plane 
of annihilation, © — 0, and the other particle is unpolarized. Then from 
eq. (30a) we find that 


do 
dQ 


{ 2 q * 
~ H°+H iH, (Ht— H.), 
which states immediately that there is no plane polarization correlations, 
since H; — H, in this description. If we conisder cireular polarization, how- 
ever, then we have 


do , 
(33) go Ho + H°+2iH.H,. 
From this equation we find that the maximum 
polarization effect occurs for 0 — x/2, and if we Gee 
define the degree of polarization in the usual 
manner as 
do,— do, 
34 (Pap i | 
i do, do” Spin 
€ * -—@< P 


where do, = 4(doy»+ do,,) and do, = $(do,,+ do,,), Hits te isa 


then one finds one of the particles is tran- 
sversely polarized. The y-ray 

D 26 given off in the direction of 
(35) i POSE 26? (1— B)] ? spin has the same helicity. 


which states that the y-ray is predominantly RCP and thus partakes in the 
helicity of the polarized particle. This situation is shown in Fig. 5. 

On the other hand, if one of the particles is polarized perpendicular to 
the plane of annihilation, no polarization correlation exists. 


4. — Matrix representation of the cross section. 


If one uses the Stokes parameters to describe the polarization states of the 
electron and positron, then it is possible to write the annihilation cross-section 


27 - Il Nuovo Cimento. 
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in matrix form (5). This form of the cross-section is very convenient for the 
description of the polarization phenomena in that the general features of the 
interaction are readily determined. 

The initial particles are described by the two four-vectors 


1 / 


1 

1 S, ii P, 

(36) (5) = 8, and (5) = = 
E 


y 


for the electron and positron, respectively, where positive values of the com- 
ponents refer to spin in the positive z, v, and y directions. With this notation, 
eq. (14) can be written in the particularly simple form: 


do = Il 
(37) io £34 (5). 
where 
pi 

38 MESIA ASIA; 
55) She? (1 — p? cos? 0)? 

re DU, Jeske) 2H, 0 

QU oat HsH,) HA vd thee it 0 1) 
MH Ho 0 Hi He 22H} 0 
0 0 0 Hj— H?+ 2H?, 


for the description of circular-polarization states of the annihilation quanta, 
and 


; 7 
a) a 8Be?(1 — pi cos? 0)? 
Hi + H? 0 0 0 
0 eS Er 2H.H, 0 
0 HE TA H?— H*+2H? 0 
i 0 0 Hi Hans 


for the description of plane-polarization states of the annihilation quanta. 


(5) W. H. McMaster: Matrix Representation of Polarization, University of Cali- 


fornia Radiation Laboratory Report UCRL-5496, 1959. (To be published in Rev. 
Mod. Phys.). 
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APPENDIX 


Spin formalism. 


The Dirac Equation, (a:p+f —e)u=0 can be written in matrix form as 


(A.1) om Se UE. 


Now let us write the four-component wave function in a form in which the 
two-component Pauli spinor v appears explicitly; 


v 
(A.2) = IN () È 


where N is a normalization factor. ‘With this notation, the Dirac equation 
reduces to the two equat ons 


(o-p)y+(1—ev =0, 


(G-p}v—(1+e)y —0, 


from which we find y = (o-pjv/(e+1). The normalization, u*u — 1, determi- 
nes the value of N, and we obtain our wave function in the form 


(A.3) UE ( È + do v(S). 


The Pauli spinor, v(S), satisfies the spin equation for electrons: 
(S-a)v(S) = o(S), 


whose solution specifies the spin direction of the electron through the equation 
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where y is the polar angle and © the azimuthal angle. Thus a particular re- 
presentation would be 


COS exp [— im] 
2 


AA % = US) = 
| sin x 
: 9 


from which we find the particularly useful forms 


il 0 
v= (0) and % = (1) 


for spins in the positive and negative directions, respectively; that is, the 
electron is said to have a positive longitudinal spin (right helicity) when the 
spin and momentum vectors are parallel. 

The Pauli spinor can be written in forms that display the three compo- 
nents of the spin explicitly; 4.e., 


1+8, | 


(A.5) OS I or, (È + iS, 


for positive spin directions, and 


v(5) = [21 — S)1* & a 


ha 


for negative spin directions. In addition to the normalization equation 
vi =1, one obtains the projection operator 

A 5 * 
(A.6) Vodo = (1+ So) 


from the adjoint equation 


From these spinors it is found that 


(A.7) 0.60; ie 


The negative energy state of the electron, 4.e., the positron, satisfies the 
Dirac equation with the signs of the energy, momentum, and spin changed (7) 
i.e 

22 


b) 


(A.8) (—a-p+ p+ eu— 0, 


(7) 8. S. Scuweper, H. A. BetHE and F. DE Horrmann: Mesons and Fields, vol. 1 
(New York, 1955), p. 142. 
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which has the solution 


SNS 
(A.9) u = ( | o-p v (— Se 


where the Pauli spinor now satisfies the spin equation 


(A.10) (S-0)0—S)= = v(—S). 
This leads to the result that 


(A.11) v*ov =—S, 
where again #, is positive when it is parallel to the momentum. For the case 
of positron annihilation under discussion, we have chosen the positron to move 
in the negative 2 direction; therefore, we must change the sign of p and S 
above so that eqs. (A.9) and (A.11) become, 


alt 
e—1\t 
(A.12) n= | = | SI »(S), 
e—1 
and 
(A.13) ose 


The solutions of the spin equa ion are then given by eqs. (A.5); we must now 
recall, however, that S,=1 means spin in the negative 2 direction. Thus we 


find that now 
Cir = D and Va = 1 
TAI ‘ AN 


are the particular solutions for positron spins in the positive and negative 
2 directions, respectively. 

Since most of our interest will center around longitudinal spins, let us now 
respecify the positron’s spin vector so that the sign of the longitudinal spin 
Mull be explicit; 4.¢., PS, "1.6. 


P,=— sin y cosa, 
(A.14) P,=-—sinysino, 
Pe = COSY. 
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Equation (A.13) then becomes 
(A.15) viov =—P. 


ESS SuU ANG Ome) 


Si presenta un metodo per risolvere l’equazione di Dirac per l’annichilazione dei 
positroni, usando gli spinori di Pauli in modo da conservare tutte le caratteristiche 
della polarizzazione. Si discutono diversi esempi specifici delle caratteristiche di pola- 
rizzazione, ed infine l’equazione di annichilazione viene scritta in forma di matrice 
per usarla con i parametri di Stokes. Quest'ultima rappresentazione è particolarmente 
utile in quanto rivela le caratteristiche generali del processo di annichilazione. 


(*) Traduzione a cura della Redazione. 
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Non-Spontaneous Disintegration of Hypernuclei in Flight (*)(**). 


H. NEUMANN and H. S. VALK (**) 


Department of Physics, University of Oregon - Eugene, Oreg. 


(ricevuto il 26 Aprile 1960) 


> 


Summary. On the basis of a binary model, expressions are obtained 
for the probability of hypernuclear disintegration in flight. Numerical 
calculations for 50 MeV *He, slowing in liquid hydrogen, oxygen, and 
xenon indicate that such processes are most prevalent for medium-Z 
media. Of the two processes considered, stripping and electric disso- 
ciation, the former appears to dominate in the medium- and large-Z range. 


1. — Introduction. 


The expansion of programs to observe hypernuclei in bubble chambers 
makes it interesting to study theoretically some of the less frequent hyper- 
nuclear processes. The recent, successful treatment of many hypernuclei as 
lightly-bound, binary associations (1) of a A,-hyperon and a nuclear core im- 
plies that such associations can be disintegrated in flight. As in the case of 
the deuteron, the two dominant disintegration processes are stripping and 
electric dissociation. These processes both produce free hyperons, whose sub- 
sequent spontaneous decays are observable. The tracks of the freed hyper- 
nuclear cores, which will retain most of the initial kinetic energy of the hyper- 


(*) A preliminary report of this work was given at the Los Angeles meeting of 
the American Physical Society, H. Neumann and H. 8. VALK: Bull. Am. Phys. Soc., 
Ser. II, 3, 402 (1958). 

(**) This paper is based on a thesis submitted to the Graduate School, University 
of Oregon, by H. NEUMANN in partial fulfilment of the requirements for the degree 
of Master of Arts. 

(*+) Temporary address: National Science Foundation, Washington 25, IDEE 

(1) R. H. Darirz and B. W. Downs: Phys. Rev., 111, 967 (1958). 
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nuclei, should also be visible. In this paper we obtain expressions for the 
cross-sections for these phenomena. In order to estimate the efficiency of 
various bubble chamber liquids in producing such events, we consider in detail 
the example of a 50 MeV *He, slowing in liquid hydrogen, oxygen, and xenon 
and compare its rate of disintegration in flight with its spontaneous decay rate. 


2. — Binary hypernuelear model. 


Among the possibilities for the origin of the A,-nucleon force are 4) the 
exchange of two pions (corresponding to a Yukawa potential with a range 
of .71fermi), b) the exchange of a K-meson (range, .40 fermi), and ce) the 
exchange of combinations of K-mesons and pions. In the present paper we 
employ b). Since we limit our attention to low-energy disintegrations (within 
an energy range for which an effective range theory is valid), the choice is 
not critical. 

For definiteness we treat 5He,, although the methods can be generalized. 
The hypernucleus *He, is one of the most frequently observed (?), and the 
tight binding of its core is especially favorable to the binary assumption. 

If we assume that both the triplet and the single A,-nucleon potentials 
have the same shape function, v(r), the average potential for the A, over 
A — 1 core nucleons is (1) 


(1) u(R) = voir — R'|) 0(R'/R,) dR’ , 


where V is the volume integral of the A, interaction with the core and o(R'/R.) 
is the density of nucleons in the core. According to electron scattering exper- 
iments (*), the density function can be taken as a Gaussian, with root-mean- 
Square radius R, — 1.44 fermi, for ‘He. However, because the A, is presu- 
mably in the hypernucleus in a (1s) orbital, if should bind the core more tightly 
than in free ‘He. The effect is to reduce R, by a small amount, taken to be 
six per cent (1). To evaluate (1) the Yukawa A,-nucleon potential is replaced 
with an equivalent Gaussian, by use of effective range theory (4). The inter- 
action u(R) thus becomes a Gaussian, with r.m.s. radius 1.54 fermi. 

For the electric dissociation calculation, it is more convenient to take the 
interaction between the A, and the core in the form of a square well. 


2) W. F. Fry: Ann. Rev. Nucl. Sci., 8, 105 (1958). 

R. Horstaprer: Ann. Rev. Nucl. Sci., 7, 231 (1957); 

a dj i BLarr and V. F. WrISSsKOPF: Theoretical Nuclear Physics (New York. 
P- 
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Applying effective range theory again, we find that the radius of the 
equivalent square well is R, —1.80 fermi. The binding energy of the A, to 
the ‘He core is taken to be 2.9 MeV (!). From the usual two-particle square 
well quantum condition, the well depth is then found to be V, = 28.0 MeV. 


« 


3. — Stripping. 


The most obvious nuclear reaction that leads to the production of free 
hyperons is stripping, in which the hypernuclear core strikes a target nucleus 
and the hyperon continues on. The reverse situation, in which the hyperon 
strikes the target and the core does not, will probably lead to another hyper- 
nucleus. If the entire hypernucleus collides with the target at the high kinetic 
energies of interest in the present paper, the target will become highly excited. 
If a hyperon is torn loose in this last process, it will probably knock out a 
fragment from the target and become attached to it (5). 

To obtain an estimate, we assume that the stripping cross-section o; is the 
product of the probability P that an incident core strike the target while the 
hyperon miss and the cross-section o, for compound nucleus formation by 
the core and target, so that os = Po: This semiclassical model (°) really re- 
quires that the collision take place in a time short compared with the oscil- 
lation period of the hypernuclear system. Since this condition is not well met 
at the kinetic energies considered here, we assume that our results give only 
an upper limit for the cross-section (*). 

The probability P can be computed from geometrical considerations and 
is given by 
(2) P =1 + (ofxa) [1 — (0/2a)2]? — (2/7) sin! [1 — (@/2a)?]' — 2o/xa , 
where a is the sum of the target and core radii and o is the projection on the 
plane perpendicular to the projectile path of the vector R connecting the two 


segments of the hypernucleus. Let Y(R) be the wave function for the binary 
hypernucleus. Then the probability averaged over all projections is 


(3) P= am fe 


|? de do = ala) | Ree aR = R2a, 


(5) Such a mechanism seems to have been involved in the production of many 
of the hypernuclei that have been observed: E. M. SILVERSTEIN: Suppl. Nuovo Cimento, 
10, 41 (1958). 

(6) K. MarHer and P. Swan: Nuclear Scattering (Cambridge, 1958), pp. 354-359. 

(*) We also neglect small diffractional effects such as those considered by A. I. 
AKHIESER and A. G. SITENKO: Phys. Rev., 106, 1236 (1957). 
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where È is the mean separation of the hyperon and core (for *He,, computed 
to be 2.24 fermi). The location of the core relative to the target is undefined 
to approximately À, the reduced wavelength of the core relative to the target. 
The estimate for P can be corrected approximately for this effect by replacing 
a with a+/. 

The cross-section for compound nucleus formation by ‘He projectiles on 
various targets is given by BLATT and WHISSKOPF (?). 


4, — Electric dissociation. 


The general theory for the electric dissociation of binary nuclear systems 
has been developed by DANCOFF (5), MULLIN and GUTH (°), and others (1°). 
Here we follow MULLIN and GUTH and consider an interaction energy 


(4) Vez 


where 7) is the sum of the r.m.s. radius of the hypernucleus (for *He,, com- 
puted to be 2.71 fermi) and the ordinary radius of the target nucleus, eZ; is 
the charge of the i-th hypernuclear segment, eZ is the target charge, r is the 
vector from the hypernuclear mass center to the target, and r; is the vector 
from the hypernuclear mass center to the i-th hypernuclear charge. The inter- 
action is cut off at 7, to exclude direct interactions of the type considered in 
Section 3. 

Since a multipole expansion of V converges rapidly, we retain only the 
dipole term, so that the cross-section per unit disintegration energy for electric 
dissociation of a binary hypernucleus described by a square well potential, 
as in Section 2, is (°) 


A do 
(5) 1 = 82[u(Z./M.)(Ze?/hv,) ? BFF, exp [2a(n, — m:)] 7 
(*) J. M. Buarr and V. F. Weisskopr: Theoretical Nuclear Physics (New York 
1952), p. 353. : i 
(*) S. M. Dancorr: Phys. Rev., 72, 1017 (1947). 
(9) C. J. MuzuN and.E. Guru: Phys. Rev., 82, 141 (1951). 
fe VEE MAMASAKHLISOV and G. A. Cumasuvrir: Soviet Physics TAB AL PES AGO 
JO). 
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where 
B; = f[2(%, + ky)? | ae f[2(%, TR k,)ro] ; 
f(x) = Ci(x) — (cos x — x sin æ — 1)/x°, 


(co) 


Ci(x) = fe cos sds, 


x 


F,= [(4%*B®V, Ro/3ap(1 + aRo)][(061)" te ey )4h 
Ve (KR)?{ (KR)? [(yBo) jo(YEo)/j.(YEo)}? sù 


+ (KR)? — (yRo) jo Ro) RP 
ni = ZZefho,, gt Que, 
B? = 2u(V.— ©), ye = 2u(Vo + ea), 
ka = 2ue,/h? , 


= 


nf RPI aR| + 125447. 


0 


Here k, and k, are the initial and final state relative propagation vectors be- 
tween the target and the hypernuclear mass center; v, and v, are the corres- 
ponding relative speeds; Z, and M, are the charge and mass of hypernuclear 
core; x is the reduced mass of the hyperon and core; & is the binding energy 
of the hyperon and the core; e, is the relative energy after disintegration of 
these hypernuclear segments; and A is the atomic mass number of the target. 
To obtain the total cross-section, expression (5) must be integrated numerically. 

The effect of using Coulomb wave functions rather than plane wave fune- 
tions to describe the relative motion of the hypernucleus and the target is 
contained (in approximation) in the exponential factor in expression (5) (°). 
The plane wave (Born) approximation consists in setting the exponential to 
unity; this yields an upper limit estimate for the cross-section. 


5. — Numerical evaluation; discussion. 


As hypernuclei traverse a detecting medium, they lose energy rapidly by 
ionizing the atoms of the medium (*). A sufficiently accurate expression for 
the slowing down time from an initial energy £, to an energy E (obtained by 

(‘) Elastic collisions may be neglected, since the cross section is small for light 
nuclei, and the energy loss per collision is small for heavy nuclei. 
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5 . . * à 5 ¢ 11 1 
integrating the classical Bohr stopping power formula ( )) is 


(E) = (m/4n2ZNe*)(MI/4m)*(2/M)'¢ , 


6 ee 
à G = Bi[(3) In (4mE,/MI)| — Hi[(3) In (4mH/ MT), 


where N is the atomic density of the medium, m is the electronic mass, M is 
the hypernuclear mass, I is a mean ionization energy for the atoms of the 
medium, taken to be 13.57 eV, and the function contained in G is 


x 


Ei(in x) = Jan TEC o 


0 


After the hypernuclei are degraded sufficiently in energy, they can decay 
only spontaneously, and at all energies considered in the present paper spon- 
taneous disintegration is much more probable than non-spontaneous. Therefore 
the probability for a particular kind of non-spontaneous event relative to 
spontaneous disintegration, taken over the entire history of a hypernucleus 
in the medium in question, is 


(7) N {exp [—t/tlovdt, 


where v is the hypernuclear speed at time t, o is the appropriate cross-section 
and 7 is the mean life-time for spontaneous disintegration (for *He,, taken 
Tokbogi=10=19g). 

These calculations have been carried out for *He, hypernuclei degrading 
from a kinetic energy of 50 MeV to 10 MeV in three liquid media: hydrogen, 
oxygen, and xenon. Characteristics of the slowing down are given in Table I. 


TABLE I. — Characteristics of the slowing down of "He, from 50 MeV to 10 MeV. 


= — E 
Stopping liquid | Density Slowing down Slowing down 
distance time - 1010 
ao 
Hydrogen .05 g/em-3 | 1.18 em 3.57 s 
Oxygen | lo | .146 .436 
Xenon | D | .148 .449 


(11) See, for example, E. FERMI: Nuclear Physics (Chicago, 1950), p. 30. 
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The stripping cross-sections and the electric dissociation cross-sections 


(both with and without the exponential factor) are presented in Table IT. 


TABLE II. — Disintegration cross sections for "He, . 
Stopping Hvd È | A 
liquid NEU Ae xygen Xenon 
Energy Ca Ge Ge ee Og CE Ga Og 
(MeV) (mb) (mb) (mb) (mb) (mb) (mb) (mb) (mb) 
10 0 0 1.14 m0 | 85 ou ~ 0 ~0 
20 014 001 2.9 160,  - 128 47 .28 11 
30 -066 OL). 4) Buss 1.6 142 53 Meal 72 
40 .10 .028 | 3.8 | 2.2 148 59 "152 122 
50 12 .053 9.9 DAT 152 66 14 167 
Cepw: electric dissociation cross section, plane wave approximation. 
de: electric dissociation cross section, Coulomb wave functions. 
GE stripping cross section. 


The use of Coulomb wave functions is seen to be important even at ener- 
gies several times the Coulomb barrier, Z,Ze?/r,. In the case of stripping, this 
is partly due to the fact that as the projectile approaches the target, some 
of its energy is transferred into Coulomb field energy, and it becomes more 
difficult to transfer energy and momentum to the hypernuclear system. For 
both stripping and electric dissociation, Coulomb wave functions describe how 
the Coulomb field reduces the probability for close approaches, while plane 
wave functions do not. 

The relative probabilities for the various disintegration events are listed 
in Table ITI. 


Tape III. — Results of disintegration of 5He,. 


, Electric dissocia- | Electric dissocia- 
Stopping elcome tion probability (*)tion probability (*) Stripping 
liquid Cue sole (plane wave (Coulomb wave | probability (*) 
threshold ||| approximation) functions) 

Hydrogen 17.8 MeV LSS Ose 4.03: 1077 — 
Oxygen 3.84 NES ral Ore ORO AIO 
Xenon 3.02 ASIE RON Ome 1.16-10-4 
(*) Relative to spontaneous disintegration probability. 


Several conclusions may be drawn: a) in the cases considered here, the 
relative disintegration probabilities are small, consistent with the fact that 
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such events have not been observed; 6) stripping appears more likely than 
electric dissociation in those media in which the release of hyperons is most 
probable; ¢) non-spontaneous disintegration appears most probable for mode- 
rate-Z targets. For low-Z targets, the Coulomb field is too weak to produce 
many electric dissociation events, and the target nuclear size is too small to 
produce many stripping events. For high-Z targets, the Coulomb field makes 
close approaches relatively improbable. These conclusions are, of course, valid 
only for those hypernuclei for which the binary model is reasonable (where 
there is not severe distortion of the core or very strong binding of the hyperon 
to the core). 

Although the relative probabilities for non-spontaneous events appear to 
be very small, it seems possible that increased use of short-radiation-length 
detectors and rapid scanning techniques could lead to their observation. The 
appearance of such Coulomb and stripping events would, of course, lend addi- 
tional support to the credibility of the binary model for certain hypernuclei. 


The authors would like to thank Mr. G. A. WILKINS for his aid with the 
calculations. One of us (H.S.V.) would also like to thank Drs. H. PRIMAKOFF 
and M. FRIEDLANDER for valuable comments. 


RIASSUNTO (*) 


Sulla base di un modello binario si ottengono aleune espressioni della probabilità 
della disintegrazione degli ipernuclei in volo. I calcoli numerici per 5He, da 50 MeV 
rallentato in idrogeno, ossigeno ed elio liquidi indicano che tali processi dominano 
maggiormente in mezzi con Z medio. Dei due processi considerati, lo stripping e la. 
dissociazione elettrica, il primo sembra predominare nel campo dei Z medi e grandi. 


(*) Traduzione a cura della Redazione. 


2630 


pre” 


IL NUOVO CIMENTO MO ex VIENE 1° Agosto 1960> 


On the Photonic Decay of the Hyperons. 


L. TENAGLIA 
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(ricevuto il 26 Aprile 1960) 


Summary. — We have evaluated the matrix element for the photonic 
decay of the hyperons: it is given as a function of the matrix element 
for the pion-photoproduction by nucleons. The results we have obtained 
(namely, the branching ratio between photonic and pionic decay of the 
hyperons) agree with the few experimental data known at present. 


Three examples of photonic decay of hyperons (precisely, of the X*-par- 
À L Y; I 


ticle) have been recently found by QUARENI (1) et al, who deduced from an 
approximate analysis of the events that the branching ratio R between, pho- 
tonic and pionic decay of the X*-hyperon has the value & 10°. 

The photonic decay rate has been evaluated in relativistic perturbation 


‘theory by KawAGUCHI and NISBIJIMA, by Iso and KAwAGUCHI and by BEHR- 


ENDS (2). The results obtained by these authors are sensibly different; they 
depend strongly on the nature of the interactions taken into account and do 
not give the experimental branching ratio R. Moreover, such results are not 
directly connected with the isotopic spin I and orbital angular momentum 
I states of the final system. An evaluation which should be more reliable, 
from a physical point of view, can be performed starting from the matrix ele- 
ments of pion photoproduction which we know from the CHEW-GOLDBERGER- 


(1) G. QUARENI, A. QUARENI-VIGNUDELLI, G. Dascora and G. Mora: Nuovo 
Cimento, 5, 1179 (1959). 

(2) M. Kawaeucui and K. NISHIJIMA: Progr. Theor. Phys., 15, 182 (1956); C. Tso 
and M. Kawacucui: Progr. Theor. Phys., 16, 177 (1956); R. E. BEHRENDS: Phys.. 
Rev., 111, 1691 (1958). 
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Low and NamBu analysis (*). Such a procedure, used in the following evalua- 
tion, is not free from some objections, since it is necessary to extrapolate 
the CGLN matrix elements, which are in good agreement with the experimental 
data at low energies, to the higher energy values which appear in the inter- 
mediate states. In addition, the matrix elements for the pionic decay of the 
hyperons are not known for any values of J and J, a fact which turns out to 
be very important in the present case: the only (very scarce) information rela- 
tive to such a process are also restricted to the low energy region. 

For such reasons, in the following we shall be compelled to make some 
rough approximations. However, even a calculation of this kind leads to a 
correct order of magnitude for the branching ratio À and, moreover, it allows 
to single out the angular momentum J and isotopic spin I states for which we 
can presume that the decay ratios dominate. 

The matrix element relative to the photonic decay of a hyperon, Y —.N ++, 
is given (4) by: 


ik 


a, Gail (27)! SP +R O4 


(1) MS = ir fay exp [—iky]<N | T(jy(y) Lx(0)) |0> wy , 


Ti) 1 (0)) = Livy), LOI (y) + Le) VY) , 


where P, k, Q stand for the energy-momentum of the nucleon, the photon and 
of the hyperon, respectively: T indicates the Wick ordered product; j,, I, are 
the current densities, sources of the photonic and of the hyperonic fields. 

With respect to that relative to the process Y > V+n—>N-+y, all the 
other matrix elements are assumed to be negligible. As a matter of fact, it 
is difficult to estimate, in the present case, the relative orders of magnitude 
of the different matrix elements of which .#, is the sum. With the only aim 
to have a tentative comparison, we have performed a relativistic perturbative 
evaluation of the matrix element for the processes 


(2) ae Nn Dae, 


(2°) DO BT > (N+P +N) So (r+) op 


(3) G. F. CHEw, M. L. GoLDBERGER, F. E. Low and Y. NAMBU: Phys. Rev., 106, 
1345 (1958); in the following indicated by CGLN. 

(*) v indicates the polarization of the emitted photon, the electromagnetic potential 
of which is A=e(2k))-? exp [ik#]. We use natural units hem 
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and we have verified that, at least in this approximation, the matrix element 
for the process (2) is much bigger (by a factor > 10) than the matrix ele- 
ment for the process (2’). 

Assuming to be sufficient to take into account only the process (2), one gets 


(3) Mi 2) 0PLk_P)d_ (Pa P, — hy) <N | 4, |n><n| Ly | Oday . 


The matrix element <N|j,|n> is known, at least for small values of 
|(P+k—P,)|, from the CGLN analysis. We can also write (5) 


j M \3 Pai 
(4) <n |Iy|0> uy = (=| > 1e [as + Woz 29) HA 


where y stands for a Pauli spinor: a,,, b,, represent the amplitudes of the virtual 
states (in the center of mass system) relative to the eigenstates of orbital 
angular momentum /=0 and /=1 (in the present case only the states with 
angular momentum J— + are possible). 

Introducing the density of the intermediate states, as a function of the 
pion energy-momentum, the matrix element of eq. (4) becomes 


i at A q? dq 
R OY eee See | 10. 
©) Mts em & He ma | 2/3 a a fi 2 


AYE Ak eo 
fe Leen) 7 (a. 427 A 4% » 


Do (20) d 


where f(P, q, k) is the photoproduction matrix element. The following symbols 
are used: 


o (gi 41), E =|[{M%||PLk=qg|}, kb=lb+DP—_; 


{ is the solid angle in the q space; 7 represents the Cauchy principal value. 
In the CGLN analysis, the nucleon recoil is approximately taken into ac- 
count in the photoproduction: in this approximation one gets 


P oq do 
Dore *(k, q)- 
qo ans + ere (&, @) 


M nee CoP 
‘p aot (ta 1D», Nye, 
0 a 


a ae 


with k= (1+ 2/M). 


(5) R. P. Dazrrz: Phys. Rev., 105, 605 (1958). 
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a,, and b,, are functions of © and of po = (1-+q2/M?)?: as a matter of fact, 
it is very difficult to determine with a good approximation the dependence of 
a,, and b,, on ©. This is important in connection to the renormalization 
corrections for the weak interaction coupling constants. In fact, .#, can be 
evaluated by means of dispersion relations, having taken into account only the 
matrix element relative to the process (2). A good knowledge of the quantities 
which appear in the higher order interactions would allow to evaluate the 
corrections to the Fermi interaction constants, both for the appearance of 
form factors and for the presence of induced terms (e.g. a PS interaction 
can arise from the V— A theory, as it was proved by GOLDBERGER and 
TREIMAN) (6). 

All these corrections are functions of the released energy-momentum, 
(P,—QY. 

The assumption that a,, and b,, are constants is equivalent, in the eval- 
uation of the branching ratio between the photonic and the pionic decay, to 
the hypothesis that such corrections have the same value both at the low 
energies for the pionic decay and at the high energies of the virtual states which 
occur in the matrix element (6). A rather approximate analysis, carried out 
by means of the Goldberger and Treiman (*) procedure, shows that a, and 
b,, can strongly decrease at high values of the energy-momentum released, 
(P,—Q). It can be also observed that the CGLN photoproduction matrix 
element cannot be corrected at high values of |P+k—P,|. In fact, if one 
evaluates the photoproduction cross-section with the preceding function f(P, q, k) 
—including only J= 4 states—one obtains a value which is increasing with 
the energy, in disagreement (*) with the experimental data. For this reason, 
in the evaluation of the matrix element (2), which involves divergent integrals, 
one must assume that the product f*(a,,+b,,) tends to zero for w > «,, where 
©, is a convenient maximum energy. This can be accomplished by introducing 
a cut-off function and, with such approximation, one can take into account 
the decrease of the amplitudes a, and b,, at high energies. 


(9) M. L. GotpBeRGER and S. B. TREIMAN: Phys. Rev., 110, 1178 (1958); 111. 
354 (1958). 

(*) In fact, from the values reported by CGLN one deduces that the photopro- 
duction cross section, in the angular momentum state J=4, is a mixture of terms. 
which are decreasing or constants with increasing q, and of terms of the type 


qÿ q 
LIA Beal’: 
kw? w?(w + W,) : 


which increase with q. They derive from those matrix elements which make divergent, 
in the present calculation, the integral over the intermediate states 


2634 


ON THE PHOTONIC DECAY OF THE HYPERONS 427 
From the matvix element (6), keeping a,, and b,, constants, one gets: 
lord 
(7) Mer" >) (4,,0; 4+ bf) - 


"A 


The term N*(m), which does not give an appreciable contribution, has been 
omitted = the CGLN matrix elements #*, FT and ¥°, the sum of which 
is f(q, k), aside from a factor 27. One ua in this way 


al F(q) do | 
= kb) + — 
7 — O° er ( ) EE MJ’ 


Go 1p sed | 
Ho ZA (K°) 4 ni) © — ky + 20/ MJ ’ 


| 
(8) 
| 


9 
(9) F(q) = fegw/[2w(1 + q2/M?)}} | E, i(0, + 63) Ps + fi Jt + 
i i - Jx | 
I Az = 9 SIAT | Mia = Si 1) — - ® 3 
F SRE ES (20, + 0:)F,— q?/3.Mc | J- + P 2/6Mo Sar ee J | 


In the function #(q), only the terms proportional to the magnetic dipole 
transitions (that is, the terms proportional to g?/3Mw) diverges linearly, but, 
owing to the factor (3W)-1, their contribution is small for © <6m_. The other 
factors are convergent or they diverge logarithmically: the value of «, has 
not practical influence on them. 

As far as G(q) is concerned, it is very small: assuming, for its evaluation, 
@) = 6m, we have verified that the contribution of @(q) is of the order of 10-3 
with respect to the contribution of i(q). We can also say that the choice 
of ‘w,, in the present case, is not very critical; a variation Ao ~m, does not 
alter sensibly the results. 

For the sake of convenience we put 


(10) a, =—a°e3 (A; — A7 + (BY + By )]d* + |A; + 4; J 

+ i(B, + B;)]} JT +i As}, 
(11) B,=—io kxe3 "(0 +iC)(Jt + 7), 
where the indices I in A+, BŸ,... refer to the I isotopic spin states. The 
matrix element for the decay of the 2*-hyperon is given by 


(12) M s+ wy, — Ms op y - M 


Z+,rta,py * 


DTT 


The two processes contribute coherently to the decay: one can then write 


(13) 3M s+ py = — 060424; + Aî + i(2B; + BY + 349)] — 
— Da, [AS + A; + iB} — B;)]}}— io-kx e(C + 10°)(3d,, + 27%,,) è 
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A not very different expression is obtained for the matrix element of the 
>° and A-decays. In the evaluation of the decay probabilities there appear 
four distinct amplitudes: notoriously, they cannot be determined from the mean 
lives only. 

By putting 


As Dh Dog as Mo 2 2 
RE TR At=1i4+@4+yi+ y 


Tia, ae ; Oey 


we get, for the evaluated branching ratios, 


[ Ry: = 0.30 -10-2(1 — 0.52x + 0.06847) A , 
(14) | Ryo = 0.29 -10-2(L1 — 0.517 + 0.0654?) A, 
| R 


A = 0.685-10 2(1 + 0.854 -- 0.21 47) 4 


The sensible difference between R,,, Ry, and EF, is due essentially to the 
much more limited value of the available energy, in the last case. 

One can at last observe that many examples of pionic decay of the A have 
been already found, but no photonic decays have been observed up to now: 
on the contrary, the only photonic decay detected is the X*-decay. Besides 
the experimental difficulties, it is obvious how a value «~1 can modify the 
branching ratio by a factor = 7~8 respect to the value at x = — 1. From 
this point of view, even if the relation (14) could be affected by an error of 
about 20%, it is always evident the clear difference of the behaviour of R,., 
E. and E, with the variation of the sign of +, for |v|~1. 


* OK OK 


I am very glad to thank Prof. M. Cini and Prof. R. GATTO who kindly 
revised the draught of this elvaluation. 


RAS SEN @ 


Si calcola il rapporto fra la probabilità di decay fotonico e pionico degli iperoni. 
L'elemento di matrice pel decadimento fotonico è espresso in funzione di quello di 
fotoproduzione di pioni da nucleoni. I risultati ottenuti sono in buon accordo con i 
pochi dati sperimentali attualmente noti. 
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Triggered Spark Counter Arrays of Large Area (Square meters) 
for Experiments on Very High Energy Cosmic Ray Particles. 


N. B. Mistry (*), G. T. Murray, P. V. Ramana MURTHY, 
and B. V. SREEKANTAN 


Tata Institute of Fundamental Research - Bombay 


(ricevuto il 23 Giugno 1960) 


Summary. — The constructional and operational features of triggered 
spark counters of large area are described. These counters are being 
used, in conjunction with nuclear emulsions, in current experiments in 
this laboratory on high energy nuclear interacting particles which arrive 
simultaneously over large areas, sometimes associated and at other times 
unassociated with air showers. For these experiments the spark counters 
and emulsions are being used in arrays covering an area of a few square 
meters. 


1. — Introduction. 


Experiments are in progress in this laboratory to study in detail the fre- 
quency, multiplicity, energy distribution and lateral separation of high energy 
nuclear interacting p rticles associated and unassociated with extensive air 
showers. The frequency of these events is so low that it is necessary to have 
large detecting areas, of the order of several square meters, in order to record 
a sufficient number of events in a reasonable time of operation. Because of 
this requirement of a large area it is not practicable to use multiplate cloud 
chambers for a comprehensive investigation. In looking for an ideal detector 
that can give information on all the aspects mentioned above and which can 
be built economically to cover large areas, it was realized that this could be 
achieved with a suitable combination of spark counters and photographic 


(*) Now at Columbia University, U.S.A. 
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emulsions (*). For this it would be necessary to construct large area spark 
counters in which a visible spark would be produced at the position at which 
a beam of extremely collimated charged particles passed through it. The 
set-up that was envisaged was the following. The nuclear interacting particles 
would interact in an absorber and produce a narrow jet of charged and neutral 
pions; the y-rays resulting from the decay of the neutral pions would, on 
passing through a lead converter below, give rise to a collimated beam of 
electrons; the beam of electrons and mesons then passes through a thin hori- 
zontal layer of nuclear emulsions followed by a spark counter immediately 
below. If the spark counter produces a visible spark at the position at which 
the beam of particles passed through it, then the co-ordinates of the beam in 
the emulsions would be known; the number of electrons in the beam can be 
counted in the emulsion and the approximate size of the interaction determined; 
a system of two spark counters one below the other, would enable a deter- 
mination of the direction of the jet and would also fix the co-ordinates in the 
emulsion more accurately; the occurrence of simultanecus parallel sparks would 
reveal the presence of multiple nuclear active particles; the system could be 
put in association with an extensive air shower (EAS) array and the size and 
core position of the associated shower determined; if the core of the EAS 
passed through the area of the spark counters, then a very precise location 
of the position of the core would be possible; it is obvious that this system 
would work only at very high energies where the particles would be produced 
with the necessary collimation. These very attractive features prompted a 
detailed experimental investigation on large area spark counters. 

In the first instance a spark counter comprising a metal plate on one 
side and a set of parallel wires on the other was constructed with a sensitive 
area 50 em x 27 em and was operated with a continuous (D.C.) voltage. This 
counter had an excellent response to «-particles but there were two serious 
disadvantages. Firstly, there were quite a number of spurious spark break- 
downs and secondly, the ozone produced as a result of the large corona dis- 
charge attacked all parts of the spark counter. To get over these disadvan- 
tages, an attempt was made to operate this type of counter on pulsed voltages. 
A polonium «-source was deposited on the plate of the spark counter and a 
mica window proportional counter fixed opposite to this source. The voltage 
was pulsed on the spark counter whenever on a-particle entered the propor- 
tional counter after passing through the gap of the spark counter. It was 


(*) This system is, in principle, similar to the one reported by GRIGOROV et al. (1) 
who have used photographic emulsions in conjunction with ionisation chambers for 
the investigation of high energy interactions. 

(1) S. I. BRIKKER, N. L. Grigorov, M. A. KONDRATEVA, A. V. PODGURSKAJA, 
A. I. SAVELEVA, V. A. SOBINJAKOV and V. Y. SESTOPEROV: Suppl. Nuovo Cimento, 
2, 733 (1958). 
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found that with pulsed voltages the spark counter did not respond to «-par- 
ticles. All possible parameters like spacing, wire diameter and clearing field, 
as also the amplitude, duration and rise time of the pulsed voltage etc., were 
tried but the response was still negative. 

After considerable experimentation a design was evolved which met with 
all the requirements stated above and which was found most suitable from 
the point of view of construction and operation and also from the point of view 
of convenience in photographing the sparks. In this paper, are given the 
details regarding the constructional and operational features of such spark 
counters developed in this laboratory. 

References to earlier work on spark counters are listed below (2). 


2. — Constructional features. 


Specially chosen cylindrical brass welding rods of good surface and linearity, 
each 60 em long and 3 mm in diameter, are placed parallel to each other on 
a plastic frame 60 cmx 30 cm with a spacing of 3 mm between the rods; the 
rods are permanently stuck to the frame with araldite. Aluminium angles 
are fixed to the frame to keep it straight. A highly polished stainless steel 
plate of thickness 3 mm and of size 50 cemx 30 cm is treated with dilute 
hydrochloric acid to remove surface contamination and then washed and dried. 
This plate is fixed to the frame by means of 2 BA screws and nuts, with a 
spacing of about 3mm between the rods and the plate; the spacing can be 
varied by changing the washers. A transparent plastic sheet is fixed over the 
frame to prevent dust from settling down in the gap of the spark counter. 
All the rods are connected together by wires. A photograph of the spark 
counter is given in Plate I. 


8. — Operational features. 


8 1. Pulsing. — The stainless steel plate is connected to earth and a nega- 
tive voltage pulse whose amplitude could be varied from 2 to 5 kV, is applied 
to the rods whenever required. The rise time of the pulse is kept less than 


(2) W. Y. CHANG and S. RosENBLUM: Phys. Rev., 67, 22 (1945); A. Ryvs: Helv. 
Phys. Acta, 22, 1 (1949); R. D. Connor: Proc. Phys. Soc., B 64, 30 (1951); N. K. Sana 
and NARENDRA NATH: Nucleonics, 15, 94 (June, 1957); M. J. SweTNICK and N. G. 
Anton: Nucleonics, 15, 93 (June, 1957); R. W. Pipp and L. MADANSKY: Phys. Rev., 
75, 1175 (1949); J. W. KEurFEL: Rev. Sci. Instr., 20, 202 (1949); E. ROBINSON: 
A 66, 73 (1953); F. BeLLa and C. FRANZINETTI: Nuovo Cimento, Proc. Phys. Soc., 
10, 1461 (1953); T. E. CransHaw and J. F. pp BEER: Nuovo Cimento, 5, 1107 (1957). 
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0.05 us to have a high efficiency of response, and the decay time kept of the 
order of one microsecond to minimize spurious breakdowns. The pulsing is 
done using a thyratron circuit similar to the one described by CRANSHAW and 
DE BEER (2). A low clearing field (~ 20 V) is maintained continuously be- 
tween the rods and the plate for removing old ions; this helps in reducing 
spurious breakdowns. The counter is operated in complete darkness to mini- 
mize breakdowns due to photoelectric effects. 

With a higher voltage of (10+14) kV, obtained from a pulse-transformer, 
it is found that the same spark counter can be operated with a spacing ~ 6 mm, 
which is almost double the previous value of 3mm. It is of importance to 
operate the spark counter with a larger spacing, since the requirement of uni- 
formity of spacing is then less stringent. This makes it possible to build indi- 
vidual spark counters of very large area about half to one square meter. 


3°2. Photography. — The chief advantage of this type of spark counter over 
parallel plate spark counters is the convenience and simplicity in the photo- 
graphy of the sparks. In parallel plate spark counters, the spark has to be 
necessarily photographed through the gap; for this the alignment of the cameras 
has to be very accurate, and also a very large depth of focus is involved. The 
latter requirement can be met by using pinhole apertures, but for this the 
intensity of the spark has to be very great. It is not desirable to increase the 
intensity of the spark, since the surface will be damaged and spurious break- 
downs will be greater. For these reasons spark photography in parallel plate 


Water tank J eS —-—-=- 


Geiger counters — FL 
spark counter 1 = aay 


nuclear emulsions 


Spark counter 2 


as Camera 


Liquid scintillator 


LLL ye Lead 2" 
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Fig. 1. — Experimental arrangement to test the feasibility of employing nuclear emul- 

sions in conjunction with double spark counter array for the study of high energy jets. 

The coincidence G, SG, selects nuclear interactions produced mostly in the water 

tank. The lead plates above the spark counter are meant for the rapid development 

of the cascades. The direction and co-ordinates of the jets are given by the position 
of the sparks in the two spark counters. 


Gelger counters 
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counters of large area poses very difficult technical problems. Moreover, to 
get the spatial co-ordinates accurately, stereoscopic photographs have to be 
taken, or else the sparks have to be photographed from two perpendicular 
directions. In the present spark counters the photography is done with just 
one camera and in a very simple way using a system of mirrors suitably inclined, 
as shown in Fig. 1. 

The sparks can be located to an accuracy of (3 mmx3 mm) which is suf- 
ficient for the experiments envisated. A typical photograph of a spark is shown 
in Plate II. 

This type of photography using mirrors imposes one limitation, namely, 
if two spark counters are used one below the other, then the spacing between. 
them cannot be reduced below a limit. However, this is not a serious dis- 
advantage since a certain minimum spacing between the spark counters is 
essential for an accurate determination of the direction of the jet and also for 
establishing the parallelism of the jets when there are multiple nuclear inter- 
acting particles incident. 


3.3. Efficiency. 


i) Response to singly charged particles at minimum ionization: To de- 
termine the operational characteristics for minimum ionizing singly charged 
particles, the spark counter was placed between two geiger counters registering 
the passage of cosmic ray 


particles. The voltage on 100 
the spark counter was pulsed 80 

: | ‘+ 
whenever there was a coin- Sa evi 
cidence between the Geiger  & 9° Cee i 
counters. It was observed  x40} 

. Cc 

that sparking took place 8 20! Cleanngficia 5 y 
almost always (99%) in the è no delay 
region covered by the Geiger 3500 3600 3700 3800 3900 4000 4100 4200 
counter telescope and very Pulse amplitude (V) 
rarely (1%) outside. The Fig. 2. — Variation of the efficiency of the spark 
percentage of spurious break- counter for single cosmic ray charged particles, as a 


downs as obtained from this function of the pulsed voltage. 


observation and by random, 

triggering was less than 1%, at normal operating voltages. The efficiency of 
the spark counter, defined as the ratio of the number of spark breakdowns 
to the number of double coincidences, was determined as a function of the 
applied voltage, the clearing field voltage and the delay in the application 
of the pulse. The results are shown in Figs. 2 and 3 for a particular spark 
counter. It was found that by careful adjustment of spacing and applied vol- 
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tage the efficiency for single particles can be pushed up to almost 90%. The 
counter was found to be equally sensitive whether the particles passed through 
the rods or in the gaps between the 


ae sus delay : rods. 
er ie = no delay ae Response to «-particles: ga 
2 20 Pulsed voltage 4100 V the same counter, when the operating 
Ò £ voltage was adjusted for an efficiency of 
è 50 100 about 5% for singly charged minimum 
x GG Uol) ionizing particles, the efficiency of re- 
Fig. 3. — Inefficiency as a function of sponse to «-particles was almost 100%. 
clearing field, and delay in the appli- = This was checked by depositing a thin 
IU ae eg RI AI film of a polonium «-source on the 


plate and pulsing the voltage by a pro- 
portional counter as before (see Section 1). This experiment shows that by a 
suitable adjustment of the efficiency, it is possible to make the spark counter 
respond preferentially to collimated jets that to some extent simulate «-par- 
ticles in ionization. 


4. — Operation with other particle detectors. 


To determine the optimum conditions of operation of the spark counter 
for the detection and location of high energy jets produced by cosmic ray par- 
ticles, the spark counter was operated in conjunction with a multiplate cloud 
chamber and with nuclear emulsions. The experiments were carried out at 
Ootcamund (2.2 km above mean sea level). 


41. Operation with the multiplate cloud chamber. — The spark counter was 
mounted above the multiplate cloud chamber. Above the spark counter was 
kept a lead block of 10 em thickness in which nuclear interactions were pro- 
duced and which also served as a medium in which the nuclear and electro- 
magnetic cascades developed. The spark counter was set to operate at an 
efficiency of ~ 0.2% for singly charged minimum ionizing particles. The cloud 
chamber and the spark counter were triggered for penetrating showers, with 
a scintillator-geiger counter selection system; since the selection system was 
not very rigorous, it responded to penetrating showers and air showers. The 
cloud chamber photographs provided visual evidence regarding the type of 
events which actuated sparks in the spark counter. It was found that the spark 
counter responded to high energy penetrating showers as expected; it was 
also found to respond sometimes to air showers. The response to air showers 
could not be explained in terms of the efficiency for single particles, since the 
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Prate III. — Cloud chamber photograph of a high energy penetrating shower pro- 

duced in the lead block above the chamber and developing further in the lead plates 

inside. A spark was seen in the spark counter in the region through which the axis 
of the shower passed the spark counter (see the superposed illustration). 
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number of such particles passing through the spark counter was small, as seen 
in the corresponding cloud chamber photographs. The response is due most 
probably to heavily ionizing fragments and collimated showers of the type 
commonly seen emerging from the plates in multiplate cloud chambers. To 
ensure that the spark counter responds. only tq high energy penetrating showers, 
and to reduce considerably its response to air showers in which the high energy 
nuclear active particles do not pass through the spark counter area, a system 
of two spark counters one below the other can be employed. Correlated sparks 
in the two spark counters can only be due to high energy penetrating showers. 
The double spark counter system also has the additional advantage of giving 
the direction of the jet and of establishing parallelism in the case of simul- 
taneous jets. 


4°2. Operation with nuclear emulsions. — A system of two spark counters 
Was operated at Ootacamund, in conjunction with a penetrating shower se- 
lection system as shown in Fig. 1. A layer of Ilford G-5 200 um thick nuclear 
emulsions was placed above the top spark counter. The spark counter voltage 
was pulsed whenever there was a coincidence of the type G, S G,. The scin- 
tillator was biased such that the G, S G, coincidence rate was about 2/hr. These 
coincidences are mostly due to penetrating showers with energies > 20 GeV. 
Air showers in which the nuclear active particles do not pass directly through 
the array will not be admitted by this triggering system. The results are given 
in the following Table I. 


Tappe I, 
| Efficiency of spark Sparks 
Ti | sounter for single - n = 
Time of | Pulsed ì particles = Number 
oes | voltage | of In top | In bottom | D 
ation | 7 c i i i, 
(h) | (V) Top spark | Bottom spark negre | Kite ten | counters 
| counter counter only only 
28 | 4070 4.7%, 8.5%, 41 10 S| Se 
42 | 3850 4.1% 4.0% 69 12 5 | 10 
22 Wis 3.740 2 |, 1.0259 1.5% (O 7 3 | 5 


It follows from the table that when the efficiency of each spark counter 
is adjusted to about 5% for singly charged minimum ionizing particles, the 
efficiency of the system for penetrating showers, 7.e. that one or both spark 
counters are triggered, is more than 50%. As already mentioned, the triggering 
rate of 2/hr corresponds to penetrating showers of energy about 20 GeV, and 
at this energy the collimation of the particles is poor and the number of electrons 
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produced by the development of y-rays is also small. An efficiency of 100% 
cannot therefore be expected. The double sparks, which correspond to very 
high energy penetrating showers constitute about 10% of the number of trig- 
gering pulses. If the latter are due to showers of energy > 20 GeV, as pre- 
viously stated, the double sparks should be attributed to showers of energy 
>100 GeV. In all cases where double sparks were recorded a search was made 
in the emulsions at the places (5 mm x5 mm) indicated by the spark counters 
and in the direction defined by the line joining the sparks, for multiple parallel 
tracks due to the collimated nuclear electro-magnetic cascades. In some cases, 
(5-15) tracks confined within a radius of ~ 50 um were observed. One can 
expect 5 tracks confined within a circle of 50 ym under about 5 radiation lengths 
of lead when the energy of the y-rays is more than 50 GeV. Such y-rays will 
in general arise from nuclear interactions involving energies of about 10” eV. 
It is difficult to detect in the emulsions events in which the number of tracks 
is <5 in a circle of radius 50 um. Thus the lower limit to the energy of the 
nuclear interacting particles that can be detected with the combined system 
of spark counters and nuclear emulsions is about 10! eV. Because of this, 
in all recorded cases of double sparks the corresponding events have not been 
detected in the emulsions. Nuclear emulsions of 50 um thickness would reduce 
the cost considerably if they could be used. On trials, however, they were 
found to be unsuitable; because of their small thickness it was difficult to 
observe the tracks and establish definitely their parallelism. It should be 
feasible to use a thickness < 200 um, possibly — 100 um. It may also be 
emphasized that only those squares of emulsions need be taken out for exam- 
ination through which jets have passed as indicated by the double sparks. 

An important experimental observation is that even though all the rods 
were connected together electrically, simultaneous multiple sparks were re- 
corded, sometimes quite close to each other (few cm). Therefore, it is not 
necessary to use many pulsing circuits. 


5. — Conclusions. 


It is feasible to construct triggered spark counters with areas of the order 
of half to one square meter, which respond preferentially to high energy jets 
(> a few hundred GeV) produced by cosmic ray particles. Double spark 
counter systems can be used to give the position and direction of high energy 
jets and will be particularly useful in the study of simultaneous high energy 
nuclear interacting particles (*4) associated and unassociated with air showers. 


(3) N. Thy GRIGOROV, V. Y. SHESTOPEROV, V. A. SOBINYAKOV and A. V. Popcur 
SKAYA: Zurn. Eksp. Teor. Fiz., 33, 1099 (1957). 

(4) S. NARANAN, R. RAGHAVAN, P. V. RAMANAMURTY, B. V. SREEKANTAN and 
A. SUBRAMANIAN: Nuovo Cimento, 16, 401 (1960). 
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Nuclear emulsions can be used in conjunction with such double spark counter 
system to obtain estimates of the energies of the nuclear interacting particles. 
Areas of the order of several square meters can be covered with an array of 
such large area spark counters and nuclear emulsions. 


‘ 


We have great pleasure in expressing our gratitude to Professor B. PETERS 
for his interest and encouragement in this investigation. We are thankful to 
Professor M. G. K. MENON for helpful discussions. We are indebted to Dr. M. 8. 
Swamy for help in the work connected with the use of nuclear emulsions. 
Our thanks are due to Messrs. F. GONSALVES and A. R. APTE who have helped 
us in preparing the spark counters and in operating the cloud chamber. 


RIASSUNTO (*) 


Si descrivono le caratteristiche costruttive di contatori a scintillazione di grande 
superficie. Questi contatori vengono usati, in unione con emulsioni nucleari, per espe- 
rimenti correnti in questo laboratorio su particelle ad alta energia da interazioni che 
arrivano simultaneamente su larghe aree, alcune volte in associazione con sciami del- 
l’aria ed altre volte senza associazione con questi ultimi. Per questi esperimenti i con- 
tatori a scintillazione e le emulsioni sono usati in disposizioni che coprono un’area di 
alcuni metri quadrati. 


(*) Traduzione a cura della Redazione. 
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ponsabilita scientifica degli scritti inseriti in questa rubrica è completamente lasciata 


(La res ques 
dalla Direzione del periodico ai singoli autori) 


Effect of Lepton Non-Conservation on ~-Decay. 


IP, IR, JAUNE 


Saha Institute of Nuclear Physics - University of Caleutta 


(ricevuto il 20 Giugno 1960) 


It has been pointed out by PauLI(!) that first-order processes throw no light 
whatever on the question of lepton conservation, a principle which is more or less 
taken for granted in almost all discussions of leptonic processes. The experiments 
on second-order processes, like neutrino capture and double B-decay, which could 
decide the question, present extreme experimental difficulties and provide only a 
very rough limit of the amount of possible lepton non-conservation. Consequently, 
it is desirable to examine other processes in which the law of lepton conservation 
can be tested. It is the purpose of this note to show that the branching ratio of 
m-e to ru decays is a sensitive test for the presence of a particular type of lepton 
non-conserving interaction. 

The branching ratio R of x-e to m-u decays is found to be (?) 


Wirt > e+ +) : 
=o Omer 


ae aa W (xt = ut +y) 


(1) 


with a quoted error of less than 20%. We know that pions interact strongly with 
nucleons and also that nucleons possess weak interactions with muons and with elec- 
trons. To understand the pion decay, therefore, it is not necessary to invoke any 
direct interaction between pions and leptons; in fact, according to the general prin- 
ciples of quantum mechanics, we should expect the pion to decay through a two- 
stage process whereby the pion first dissociates through its strong interaction virtually 
into a nucleon-antinucleon pair which trasforms into a lepton and a neutrino through 
the weak interaction. 

Since there is no experimental evidence which requires us to postulate a direct 
decay coupling between pions and leptons, we shall assume that there is no such 
direct interaction. Assuming that the pion-nucleon interaction conserves parity 
and that the pion field is pseudoscalar, it follows from general invariance arguments 
that the decay proceeds only through the P and A forms of the Fermi interaction. 


(4) W. PauLI: Nuovo Cimento, 6, 205 (1957). 
(@) H. ANDERSON, T. Fugu, R. H. MILLER and L. Tau: Phys. Rev. Lett., 2, 53 (1959). 
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Unfortunately, we cannot in general calculate the decay rates of either of the two 
modes because of the presence of the strongly interacting particles in the intermediate 
state (*). However, for the case of either pure P or pure A interaction, the branching 
ratio becomes independent of the strong interactions because the effect of the strong 
interactions may then be represented by a constant factor which is the same for 
both modes. Assuming lepton conservation and equal coupling constants for 8-decay 
and muon absorption, the calculated ratio for the first case is far in excess of that 
observed, (1), but for the A case, the ratio is (*) 


2 2/02 
m, Imm 

Re oe aL LO ts 
m,| 1 —m,/m,, 


which is in surprisingly good agreement with (1). Further, inclusion of electrodynamie 
radiative corrections (°) reduces R by about 11% for an energy resolution of 1 MeV. 

It is of course possible that the agreement between (1) and (2) is purely acci- 
dental but we shall not adopt that view, particularly since we note that 8-decay 
certainly contains A interactions (°) and it appears that the Hamiltonian for muon- 
absorption does too (7); the strengths are known to be comparable. Also, there is 
no positive evidence for P interaction in either case. From the agreement of (1) 
and (2), we may exclude the possibility of any P interaction, or more accurately 
any P interaction with a coupling constant of the order of, or exceeding, 10-* times 
the A coupling constant (8). Strictly speaking, we cannot calculate the ratio À for any 
mixture of interactions but we quote this figure because a perturbation-theory esti- 
mate with a cut-off at the nucleon mass indicates that, for equal coupling constants, 
the amplitude for decay due to P interaction exceeds that due to A by a factor of 
the order of 10%. To understand how the introduction of lepton-non-conserving inte- 
ractions affects R, we recall briefly why the old lepton-conserving theory gave the 
small value (2) instead of the much larger number which we should expect from 
phase-space arguments. The expression in square brackets in (2), which is greater 
than unity, is essentially this factor. 

It is a feature of the A interaction (and the opposite for P) that if we describe 
the neutrino, which we assume to be massless (°), in terms of handed states, negatous 
are emitted preferentially with the same handedness as the neutrino (which we 
define by n->p+e-+v) and positons the opposite; the probability of the favoured 
spin-state being 3(1+v/c) as compared to 3(1— v/c) for the other. Therefore, since 
the decay x*+->e*-+y produces extremely relativistic electrons (electron energy 


TC 


— 70 MeV) the probability of emission of positons with handedness opposite to that 


(?) Nevertheless, M. L. GoLDBERGER and S. B. TREIMAN have succeeded in obtaining a remarkably 
good estimate of the absolute decay rate on the-basis of certain plausible assumptions: Phys. LRewv., 
110, 1178 (1958) and 111, 354 (1958). 

(4) M. RuDERMAN and R. J. FINKELSTEIN: Phys. Rev., 76, 1458 (1949). 

(5) S. BERMAN: Phys. Rev. Lett., 1, 468 (1958). 

(6) C. S. Wu: Rev. Mod. Phys. 31, 783 (1959). 

(7) V. TELEGDI: see J. POLKINGHORNE: Phys. Today, 12, 12, 36 (1959). 

(8) It is possible for a small P admixture of this order to reduce the ratio R through destructive 
interference, and this was pointed out by HUANG and Low at a time when it appeared that m-e decay 
was even more rare than (1). We do not consider this pathological possibility. 

(°) This is only to simplify the discussion, and is not essential to the argument. The assumption 
of a finite rest-mass (compatible with other experimental evidence) would not affect the principal 


conclusion. 
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of the neutrino, is almost unity. But consideration of angular momentum conser- 
vation in the rest system of the pion (J=0), in which the positon and the neutrino 
go in opposite directions, requires both particles to have the same handedness ; 
since the probability of emission of such positons is so small, the decay occurs very 
seldom. In the case of the muon-decay mode, the muon is quite non-relativistic 
(muon energy = 4 MeV) and the probability of either handed state is almost equal. 
There is therefore no inhibition due to the angular momentum condition and the 
decay proceeds at a «normal» rate. Notice that no assumption was made about 
a particular handedness for the neutrino, thus the argument does not depend on parity 
non-conservation or on the two-component theory of the neutrino. 

But now we see that the argument rests on the assumption that positive pion 
decay involves the emission of a neutrino and not of an antineutrino, for which 
the electronic decay mode would again occur at a «normal» rate. Thus a slight 
relaxation of the law of lepton conservation would have a large effect on the 
branching ratio R. The general theory of Fermi interactions including possible lep- 
ton non-conservation has been given by PAULI (!); the B-decay Hamiltonian, for in- 
stance, takes the form 


(3) Ho => ali) [fg T + gin Tv + (vy OU Ls + fers lave} + hc. , 


where the first term in the square brackets represents the conventional Hamiltonian, 
including parity-non-conserving terms. Taking the pseudoscalar coupling constants 
as zero (‘°), the branching ratio calculated from interactions of the form (3) for 
muon-absorption and B-decay is also independent of the strong couplings: 


Bia 
(4) pp as Mame DI 
0% + m2 /m2 DI’ 
(3 PTT SI 
where 


= ge |gijî, - Di=|L4 |. 


If we take the muon and electron coupling constants as equal and consider the rela- 
tive strength e = D?/0? of the lepton-non-conserving interaction as small, (4) can 
be rewritten as 


RONN N= — = 745-104. 


Thus we see that an admixture of 10-5 of lepton-non-conserving A interaction would 
almost double the branching ratio R, and we may conclude from the agreement be- 
tween À, and R.., that any such admixture must be considerably smaller or non-existent. 
The selection rule which prevents S, V and T interactions from contributing to pion 


(°) The argument given earlier for the smallness of the P coupling constant applies equally to 
lepton-conserving and non-conserving interactions, However, in the lepton—non-cunserving case 
a small P admixture would manifest itself as much in enhancing the total decay rate as in char 
the branching ratio R. i ù 
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decay also prevents us from drawing any conclusions about possible lepton- 
conserving terms of those forms. Similar arguments apply to the Keo/Kyy branching 
ratio, although the greater experimental uncertainty in this case permits only the 
weaker restriction e < 1078. ; 

It should be remarked that another requirement for the agreement of (1) 
with (2) is the exact or near-exact equivalence of muon and electron with respect 
to weak interactions with all baryons. The nucleon pair which we assumed in the 
virtual intermediate state may equally well be replaced by some heavier hyperon 
pair (consistent with the strong interaction conservation laws) and if these particles 
contribute appreciably to the decay process, we must require their couplings to 
electron and muons to be the same if we are to leave R unchanged. In particular, 
if the interaction constants for hyperon B-decay are smaller than the nuclear B-decay 
constant, as some authors have suggested, than the coupling constants for muonic 
decay of hyperons must be smaller in the same ratio. 


* K OK 


This investigation arose from a remark by A. SIRLIN. I would also like to thank 
Prof. R. E. Prrerts for helpful comments, and H. BANERJEE for checking the 
calculations. 


29 - Il Nuovo Cimento. 
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High-Energy Pion-Nucleon Collisions and Isotopic Properties 
of the Pion-Pion Interaction. 


A. KRZYWICKI 


Institute for Nuclear Research - Warsaw 


(ricevuto il 23 Giugno 1960) 


The aim of this note is to present a short account of obtained results, con- 
cerning the high-energy (H% > 4 GeV) pion-nucleon interaction. The detailed 
explanation and justification of approximations made will be contained in a (larger) 
paper, being in publication. 


1. — Our essential assumptions are as follows. We adopt the idea of strong 
pion-pion interaction. We limit ourselves to peripheral pion-nucleon collisions, 
that is the quantities whose behaviour will be discussed correspond to certain 
selected events. We identify peripheral collisions with such collisions where only 
a small amount of momentum of one incident particle is transferred to the other 
by a single intermediate pion (1). At the same time we assume that by selecting 
experimentally observed stars according to properties which are typical for periph- 
eral collisions, we really get a sample of stars in which the contamination of stars 
due to non-peripheral collisions is small. Concerning the interaction mechanism we as- 
sume that the incident pion interacts with one and only one pion of the virtual cloud 
surrounding the target nucleon, but this interaction can lead to the excitation of the 
whole nucleon. Moreover we assume that only the pion-pion interaction and the decay 
of the produced isobar are responsible for the production of additional pions in periph- 
eral pion-nucleon collisions. In calculations we use Dremin and Cerniavski’s method 
and assume that the peripheral pion-nucleon collision is characterized by the scat- 
tering of the incident and virtual pions in a low relative angular momentum state 
and that the pion-pion scattering amplitude weakly depends on energy. The upper 
limit of the energy to which our calculations are applicable is fixed by the assumption 
that we can neglect the probability of production of more than two additional 
pions in a peripheral pion-nucleon collision. The following notation is used (we 


(1) I, M. DREMIN and D. S. OERNIAVSKI: Zurn. Eksp. Teor. Fiz., 38, 229 (1960). 


2650 


tomy 


HIGH-ENERGY PION-NUCLEON COLLISIONS AND ISOTOPIC PROPERTIES ETC. 443 


work in cm): 


P> Po momentum and energy of the incident proton; 
; o momentum and energy of the incident pion; 
BAP. momentum and energy of the isobar; , 

k, w momentum and energy of the virtual pion; 

E, ty momentum and energy of the final nucleon; 
k,;, w, momenta and energies of final pions; 


M,m, isobar, nucleon and pion masses. 
2. — First we examine the problem, whether one can obtain the branching ratio 


(n+p >n+r +p) 


k= ; 
(rap mers En) 


from the experimental branching ratio 


Re) (m + p TT) 
ay == . 


(r-+p-=r +rt+n+?) 


Let us consider an experiment which consists of counting the relative frequency 
of occurrence of two kinds of stars (in the lab. system), first with two energetic prongs 
and second with one energetic and one relatively slow prong identifiable as a proton. 
The following restrictions are imposed on the momenta (in em) of final charged 
particles, in order to ensure the peripheral character of selected events. 


[K.]> (kg), for final charged mesons, 


1 
©) |é| > h(E) , for final protons, 


h; are chosen in such a way that (1) excludes k;-q <0, t-p<0. (1) must also 
imply that the charged meson did not result from the decay of the isobar. hh, is such 
that we do not take into account protons which could be produced by the isobar 
and emitted backward in the isobar rest frame. We assume here that |p — P|/w is 
limited and its maximum value is of the order of magnitude of unity. The simplest 
choice of h, and hg is: 


fr — const, for k;-q> 0, 
(2) h;(k;) = ( 
| eee | A for k;:q (Ne 


where r is larger than the maximal momentum of a pion flying forward (in em) 
and emitted by the decaying isobar flying backward (in em). 
For events selected in such a way, one can calculate R%/R. One obtains: 


(3) R&R = (1+ )/(1+ @) 
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where 
EJ (ooo: pa? > pa CP + wR) FP) 
1° 2r8(g*[4n) [ast [(m? + (pt)) fot — py? + 4809) LO 
{tlh 
(4) 
— il far ([po® IL. (OT ne?) \/[(P — p) + Fa) UP) 
27 An°(g°/4r) [dot ((m2 + (pt)) tol (t — p)? + 4°) F(t) i 


(the restrictions imposed on the integration with respect to P are not written expli- 
citly here) F, and F, denote certain kinematical factors corresponding respectively to 
cases when in the final state we have proton or neutron. Æ’s are defined by the 


expression: 


fale. CE: 
(5) F(x) | È il 3 (Pare D he), | ky, |> hy, 


oO; 


RER 


(in the case when k, denotes momentum of a neutral pion h,= 0). We have also (1) 


(6) POL qn o(@) = EME — m? — p?)? — Am? o( M?) . 


Tu 


whe 


A very rough calculation gives (we assume that only a 3 isobar is produced): 


È 9 9 pe bua: 9.3-u-m 
(7) = Cy RE (oi a LAND 3.7 AAP ne = a ae 
8 2 (9°/470)(AP)max Pot Do 


(AP)max denotes the maximal momentum transfer in a peripheral collision, is 
of the order of # and plays here the role of a parameter. One can verify that for 
Walker’s experiment (?), at 4.5 GeV, c, is of the order of 10, so that R/R is 
not far from 4. 

Using the approximate formula (A denotes pion-pion scattering amplitude) 


(8) aa 


2 |A(ntn- > Tt) |? 


we calculated R for seven extreme cases (the index at A refers to the isotopic spin 
of the two-pion system): 


| Aso | > | A | , (PAR | : | AM | > | A | ; PACO | : | A | > | A | x | AM | ; 
| AO | 28; | AW | > |A® | : | AM | mS | A | > ED | ; 
| ACD Less | A®) | > [AO | : |A® | Co | AW | =; | A®) | 3 


Table I gives R. The first line of the table shows isotopic spins which characterise 
«large» A™’s, e.g. to the case [AM | ~|A®|>> | 4] corresponds the symbol (182) 
in the first line. ; 


(*) W. D. WALKER: Phys. Rev., 108, 872 (1957). 
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(ABE el: 
i 
if | (0) (1) (2) (OLIS) | (1, 2) (0, 1, 2) 


R | 0 0.5 4.5 0.18 | 0.5 eee 25 0.5 


T 


3. — Consider the process 7-+p >n-+n-+x*+p occuring in a bubble chamber. 
We take into account only 4-prong stars where three particles are very energetic 
and one relatively slow (in the lab. system) where the momentum is balanced and the 
slow particle is identifiable as a pion. This slow pion was emitted backward in em. 
We assume that at least two fast particles are collimated forward in em (periph- 
erality condition). Relative frequency of occurence of stars with slow pion being 
positive or negative can be measured (we denote this quantity by R, and assume 
that only one 3, $ isobar was produced) 


9A (an) |e 


‘| A(nort > nt) |? È 


(9) RE 


Table II shows R, in seven extreme cases, when the dependence of A‘ on other 
variables than iso-spin is neglected. 


(UAB IE lls 
| | | 
T | (oye PPR A) | (2) 1) | (02) (12 ec(On te) 
a 


The cross section for the peripheral process 7-+p>n+7x+©* Lp when n° is 
emitted backward is 


A2 
p° 9.3-u-m 
(10) Operiph. (TP 27 xtp) = const —— |3 MAP) 
T— back. Poo Po — Wo 


One can show that cpenpn. (mp >2r-r*p) corresponding to E — 4.5 GeV is of the 
T— back, 


order of, or somewhat smaller than 6,8 (TP 7 7) 80 that it is probably of 
the order of 1 mb. 


ok * 


I should like to offer my sincere thanks to prof. Z. KoBa for advice. 
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Die elastische und inelastische Streuung von negativen 
K-Mesonen an freien Protonen bei Energien bis 90 MeV. 


H. GOING 


Max-Planck-Institut fiir Physik und Astrophysik - München 


Nuovo Cimento, 16, 848 (1960) 


ADDENDUM 


In der genannten Arbeit wurde erwähnt, daB 55% der wom Verfasser analy- 
sierten Ereignisse von den Instituten in Bologna, Paris un Parma beigesteuert 
worden waren; es war jedoch versehentlich der Hinweis darauf unterlassen worden, 
daB die analysierten Ereignisse zum groften Teil mit denen identisch sind, die auch 
in der Gemeinschaftsarbeit von ALLES et al. (') behandelt wurden. Der Verfasser 
dankt seinen Kollegen in Bologna, Paris, Parma und München für ihre Hilfe und 
Mitarbeit. 


(1) W. ALLES, N. N. Biswas, M. CECCARELLI, R. GESSAROLI, G. QUARENI, H. Gd6oING, 
K. GoTTSTEIN, W. PUSCHEL, J. TIETGE, G. T. Zorn, J. CRUSSARD, J. HENNESSY, G. DascoLa und 
S. Mora: Nuovo Cimento, 11, 771 (1959). 
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LIBRI RICEVUTI E RECENSIONI 


M. R. WEHR and J. A. RICHARDS jr. 
— Physies of the atom. Addison- 
Wesley Publishing Company, Inc. 
(Reading, Mass., and London), 
1960, pp. x1-420, prezzo dol- 
lari 8,50. 


Ii volume, si dice esplicitamente ma 
forse un po’ troppo genericamente nella 
prefazione, non vuole essere una ras- 
segna, nè un trattato, ma un vero e 
proprio libro di testo destinato a col- 
mare la lacuna fra la fisica classica ed 
il fronte attuale della ricerca. Esso viene 
quindi ad aggiungersi a parecchie altre 
opere pubblicate, particolarmente negli 
Stati Uniti, con simile, e non facile, 
dichiarato intendimento. 

La mole di materia trattata è vera- 
mente grande. I primi tre temi sono 
cosi enunciati dagli autori, nell’ordine: 
concezione atomica della materia, del- 
l’elettricità, della radiazione. Seguono, 
sempre traducendo letteralmente i titoli: 
modelli atomici di Rutherford e Bohr, 
relatività, raggi X, onde e corpuscoli, 
concezione atomica dei solidi, radio- 
attività naturale, reazioni nucleari e 
radioattività artificiale, energia nucleare 
ed infine raggi cosmici e particelle ele- 
mentari. Ogni capitolo è corredato da 
ampi riferimenti bibliografici e da una 
ricca serie di problemi, per metà dei 
quali è offerta la soluzione: notevole 
l’uso quasi esclusivo del sistema di 
misura MKS. 

Il piano generale dell’opera non sod- 
disfa del tutto sia per numerose diva- 
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gazioni dal fine prefisso, sia per il cri- 
terio con cui le nozioni si susseguono, 
ben spesso accavallandosi: inoltre, se é 
vero che una certa superficialita nella 
trattazione è annunciata e scontata, essa 
si spinge non di rado ad un punto tale 
da chiedersi se non sarebbe stato meglio 
escludere del tutto alcuni argomenti, 
specie di natura teorica. 

A facili critiche si prestano poi le 
fantasiose analogie, di cui il testo ab- 
bonda, troppo spesso suscettibili di in- 
durre in equivoci anche gravi un lettore 
inesperto; così, mentre il libro è senza 
dubbio elementare, non se ne può con- 
sigliare la lettura che a persone di una 
certa esperienza, come vivace riepilogo 
di nozioni già maturate. 

Piacevole lo stile, americano come 
non mai, dell’esposizione e certamente 
lodevoli così l’attenzione per i protago- 
nisti della ricerca, che rivivono spesso 
in aneddoti, come la sollecitudine dimo- 
strata per l'avvenire scientifico dei let- 
tori. Risentono di questa singolare atti- 
tudine umanistica già la prefazione in 
cui si ammonisce che le conquiste mag- 
giori Si devono a ricercatori in giovane 
età, e perfino le prime due delle otto 
appendici: in una si riepilogano in ordine 
cronologico le ‘principali scoperte nel 
campo dell’atomismo dando particolare 
rilievo tipografico ai nomi degli autori 
riportati con completezza insolita, nel- 
l’altra si elencano tutti i premi Nobel 
per la Fisica e per la Chimica. 


A. LORIA 
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S. S. PENNER — Quantitative Molec- 
ular Spectroscopy and Gas Emis- 
sivities. Addison-Wesley Publ. 
Co., Reading, Mass. Pergamon 
Press, London-Paris, 1959. Pa- 
gine xv-587. 


Trattasi di un libro particolarmente 
destinato a coloro che desiderano specia- 
lizzarsi nello studio dei moderni motori 
di propulsione a reazione, e richiede già 
una certa preparazione di fondo di mate- 
matica, di spettroscopia e di meccanica 
quantistica ed ondulatoria. Esso esamina 
in generale problemi, relativi all’energia 
radiante, che di solito non sono svolti 
nelle opere classiche di spettroscopia © 
di ottica radiante, in quanto il campo 
di temperatura particolarmente consi- 
derato è intermedio tra quello della 
spettroscopia molecolare abituale e quello 
peculiare agli studi astrofisici. 

La materia esposta è suddivisa in 
19 capitoli i quali trattano nell’ordine: 
la radiazione termica; i coefficienti e le 
intensità di Einstein; i profili spettrali; 
la radiazione da righe spettrali isolate; 
la misura. della trasparenza e dell’assor- 
bimento per righe spettrali isolate; la 
determinazione sperimentale delle inten- 
sità assolute; calcoli teorici di intensità 
assolute e relative; la misura sperimen- 


tale in spettri IR della semi-ampiezza. 
di dispersione; la misura diretta della 
larghezza delle righe; la dispersione 
ottica; calcoli teorici del potere di emis- 
sione IR a temperatura moderata da 
parte di gas in stato di equilibrio; calcoli 
teorici di assorbimento isotermico delle 
radiazioni del corpo nero; la formula- 
zione matematica per problemi di tra- 
sferimento radiante di gas non uni- 
formi; il potere di emissione di aria 
riscaldata; stima del potere di emissione 
da parte di atomi di idrogeno; la tem- 
peratura di popolazione e la temperatura 
di traslazione di radicali liberi in fiamma; 
l’effetto della radiazione sulla combu- 
stione e sulla velocità di propellenti 
solidi; stima del trasferimento di calore 
radiante e delle variazioni di tempe- 
ratura in propellenti liquidi per missili. 

La trattazione della materia, tenendo 
presente il livello medio presunto dei 
lettori interessati, è abbastanza piana e 
non difficile. Im ogni caso questo libro 
raccoglie una notevole massa di infor- 
mazioni non facilmente reperibili nei 
trattati classici e, pertanto, esso è rac- 
comandabile nella biblioteca di ogni 
spettroscopista. 

La presentazione tipografica è ottima 
per chiarezza sia del testo che delle 
figure. 

G. Minazzo 
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